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Exploratory Wind-Tunnel Investigation of 
Wings and Bodies at //=6.9 


CHARLES 


H. McLELLAN* 


Langley Aeronautical Laboratory, N.A.C.A. 


SUMMARY 


Results of an exploratory investigation of wings and bodies 
ata Mach Number of 6.9 in the Langley 11-in. hypersonic tu inel 
are presented data obtained on a square and a 
triangular plan-form wing having 5 per cent thick diamond 
and on a semicircular 


Force were 
sections, on a circular cone-cylinder body, 


cone-cylinder body Pressure-distribution data were also 
obtained on the square plan-form wing and the circular body 
The results have been compared with calculations made by 
various theoretical methods. The 
ments gave overall force values in good agreeme1it with experi 
Measured pressure distributions departed appre 
of the in- 


more exact theoretical treat 


mental results 
ciably from calculated results in some 
fluence of boundary layer. 


cases because 


INTRODUCTION 


PERFORMANCE of contemplated 
Mach 


y EVALUATING THE 
high 
Numbers, it has been necessary up to this time to use 
theoretical results without experimental verification. 
An 11-in. hypersonic wind tunnel capable of providing 


missiles at supersonic (hypersonic) 


the necessary experimental data has recently been 
developed and put into operation by the N.A.C.A. at 
present 


its Langley laboratory. The purpose of the 


paper is to discuss the results of an exploratory in- 
vestigation of wings and bodies in this new facility at a 
Mach Number of 6.9 and to compare these results with 


available theories. 


SYMBOLS 


= wing chord 
wing chordwise force coefficient 
wing or body drag coefficient based on plan-form area 
wing or body lift coefficient based on plan-form area 


C.P. = center of pressure 


Presented at the Hypersonic Aerodynamics Session, Nineteenth 
Annual Meeting, I.A.S., New York, January 29-February 1, 
1951. 

* Head, Tunnel Section. 


11-In. Hypersonic 


D = drag 
L = lift 
= free-stream static pressure 
= surface pressure 
= angle of attack 
= radial position measured from the top of the body 


DESCRIPTION OF APPARATUS AND MODELS 


Tunnel 


hypersonic tunnel is shown 
schematically in Fig. 1. Air is discharged from a 50- 
atmosphere high-pressure tank through a_ pressure 
regulator to a heater, where it is heated to about 700°F. 
After leaving the heater, the air is expanded through a 
two-dimensional single-step nozzle designed for a Mach 
Number of 7.08 and then passes through an adjustable 
second throat to a vacuum tank. The tunnel is de- 
scribed more fully in reference 1, and the results of 
surveys of the flow through the nozzle are presented in 
reference 2. The flow surveys showed that the Mach 
Number at the test section was approximately 6.9 and 
that a central core of the stream having a cross section 
had sufficiently uniform flow for 


The Langley 11-in. 


about 5 in. 
model testing. 


square 


The high stagnation temperature used in the present 
tests is required to maintain the air temperature in the 
test section above the liquefaction temperature of air 
at the test-section pressure. The presence of liquid air 
particles in the nozzle has been observed when a suffi- 
high stagnation temperature used. 
condition, pressure 


ciently was not 


Furthermore, measure 


for this 
flow have shown that large 
of liquefaction. 


ments in the increases in 
entropy occur which are indicative 
It is therefore believed essential that the high stagna 
tion temperature be used. This problem of liquefac 

2 


tion is discussed in greater detail in reference 3. 








































































642 FOURNAL OF THE ABRONAUTICAL 
50-ATMOS __ VACUUM TANK 
PRESSURE ; 
TANK QUICK 
OPENING 
VALVE 
TEST 
= SECTION 
ihe + Woon HP 
/ “HEATER ADJUSTABLE \ COOLER 
| PRESSURE pn 
REGULATOR THROAT 
Fic. 1. Schematic diagram of the N.A.C.A. Langley 11-in. 
hypersonic tunnel. 
Le g" ~ 
- I 
4 uro(_ ) 
t 
(a) CIRCULAR CONE CYLINDER BODY 
a 8" — 165d- 
—e | /_ 
(b) SEMI-CIRCULAR CONE CYLINDER BODY 
i 2 + ne 
eS I i 
4" 4" 
=a es ! ! 
MACH ANGLE te 1.173" 
= k i 


(C) SQUARE WING WITH 5% (d) TRIANGULAR WING WITH 5% 
DIAMOND AIRFOIL SECTION DIAMOND AIRFOIL SECTION 


Fic. 2. Wing and body configurations tested 
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Fic. 3. Square plan-form force and pressure models. 
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Models 


The four models used in this investigation are showy 
in Fig. 2. The two body shapes consisted of a circular 
cone-cylinder body having a 10° half-angle circyla; 
cone with a cylindrical afterbody four diameters jy 
length and a semicircular cone-cylinder body having 
the lower surface flat. The cross-sectional areas of the 
two bodies were equal at corresponding streamwig 
stations; thus the cone angle was increased for the flat. 
bottomed body. 
had a square plan form, and the other had a triangular 
plan form with 60° swept leading edges. 
had 5 per cent thick diamond airfoil sections. The test 


Two wings were also included. One 
Both wings 


Reynolds Number was approximately 1,000,000 based 
on the 4-in. chord of the square plan-form wing 
This Reynolds Number corresponds to a wing of 4-ft 
chord flying at a Mach Number of 7 at an altitude of 
about 120,000 ft. Fig. 2 also shows the small portions 
of the wings which lie within the regions of three- 
dimensional flow from the tips of the square plan-form 
wing and from the apex of the delta wing. 


Instrumentation 


Both pressure and force measurements were made 
on the circular body and the square plan-form wing, 
while only force measurements were made on the semi- 
circular body and the triangular plan-form wing. 
The forces were measured by means of a strain-gage 
balance located in the sting support of the model. 
The wing force models were attached to the balance by 
the conical support shown in Fig. 3a. Corrections to 
the lift and drag due to the aerodynamic forces on the 
unshielded portion of the conical support have been 
applied to the wing force test results. The two body 
force models were mounted on a support that was 
shielded from the stream. Two balances were used 
one for low angles of attack where the forces were 
extremely small, and one for the high angles of attack 
where the forces were relatively large. 

The pressures on the circular cone-cylinder body were 
measured by one row of orifices running the length of 
the body. Pressure measurements around the model 
at angles of attack were obtained by rotating the model 
about its longitudinal axis. 

Fig. 3b shows the chordwise row of orifices that were 
used to obtain the pressure distributions on the square 
plan-form wing. It can be seen that the area affected 
by the three-dimensional flow at the tips is small and 
that the distributions obtained by these orifices will be 
essentially two-dimensional. Since the model was 
extremely thin, it was impractical to locate the con- 


within the model. 


necting tubes entirely 
Tubes were therefore brought out through the surface 


opposite the one in which the measurements were being 


pressure 


made. Experimental studies have indicated that the 


presence of these tubes and the model support does 
flow over the 


not have a measurable effect on the 
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WIND-TUNNEL 
top of the wing. The pressures were measured by 
means of optically recording aneroid units described in 
reference 1. 


RESULTS AND DISCUSSION 


Wing Pressure Distributions 


With the aid of existing supersonic flow tables, the 
pressure distribution in the two-dimensional flow region 
of the square plan-form wing can be easily determined 
by the use of the characteristics system, and there is 
little reason to use the approximate methods. A com- 
parison of the experimental and theoretical pressure 
distribution for the square plan-form wing is shown in 
Fig. 4 for angles of attack of 0° and 10°. The pressure 
distributions are presented in the form of a local static 
pressure P2 minus the stream static pressure p, divided 
by the stream static pressure. This ratio has a value 
of zero when the local surface pressure is equal to the 
stream pressure and is —1 at zero absolute surface 
pressure. Theoretical pressures obtained by the char- 
acteristics system are presented as dotted curves. The 
results indicate large deviations of the experimental 
pressures from the theoretical values near the leading 
edge of the wing. A similar deviation occurs just back 
of the break in the surface at the maximum thickness. 
The deviation results from the rapid growth of bound- 
ary layer near the leading edge and just back of the 
maximum thickness at the high test Mach Number. 

A study of the growth of boundary layer at this 
Mach Number showed that, at the test Reynolds 
Number of 1,000,000, the calculated pressure rise near 
the leading edge due to a laminar boundary layer was 
in good agreement with the experimental deviation of 
the pressures from theory. At a Mach Number of 7, 
the boundary-layer displacement thickness is about 
ten times as great as at a Mach Number of | at the same 
Reynolds Number. This thick boundary layer can be 
seen in the schlieren photograph of the square plan- 
form wing presented in Fig. 5. 

At 10° angle of attack, the flow over the rear surface 
of the wing is separated, and the pressures remain 
constant at the same value as for the forward upper 
surface. It is evident from the pressure diagram (Fig. 
4) that this flow separation has only a small effect on 
the lift because of the small lift contribution of the 


upper surface at high Mach Numbers. 


Wing Force Characteristics 


There is no available theoretical method for rigorous 
calculation of the characteristics of the complete wings. 
The linearized theory, however, may be used for an 
approximate evaluation of the tip effects. In the case 
of the square wing, the tip effect reduces the lift-curve 
Analysis of the pressure data 


slope by 7 per cent. 
(Fig. 4) indicates that the section lift-curve slope is 
reduced by about the same amount because of viscous 


effects. Thus, a total reduction in slope of about 15 
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per cent below the two-dimensional values obtained 
by the characteristics method would be anticipated. 
The test data (Fig. 6) indicated a reduction of about 10 
It would appear from this result that the tip 
linearized theory. 


per cent. 
effect is overestimated by the 
Further investigation is required to evaluate the tip 
effect more accurately. In Fig. 6, the lift curve, calcu- 
lated by Linnell’s hypersonic approximation for two- 
dimensional flow‘ is essentially the same as for the 
characteristics theory at low and moderate angles of 
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Fic. 4. Typical pressure distributions on the square plan-form 


wing, 1J = 6.9. 





Schlieren photograph of flow about the square plan-form 


wing, a = —0.5°, MV = 6.9 


Fic. 5. 
























































644 JOURNAL OF THE AERONAUTICAL 
A o O EXPERIMENTAL VALUES 
48> 60° f 
7 oS o *Y 4 
C.R, TWO-DIMENSIONAL 
324 CHARACTERISTICS L C.P. 
THEORY . 
CL 4 CL, LINNELL’S — SY D -2 
APPROXIMATION 
164 - CL, TWO-DIMENSIONAL a 
F B CHARACTERISTICS 
THEORY 
4 4 / J . 
“C,, LINEARIZED THEORY 0 
T T T T T T ' 
O 8 16 20 
a, DEG 
Fic, 6. Lift of the square and the delta plan-form wings, 
M = 6.9. 
» EXPERIMENTAL VALUES 
_ (| | | 
20 
| Cp, TWO-DIMENSIONAL / 
i6- CHARACTERISTICS THEORY 
Cp, TWO-DIMENSIONAL 
Cn. LINNELL’S CHARACTERISTICS THEORY 
124 0» APPROXIMATION 03-\ WITH Gf =.0028 
Cp Cp, TWO-DIMENSIONAL 
CHARACTERISTICS THEORY 
084 02 
o Af 
044 o1 hia 
= AC 
men wean 
O 4 8 2 6 20 24 O 4 8 
a, DEG a, DEG 
Fic. 7. 


Drag of the square and the delta plan-form wings at 
M = 6.9. 


attack, while at high angles, which are actually beyond 
the range in which it should be applied, it is appreciably 
greater. 

In the case of the triangular wing, only the two- 
dimensional flow area ahead of the Mach line from the 
apex can be analyzed by exact theory. However, with 
the 60° swept leading edges and the high test Mach 
Number, this area constitutes approximately three- 
quarters of the wing area. The lift of this area is 
slightly greater than for a two-dimensional wing. 
The linearized theory, which can be applied at low 
angles of attack, shows that the total lift of the tri- 
angular wing is the same as the two-dimensional lift. 
The lift has therefore been compared over the whole 
range of angle of attack with the exact two-dimensional 
theory in Fig. 6. It can be seen that the experimental 
slope is approximately the same as predicted by the 
two-dimensional theory and slightly greater than that 
for the square plan-form wing. The difference between 
the two-dimensional theoretical results and the test 
data from the delta wing is probably of the same mag- 
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nitude as the experimental errors because of the ac- 
curacy of the measurements and the errors in evaluating 
the effect of the support sting on the forces on the delta 
wing. 

The center of pressure as obtained from pressures 
along the centerline of the square wing is also presented 
in Fig. 6. The measurements indicate that the center 
of pressure is slightly more forward at low angles of 
attack than theory predicts and tends to become equal 
to the theoretical value at high angles. 

The drag characteristics from the force measurements 
for the two wings are presented in Fig. 7. The solid 
curve shown in this figure is the drag coefficient calcu- 
lated by the two-dimensional characteristic theory, 
The linearized theory gives essentially the same drag 
coefficients for the two wings. At zero angle of attack, 
the linearized theory gives drags that are nearly equal 
to the two-dimensional theoretical drag. However, as 
the angle is increased, the rate of increase of the drag 
coefficient from the linearized theory is appreciably 
less than for the characteristic theory, since the lift is 
underestimated by the linearized theory. 

Linnell’s hypersonic approximation gives essentially 
the same drag as the characteristics theory at low and 
moderate angles and only slightly higher drags at high 
angles of attack, although they are beyond the angles 
at which this theory would be expected to apply. 

At low angles, the drag of the square plan-form wing 
measured by the balances is appreciably greater than 
the inviscid theory predicts. However, when a fric- 
tion drag coefficient of 0.0028 calculated for laminar 
flow at the test Reynolds Number is added to the in- 
viscid theoretical drag, the total agrees with the force 
measurements. The value of the friction drag coeffi- 
cient is twice the pressure drag coefficient at zero angle 
of attack. Of course, the 5 per cent diamond section 
used in this investigation has a relatively low pressure 
drag. Only moderate percentage reductions in total 
drag could therefore be obtained for this section by 
reducing the wing thickness below 5 per cent at the test 
Reynolds Number. From the same _ standpoint, 
changing to a wing section that might be expected to 





have lower pressure drag, such as a wedge with the | 


point of maximum thickness at the trailing edge, will 
have only a secondary effect on the total drag. If it ts 
assumed that the boundary layer remained laminar, 
the effect of the thickness and wing section would be 
relatively greater at higher Reynolds Number. How- 
ever, the test Reynolds Number is within the range 
of interest in missile applications. 

Fig. 7 shows that the drag of the delta wing is slightly 
higher than that of the square wing at high angles 
of attack, mainly because of the higher lift. At low 
angles, the drag is about 30 per cent higher than for the 
square wing. The linearized theory indicates that the 
pressure drag of the delta wing is only about 2 per cent 
greater than for the square wing. The increase is pri- 
marily a result of a decrease in the effective Reynolds 
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Number of the delta wing. Based on the local chord, 
the delta-wing Reynolds Number is about 90 per cent 
of that of the square wing at the centerline and de- 
creases to zero at the tip. 

Lift-drag ratio data for the wings are presented in 
Fig. 8 as a function of angle of attack. The theoretical 
curve obtained using the inviscid two-dimensional 
characteristic theory is compared with the values 
given by the diamond symbols that were obtained from 
an integration of the pressures over the centerline of the 
square wing. It can be seen that the pressure results 
are in reasonable agreement with the inviscid theory 
and indicate a maximum L/D of over 10. 

When friction is included, the theoretical curve of 
L/D drops appreciably at low angles, and the maximuin 
L/D is about 6, which is about 25 per cent less than that 
expected at a Mach Number of 2. Since considerable 
scatter exists in the force data for the square plan-form 
wing (the square symbols), it appears that the theoret- 
ical curve with friction can be considered to represent 
approximately a faired curve of the data. Considerably 
less scatter is evident in the force data from the delta 
wing (the triangular symbols) because of improved 
instrumentation and techniques. The L/D of the delta 
wing is slightly lower than for the square wing, pri- 
marily as a result of the lower effective Reynolds 
Number of the delta wing. At high angles, there is no 
noticeable difference in lift-drag ratio of the two wings; 
in fact, they closely approach the calculated curve 
corresponding to a flat plate with no chordwise com- 
ponent. This curve is independent of Mach Number 
and represents an upper limit for any airfoil regardless 
of the plan form, provided no leading-edge suction is 
obtained. At the higher angles of attack, the pressure 
L/D values are close to the limiting curve. 

Because of the large frictional component of the drag, 
it appears that an appreciable increase in maximum L/D 
could be effected for a wing of given area and thickness 
ratio by reducing the aspect ratio, thereby increasing 
the Reynolds Number. At Mach Numbers of the order 
of 7 or higher, it appears from this investigation that 
aspect ratios considerably lower than unity can be 
used without appreciable loss of lift-curve slope. 

In order to illustrate the nature of the flow about the 
wing, a schlieren photograph of the square plan-form 
diamond airfoil wing mounted on the sting balance at 
an angle of attack of —16.8° is shown in Fig. 9. In 
this photograph, it can be seen that the shock is close 
to the high-pressure surface of the model. On the low- 
pressure side of the model, an expansion to a Mach 
Number of nearly 11 just back of the leading edge 
would be expected in the absence of viscous effects; 
however, a rapid growth of boundary layer results in 
the occurrence of a shock at the leading edge. A 
maximum Mach Number of about 8 occurs well back 
of the leading edge. Separation occurs at about the 
maximum thickness on the low-pressure side at this 
angle of attack. 
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Cone Pressure Distributions 


One of the simpler cases of bodies of revolution to 
analyze is the cone-cylinder body at zero angle of attack. 
In this case, the pressures on the forward or conical 
portions of the body can be obtained from the work of 
Taylor and Maccoll,® which has been calculated and 
tabulated by the M.I.T. group in reference 6. The 
flow over the rear or cylindrical portion can be calcu- 
lated using the flow from the conical portion and the 
three-dimensional For the 
particular body used in this investigation, it is not 


characteristics system. 
possible to apply the linearized theory because the 
half-cone angle of 10° is greater than the Mach angle. 
A comparison of the theoretical and the experimental 


results is presented in Fig. 10. On the conical portion 
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and on the rear of the cylindrical portion of the body, 
the experiment and theory agree, while on the cylin- 
drical surface near the cone-cylinder juncture, the 
pressures are considerably higher than predicted by 
theory. This is believed to be largely due to the 
boundary layer altering the effective body shape at the 
juncture. 

At angles of attack, several approximate theories are 
available for analyzing the flow about the conical por- 
tion of the body. The radial distribution of the pres- 
sure ratio (Pp. — p,)/p; around the cone measured in 
degrees from the top of the cone is compared with the 
available theories in Fig. 11 for an angle of attack of 
6.7°. Experimental points at each radial station are 
presented for stations. The 
points nearest the leading edge should probably be 
ignored because they are affected by local conditions 
such as accuracy of the cone at the point and rapid 
boundary-layer growth near the leading edge. A first- 
order theory, shown as the solid line, was developed by 
Stone’ and was tabulated by an M.I.T. group under 


several longitudinal 
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the supervision of Kopal.* This work, which will p¢ 
referred to as the Stone-Kopal theory, is a first-order 
theory that assumes that the entropy is constant in any 
meridian plane but varies with position around the cone. 
In addition to the usual angle-of-attack limitations oj 
a first-order theory, the Stone-Kopal theory is limited 
to angles of attack less than the half-cone angle, singe 
the coordinate system is based on the wind axis, 4 
more recent theory has been developed by Ferri! which 
uses a coordinate system about the body axis instead 
of the wind axis and which introduces the concept of q 
vortical layer next to the cone surface. The entropy 
is shown by Ferri to be constant on the cone surface 
except for a singular point at the top of the cone. [py 
the work of Stone, the entropy was allowed to vary 
around the surface. The change in the axis of the 
coordinate system allows for solutions of the equations 
at angles of attack greater than the surface cone 
angle, although the theory is not expected to be valid 
at high angles of attack. <A third theory, presented in 
Fig. 11, was developed by Grimminger'’ and is based 
upon the Newtonian corpuscular concept in which the 
gas stream maintains its speed and direction unchanged 
until it strikes the solid surface, whereupon it loses the 
This 
approximation of the flow is valid at high supersonic 
Mach Numbers, particularly when the angle of attack 
However, this method does not predict 
Re- 
finements of this theory, which include the centrifugal 


component of momentum normal to the surface. 


is appreciable. 
the pressures on surfaces shielded from the flow. 


forces of the air as it flows about the body, were also 
presented by Grimminger. The case included in Fig. 
11 assumes the velocity over the surface to be the com 
ponent of the free-stream velocity parallel to the sur 
face. 

At 6.7° angle of attack, the results from Stone's 
work give approximately the same values as experiment 
at the sides of the models but appreciably lower values 
at the top and bottom. The Ferri values are in con- 
siderably better agreement; however, they tend to be 
slightly high near the top of the cone. Grimminger's 
results, without allowance for centrifugal force, agree 
reasonably well with the experimental data on the high- 
pressure portion of the cone but are slightly low on the 
low-pressure or upper surface. The inclusion of the 
centrifugal force in Grimminger’s theory increases the 
discrepancy at the test Mach Number. Grimminger" 
also showed that the method used overemphasized the 
effect of the centrifugal forces and suggested several 
methods for reducing the calculated centrifugal forces. 
These reduced values would have given better agree- 
ment but not so good as without the centrifugal force 
corrections. 

At an angle of attack of 14°, the results of Stone’s 
work can no longer be applied. Although this angle of 
attack is beyond the range for which first-order theories 
in generalsuccessfully approximate the physical phenom- 
enon, nevertheless, the application of Ferri’s theory 
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WIND-TUNNEL 


Ids good agreement over the high-pressure portion 


yie 
of the cone but gives slightly higher pressures than 
However, 


experiment over the low-pressure surface. 
it should be pointed out that, over the top of the cone, 
the pressures are probably affected to a considerable 
extent by boundary-layer flow about the cone. Grim- 
minger’s hypersonic approximation without centrifugal 
force agrees extremely well with the experiment; how- 
ever, it cannot be applied on the upper 28 per cent of 
the model. The inclusion of centrifugal force again 
increases the discrepancy. 


Body Force Characteristics 


Of the various approximate theoretical methods 
used in the analysis of the flow about the cone at angle 
of attack, only Grimminger’s method can be readily 
applied to the flow about a complete body of circular 
cross section. The work of Ivey'! can only be applied 
on the cylindrical portion of the body and must be 
used in conjunction with other theories for the conical 
portion of the body. The approximate method of 
Grimminger applies to both the cylinder and the cone. 

Lift coefficients based on the projected area of the 
circular cone-cylinder body are presented in Fig. 15 
as a function of angle of attack. The theoretical values 
based on Grimminger’s approximation are included for 
the circular cone-cylinder body. To permit a compari- 
son of the lift of the bodies with that of the wings, the 
theoretical curve of lift versus angle of attack for the 
5 per cent diamond profile has been included. The lift 
predicted by the hypersonic theory is slightly higher 
than the experimental values for the circular body. It 
is interesting to note that the circular body has a lift- 
curve slope of about 60 per cent as great as that for a 
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Fic. 12. Pressure distribution around the conical portion of 
the circular cone-cylinder body at an angle of attack of 14°, 
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Lift data for the flat-bottomed body with semicir- 
cular cross section are also included in Fig. 13. As 
shown by Sanger,'® bodies of this type are of interest in 
hypersonic applications because they have the possibili- 
ties of increased lift-curve slope, increased lift at 
maximum L/D, and higher maximum L/D than com- 
parable round bodies. It can be seen that the body 
having the semicircular cross section has a lift-curve 
slope considerably greater than the circular cross-sec- 
tion body; in fact, its slope is about the same as that 
for a two-dimensional wing. It should also be noted in 
Fig. 13 that, as expected, the slope of the lift curve at 
negative angles (the curved portion on the high- 
pressure side) approaches that of the circular body more 
nearly than that of the wing. 

The drag coefficients based on the plan-form area for 
the two bodies are plotted against angle of attack in 
Fig. 14. Comparison of the results of force and pres- 
sure measurements on the circular body shows that, 
for this relatively low fineness ratio circular body, only 
15 per cent of the drag at zero angle of attack is due to 
friction. This compares favorably with the calculated 
friction drag assuming the flow is entirely laminar. 
The theoretical values of total drag for the circular 
body are slightly higher over most of the range than 
experimental values, even though friction is not in- 
cluded in the theory. At high angles, some of the loss 
in drag can be accounted for by the lift being lower 
than predicted. The differences in drag coefficients of 
the two bodies based on the projected area in this 
figure are small; however, if the drag coefficients are 
based on maximum cross-sectional area, the drag co- 
efficient of the semicircular body would be about 40 
per cent larger than for the circular body. The higher 
drag is largely due to the blunter-cone angle of the 
semicircular body. 

The lift-drag ratios of the bodies are presented in 
Fig. 15. A curve of L/D based on the theory of 
Grimminger is included for the circular bodies. The 
experimental lift-drag ratios of the circular body are 
only slightly lower than those in this nonviscous theory. 
This is largely the result of the fact that the pressure 
drag constitutes the major part of the drag because of 
the low fineness ratio of this body. 

Although early data indicated that the maximum 
lift-drag ratios of the two bodies were essentially the 
same, the more recent test results of Fig. 15, in which 
the experimental scatter has been greatly reduced by 
the use of improved instrumentation and techniques, 
show that the maximum lift-drag ratio of the flat-bot- 
tomed body is about 15 per cent greater than for the 
circular body. This trend was anticipated from the 
work of Sanger.'?» The maximum lift-drag ratio also 
occurs at a higher angle of attack and lift coefficient for 
the flat-bottomed body than for the circular body. 

At high angles of attack, the L, D values approach 
the values for a flat plate, indicating that at these high 


angles of attack the bodies will give about the same 
L/D values as any other good lifting surface. The 
theoretical L/D curve for the 5 per cent diamond wing 
section with friction is also included for comparative 
purposes. It can be seen that, at the test Reynolds 
Numbers, the maximum L/D of the bodies is about 
half that of the wing section. The L/D ratios of these 
bodies will not be affected to nearly the same extent by 
Reynolds Number as the wings, since the friction is 4 
considerably smaller portion of the overall drag. 


CONCLUDING REMARKS 


From the results of this exploratory investigation, it 
can be concluded that the aerodynamic characteristics 
of wings and bodies at a Mach Number of the order of 7 
can be predicted by available theoretical methods with 
the same order of accuracy usually obtainable at lower 
speeds, at least for cases in which the boundary layer is 
laminar. Allowance should be made for the effects of 
the relatively thick boundary layers, which occur at the 
high test Mach Number, in altering the effective shape 
of the wing or body, especially near the wing leading 


edge. 
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The Evaluation of Theoretical Critical 
Compression in Sandwich Plates’ 


J. N. GOODIER?# anv I. M. NEOU?t 
Stanford University 


ABSTRACT 


Formulas for the critical compression of sandwich plates, 
which take into account everything that might influence the 
When simplicity is 
achieved by concessions in accuracy, the range of validity is not 


buckling load, are extremely complex. 


easily distinguished. The evaluations presented here use check 
points obtained from an ‘‘exact’’ analysis to test the accuracy 


of simpler formulas, with the result that one of the latter (pro- 
posed in reference 5) is found to be remarkably accurate over a 
wide range of geometrical and elastic parameters. 


(1) INTRODUCTION AND SUMMARY OF RESULTS 


b lise COMPRESSIVE STRENGTH of thin sheets can be 
realized only if they are stabilized against buck- 
ling. In sandwich construction, two such sheets are 
bonded to a core slab of different (light) material. The 
core is sometimes isotropic, sometimes anisotropic. 
Only the isotropic case is considered here. 

Panels of this construction give rise to a set of prob- 
lems of strength, stiffness, and stability analogous to, 
but by means identical with, the well-known 
problems of the ordinary homogeneous and isotropic 
One of these is “cylindrical buckling,”’ 
The panel is so wide that lines 
In the ordinary 


no 


elastic plate. 
illustrated in Fig. 1. 
such as EF can be taken as uncurved. 
plate, a unit width is treated as a Euler column, but, 
since anticlastic curvature is prevented, the Young's 
modulus FE is replaced by E/(1 — v?), v being Poisson’s 
ratio. ** 

If it could be assumed that plane sections remained 
plane, the corresponding sandwich problem would 
require merely a straightforward modification to take 
account of the fact that the section is composite rather 
than homogeneous.¢+ This assumption will be valid 
when the modulus of the core (Z,) and the modulus of 
the face sheets (Z,) are not too far apart. The sand- 
wich constructions of interest to aeronautical engineers 
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** Actually a wide strip must begin to bend with anticlastic 
curvature, but this curvature does not develop as bending 
Proceeds. There is a transition from E to E/(1 — The 
critical load for infinitesimal buckling would appear, therefore, 
not to require the factor 1/(1 — v?). 

tt See Eq. (11). 
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have light cores with low moduli—not so low that they 
leave the sheets practically free to buckle as inde 
pendent sheets and not so high that plane sections 
remain plane. The actual buckling is either like 
column buckling (quasi-Euler buckling) or a_ short 
wave ‘wrinkling’ of the face sheets. In the former, 
the core may exhibit a substantial shearing deforma- 
tion; in the latter, it acts like an elastic foundation, 
and the buckling deformation may be confined to the 
layers adjacent to the face sheets. These statements 
rest on the numerous investigations published during 
the past 10 years. (A few of these are cited in the 
list of references.) 

One group of these deals only with the quasi-Euler 
buckling, the procedure simplest in principle being 
‘correction due to shear’ in 


the introduction of the 
the elementary bending theory. 
with both quasi-Euler buckling and with wrinkling, 
treating the core deformation by the methods of the 


Another group deals 


ordinary theory of elasticity. This involves the as 
sumption that the (usually small) compressive load 
borne by the core does not influence its response to 
lateral forces. A third group goes further by taking 
this influence into account. The theory used for the 
core is an extension of the ordinary theory of elasticity 
developed in the past 20 years which can deal, in 
general, with elastic the 
of initial stress severe enough to exert a destabilizing 
effect; a simple and well-known instance is the lateral 
loading of an initially compressed bar. 

The solutions of the sandwich problem so obtained 
Prior to their appear 


deformations in presence 


are necessarily complicated. 








Fic. 1. Cylindrical buckling of a wide plate in compression. 
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Fic. 3. Reactions between core and faces induced on buckling. 


ance, solutions of an approximate character were pro- 
posed, having the great advantage of simplicity. 
However, as with all except the extremely complicated 
exact theories, the effects of the approximating and 
simplifying assumptions on the range of validity of the 
results are not evident, even when attention is entirely 
confined to the idealized problem of the geometrically 
and elastically perfect sandwich. 

They can be made evident, attention being so con- 
The 
purpose of this paper is to put forward a conclusion, 
evaluations of 


fined, by comparison with the exact results. 
reached after extensive numerical 
several approximate solutions and one exact solution 
of the problem. It is that one of the simpler theories, 
which introduces a simple correction factor to account 
for the destabilizing effect of the compressive load in 
the core, is accurate over an extremely wide range of 
sandwich variables—from extremely thin sheets and 
extremely deformable cores to thick sheets and cores 
with moduli comparable to those of the sheets. This 
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range of course includes, but goes far beyond, that of 
interest in aeronautical construction. 

THE CORRECTED PLANE-STRAIN THEORY 


(2) 
The simpler theory referred to uses plain-strain equa- 
tions for the core. It is a modification of the theory of 
Williams, Leggett, and Hopkins,*® which in turn is a 
modification of that of Gough, Elam, and de Bruyne.! 
In both of these papers the core is treated by the 
equations of plane stress rather than plane strain, 
Plane stress implies anticlastic curvature in the core 
when it bends. As indicated previously in this paper, 
it may not be incorrect to admit such curvature at 
incipient buckling. However, in the interests of clarity 
and because the exact theory we use for comparison 
is in terms of plane strain, it seems preferable to restate 
concisely the simpler theory in these terms also. The 
conditions of continuity at the interface between core 
and face sheets are affected, as well as the equations 
of the core itself. 

We need first an account of the prebuckling state of 
compression. Let the whole sandwich be compressed 
vertically, as by a pair of rigid testing machine heads, 
so that a uniform vertical strain is imposed throughout 
(Fig. 2,a and b). No forces are applied horizontally— 
on AC or BD—but, if the Poisson's ratios of faces and 
core are different, the differing tendencies to horizontal 
expansion will require reconciliation. This is con- 
sidered to be brought about by some horizontal com- 
pressive stress in the faces and some tensile stress in 
the core (or vice versa) in such a ratio that the resultant 
force per unit run vertically is zero for the faces and 


core together. Using the subscripts f for the faces 


and c for the core, this means 


2tose + hoe = 0, Zor, + hore = P (1 


the notation being indicated on Fig. 2. P is the load 
per unit run. 
Equality of strains in the x- and z-directions means 


(1/ Ey) (627 — veer) = (1/ FE.) (Ore — Ve zc) (2) 


(1/Ey)(o2z7 — veors) = (1/E.) (G22 — VeOre) (0 


where / and vy mean Young’s modulus and Poisson's 
These four equations determine 
The ratio of the two vertical 


ratio, respectively. 
and Gz. 


|; — ve + 2rp(1 — »,”) 
=p = (4) 
1 — vp + 2rp(l — vy) 


Orfs Ozf, Trey 


stresses 1s 


Ore 


where p = E,/E, andr = t/h. 

All critical conditions for buckling will be expressed 
in the form of a critical value for o,;/E,. When this 
is known, o,, can be found from Eq. (4), and the 
critical load can be found from Eqs. (1). In_ the 
numerical evaluations, we have not attempted to cover 
a range of values of the Poisson’s ratios but have taken 
vy, = '/3, and », = Oat once. 
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CRITICAL COMPRESSION 
When buckling occurs, it is assumed that the faces 
nd into sinusoidal form with one or more half waves 


be 
For convenience, 


in the length (Fig. 2, c, d, and e). 
we suppose at first that the plate is infinitely long 
rtically and that an arbitrary wave length is im- 
posed. We then inquire what critical stress is required 
for its maintenance. In an actual panel of finite 
tength, we have to identify the half wave length with 
2, or with a fraction 1/n of 


a9 


ve 


the dimension AC, Fig. 
it, and select the integer to give the smallest critical 
stress. When this integer is 1, the buckling is of the 
quasi-Euler type. When it exceeds 1, it is usually 
large, and the buckling is described as wrinkling or 
crinkling. 

The derivation of the “plane-strain formula’ for 
the critical stress in the faces is as follows: The dif- 
ferential equation of the left-hand face sheet in Fig. 3 
is 


t dq’ 
1 =0 


d 1, 


1 
dx* 


dv, 
2 


l a ( 5) 
dx? 


D + p’ - 


2 dx 


which is Eq. (14) of reference 3; v is the deflection, 
P, is the compressive force per unit width in the face 
only, and p’ and g’ are core reactions as indicated in 
Fig. 3. 

The sinusoidal plane-strain displacements of ap- 


propriate type are 


= + [Ai + A.(1/a)(3 — 4»,)] sinh ay + 


Avy cosh ay} cosax? (6) 
v = (A, cosh ay + Ay sinh ay) sin ay 


and these are adopted as the displacements in the 


core. The constants A;, A» have to be related to 
boundary forces p’, qg’ at y = +h/2. We find 
/ E. f . 7 

p = ) [aA, + (1 — 2y,)Ao] S + Ace,C} sin ax 

1+», 
(7) 

, Ee { . ? 

q = ’ [aA + 2( l v-)A2]| C + Ave,S { COS QX 
l+»y 


cosh (ah /2) 


where S = sinh (ah/2) = sinh e, and C = 
= cosh @;,. 


The conditions of continuity between core and face 


(v)y = (h/2) = 1% | 


are 


(S) 


t dv 1 — »,? q’ (**) 
2 dx? ky at Ox y (A , | 
corresponding to those adopted and explained in 


reference 3, Eq. (43). The second equates the tensile 
strain in the x-direction, in the core and the face, at 
the junction. 

When Eqs. (6) and (7) are inserted in Eq. (5) and 
the second equation of Eqs. (8), observing the first 
equation of Eqs. (8), the result is two homogeneous 
equations in A; and A», and the condition of consist- 


ency yields the critical condition. This reduces to 
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oy e,” F+ GT + HT? 9) 
- an rrr ( 
FE, 3(1 — vf) I+ JT + KT? 
where o, denotes the critical value of o,,, and 
T = tanh (af/2) = tanh e, = tanh (e,/r) 
(en, €, are defined by e, = e,/r = ah/2) 
, ; l — pv 
F = [1 + (1 — »)rje. — 
4(1 + vy.) pe, 
; (1 — v,*)r 
G = (1 — 2»,)r+ nes 
1(1 + v,) pe," 
r l o— vy 
H = (1 — ») —e,+ 
er 4(1 + vy.) pe; 
IT = (1 + »)pe, + (1 — »,)(1 — vf*)(r/er) 
J = (1+ »,)(8 — 4y,)pr 
K = —(1 + »)pe; 
Eq. (9) is an explicit formula for the critical stress. It 


is the plane-strain analog of one given by Williams, 
Leggett, and Hopkins. When r (i.e., ¢/h), p (ie., 
E,/E.), vy, and vy, are given, Eq. (9) yields o;/Hy as a 
function of e,, which is at/2 or mt/X, ¢ being face thick- 
ness and \ wave length of buckles. 

Curves representing Eq. (9) (with vy = '/3, ». = 0) 
in this manner are shown in Figs. 4a—4g. In Fig. 4f, 
0.002. Eq. (9) is repre- 


wo curves are shown for 7 
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sented by the broken line. This prevails wherever 
two curves (one broken, one solid) appear for one value 
of r. Otherwise, the solid curve corresponds to Eq. 
(9). 

The crosses appearing on these figures are points 
calculated from the exact formula specified below. 
They serve as checks on the accuracy of Eq. (9). It 
will be observed that these check points lie on the solid 
curves, not on the broken ones. Thus, the theory just 
given, resulting in Eq. (9), is found correct only where 
the broken and solid curves coincide, so that only the 
latter appears. 

The solid curves in Figs. 4a—4g are obtained by 
applying to Eq. (9) a correcting factor proposed by 
Williams, Leggett, and Hopkins page 
12) to express approximately the destabilizing effect 
This factor, 


(reference 5, 


of the compressive stress in the 
to be applied to the right-hand side of Eq. (9), is 


core. 


1/{1 + (@re/ ory) (h/2t) | (10) 


where o;,, or, are the prebuckling stresses in core and 
faces and their ratio is given by Eq. (4) independently 
of the stability calculation. 

The factor is obtained by considering the compressive 
force P; = oy (per unit width) in one face and the 
compressive force o,,4/2 in the half core thickness h/2, 


SCIENCES—OCTOBER, 1951 
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The effect of the former on the 
the term 


as indicated in Fig. 5. 
deflection of the face appears in Eq. (5) as 


P,(d*v,/dx*). The suggestion of Williams, Leggett, 
and Hopkins was to include the core force with P,, 
so that Eq. (5) is changed to 
D d'y, + ( de *) dv, , t dq’ se 
o Tas = ies = (9) 
* dx? - * O) dx? 2 dx 


Since this is the only equation affected and is the 
only place where P; or o,,¢ appears in the set of Eqs. 
(5)-(8), the outcome is simply the replacement of 
ore DY or{1 + (o2¢/ors)(h/2t)| in the critical condition 
(9), where the subscript x in o,, has been dropped 
Thus, the corrected value of oa, is obtained from Eq. 
(9) by multiplying the right-hand side by the factor 
(10). 

This procedure implies two assumptions: 

(1) 
on the curvature d*v,/dx* of the face throughout the 
core thickness. This would be correct for the Euler 
type of buckling which prevails for long wave lengths. 

(2) The corresponding transverse force a,,(h/2) X 
(d*v,/dx*) is communicated directly to the face. Actu- 
ally there will be shear forces in the core which will 


The core stress o,, may be considered to act 








partially balance this force. 
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CRITICAL 


The factor therefore represents an approxi:nation 
and probably, in view of Eqs. (1) and (2), gives an 
evercorrection to Eq. (9). In reference 5, the core 
force Ore(h/2) [see Eq. (5’) above] was multiplied by 
a coefficient A less than 1 and, of course, not known 
precisely. We not this coefficient 
that is, we have taken it as unity—because the curves 
given by applying (10) to Eq. (9)—the 
curves in Figs. 4a to 4g-—-show remarkable agreement 
with the exact check points indicated by the crosses. 
This 


( = 


have included 


solid 
The crosses lie almost exactly on the curves. 


values of p 
4d), where the 


for the low 
300 (Figs. 4a 
(10) may be 
hundred per Eq. (9) 
(10) yields remarkably accurate results over the wide 
range of p from 1 to 10,000 and 7 from 0.002 to 0.2 


accuracy prevails even 
E,/E.) of 1, 100, 200, 
the 

Thus, 


correction due to factor several 


cent. corrected by 


that is, pry between 0.002 and 2,000. 


It may be observed that, where the correction is 
significant, the curves are either straight, indicating 
Euler buckling with plane sections remaining plane, 
or nearly straight so that a shear correction is adequate. 
In particular, the correction factor (10) is not needed 
for any of the minimum points, such as M in Fig. 4f. 
It is indicated below that these points give the critical 
stress the buckling. Thus, 
critical wrinkling stress is not affected by compressive 


for wrinkling form of 


stress in the core. 


The factor is merely the ratio of thrust in the faces 
to the total the state. It is 
therefore unity and superfluous when the core takes 


thrust in unbuckled 


only an insignificant part of the thrust. 
With the help of Eq. (4), the factor (10) may be 


transformed into 


2rp[1 — vy. + 2rp(1 — »,?)] 


a = ir 
lL — vf + 4ral(l — vp.) + Sr*p"(l — »,*) 
and further, with »y, = 0, vy = 1/3, into 
2rp(1 + 2rp)/[(1 + 2rp)? — (1/9)] (10”) 


This differs from unity by less than 5 per cent so long 
as pr exceeds 10. Thus, the thrust taken by the core 
and the correction factor can be neglected (with errors 
that is, (E,/hE, 


less than 5 per cent) if pr exceeds 10. 


(3) INTERPRETATION OF A TYPICAL CURVE 

To explain in greater detail the significance of the 
curves presented, we select as an example the curve 
0.02 in Fig. 4f. The plotting is double 
The sandwich plate to which this curve 
applies has, since r = t/h = 0.02, a core thickness 50 
times the face thickness. On Fig. 4f we see that 
p = 1,000—that is, the Young’s modulus of the faces 
is 1,000 times that of the core. The different curves 
have different values of r. 


marked r = 
logarithmic. 


COMPRESSION 
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The horizontal coordinate e,, or rt/X, is r times the 
ratio of face thickness to wave length of the buckling 
Concerning ourselves with the one 
curve ry = 0.02, the /eft-hand (small e,) end of the 
horizontal axis means Jong wave length; the right- 
hand end means short wave length. The 
a,/E, gives the ratio of the critical compressive stress 


deformation. 


ordinate 


in the faces, in the direction of load, to the Young's 
modulus of the faces. Thus, the curve gives the critical 
stress when the wave length is assigned. 
length \ for a given length / of plate (Fig. 2b) is 2/ if 
the plate buckles in the quasi-Euler manner (Fig. 2c) 
several half 


The wave 


in one half wave length. If there are 
wave lengths in the length / as in Fig. 2d, A is 2 //n 
where 7 is the integer number of half wave lengths in 
!. The appropriate value of » is that which gives the 
smallest ordinate. 

Beginning at the left-hand end of the curve for r = 
0.02 in Fig. 4f, we observe that the curve rises for a 
while following a straight line ACB. This corresponds 
to a simple Euler buckling for long wave lengths, plane 
sections remaining plane in the bending. The equation 
for this line (the Euler formula for the composite 


plate) is 
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Oy 1 + [2p/(1 — v/*)](3 + 6r + 4r’?) ( l ) (e,2) 
Ey, 1+ 2rp{(1 + 2rp)/(1 — v2 + 2rp)] \3rp/ ** 
(11) 


which on log-log paper gives a straight line with slope 
2. 

Proceeding to the right along the curve, we observe 
that it begins to droop away from the straight line at 
C at a wave length corresponding to e, about 6 X 
10-4. This, at first, can be attributed to additional 
deformation due to a shear and can be accounted for 
by a corresponding correction* to the Euler formula 
[Eq. (11)]. Going further along the curve in Fig. 4f, 
we find it becomes horizontal (point E). Plots of 
the shear correction, not reproduced here, show that 
the range of validity of such a correction reaches but 
ends at point E, in thisinstance. Passing the maximum 
point E, the curve dips to a minimum at M. This 
occurs at about e, = 8 X 10~°, or when the half wave 
length \/2 is about 20 times the face thickness /, 
which would be “wrinkling’’ or “corrugation”? with a 
considerable number of waves in the length of a plate 
or panel of practical proportions. A fairly fine adjust- 
ment of wave length is then effected by addition or 
subtraction of one half wave, and, therefore, the wave 
length of the minimum point can always be closely 
approached when it is small compared with the plate 
length /. Accordingly, the part of the curve which 
lies above the minimum, the broken portion DEM, 
must be discarded for practical evaluations. This 
leaves us with the minimum point M itself, correspond- 
ing to the wrinkling stress and wave length, together 
with the left-hand part ACD of the curve, which in 
this case could be adequately represented by a shear 
correction to the Euler formula, Eq. (11). 

Returning to the minimum point and proceeding to 
the right, we find the curve rising toward an asymptote. 
This asymptote, marked PN, represents Euler buckling 
of the two faces as independent sheets without a core. 
The approach of our curve, therefore, means that at 
extremely short wave lengths the core ceases to in- 
fluence the faces. 

The amount that the curve stands above the asymp- 
tote PN shows how much the presence of the core 
increases the critical stress in the faces. On the other 
hand, the amount that the curve falls below the initial 
asymptote ABC shows how much the core falls short 
of compelling the sandwich to bend in the desirable 
“‘plane-sections-remain-plane’’ manner. These com- 
parisons presuppose a common wave length. 
face sheets of length /, without core, will have \ = 2/. 
With the core we shall have either ) = 2/ or a much 
shorter wave length \ = 2//n corresponding to wrin- 
kling at the minimum point M. 


Two 


* Several such corrected theories for this ‘‘quasi-Euler’’ 
range have been proposed (see, for example, references 4-8, 
10, and 11). These are based on various simplifying assump- 


tions, and the ranges of validity are not easily defined. 


AERONAUTICAL 


SCIENCES—OCTOBER, 1951 

The line QR on Fig. 4f represents the Euler buckling 
of a monolithic plate of thickness 2/—that is, the 
single plate formed by combining the faces with no 
core between. The critical stress is then four times 
that of the plate of thickness ¢ represented by PN, 
Thus, the difference between the curve ACDEMN 
and QR shows the net effect of using the metal in the 
form of two faces (plates) separated by a core 50 times 
as thick as one face (r = 0.02). To the left of QR 
the ordinates of the curve are greater and show the 
advantage of the core for the corresponding wave 
lengths. 

Again, this immediate comparison is only valid 
when the sandwich has the same wave length as the 
monolithic plate—i.e., \ = 2/—and not when the 
sandwich wave length is that of wrinkling. 

It will be observed that several other curves on 
Fig. 4f come close to the minimum point M of the 
curve ACDEMN for r = 0.02. One of these other 
curves (marked ‘.02 sym.’’) is valid for the same 
sandwich but buckling in the symmetrical mode of 
Fig. 2e. Another is marked ‘.002 sym. + asym.” 
to indicate that it is valid for both modes for a sandwich 
having ry = 0.002. The merging of these curves at 
point M can be explained by considering first the 
sandwich with r = 0.02—that is, h/t = 50—and then 
an increase in core thickness by a factor of 10, resulting 
in r = 0.002, other things remaining unchanged. 
The merging shows that this change has not altered 
the critical stress for wrinkling. The interpretation is 
that the short wave length results in a deformation 
in the core which dies out rapidly going in from either 
face, leaving a middle zone unaffected by the buckling 
deformation. Increase of core thickness then merely 
increases this idle zone and makes no difference to 
the critical stress in the faces. Whether the deforma- 
tion is asymmetrical as in Fig. 2d or symmetrical as in 
Fig. 2e is of no consequence, since the two active 
zones are rendered essentially independent by the idle 
middle zone. Under these conditions, the wrinkling 
stress can be found by considering a single face on a 
semi-infinite core. 

Next, examining the curve marked ‘‘.05 sym.,’’ we 
find its minimum a little higher. The value 0.05 for 
r would represent a thinner core, other things being 
unchanged, and the higher minimum means that the 
core deformation in the wrinkling is no longer in two 
separated layers close to the faces but affects the whole 
thickness. Each face and the core close to it can no 
longer be considered as independent and equivalent to 
a single face on a semi-infinite core. The curve higher 
up, marked “.2 sym.,”’ 
The curve “.2 asym.”’ lies below and 


shows the same thing to a 
greater extent. 
has no minimum; therefore, any sandwich for which 
r = 0.2 and p = 1,000 will have no wrinkling mode. 
It indicates buckling with one half wave in the length /. 
The same applies to the curve marked ‘.05’’ below 


ACDEMN. 
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CRITICAL COMPRESSION 


(4) Exact RESULTS 


The theory relied on for ‘“‘exact’’ results is that of 
This takes into account the destabilizing 


reference 3. 
The paper of 


effect of compressive stress in the core. 
reference 2 proceeds from essentially the same basic 
equations but is more complex in not accepting ordinary 
plate theory for the faces. 

Following Eqs. (26) and (27) on page A-256 of 
reference 3, it is indicated that the critical stress may 
be found by solving a transcendental equation. Since 
this equation is intractable, a different procedure is 
followed in reference 3, resulting in an explicit but com- 
plicated formula for the critical stress. This is Eq. 
(46c) on page A-260 of reference 3 and is obtained only 
by neglecting certain terms expected to be small. 
On evaluation of this formula, indications’ were en- 
countered that the neglected terms were not, in fact, 
insignificant over the whole range of practical values. 
It was therefore necessary to resort to the trans- 
cendental equation for exact results. This equation 


is 
Of _ e, ~" ef + = +- Gal, + ae 8 (12) 
Ee Ae = 99°) L— vw + AT + Joli 
where 
: (1 — »f)[1 + (1 — 2»,)s?] 
G = - 7 + 
1(1 — v,“)pse, 
(1 — 2v.)(1 + s) l l 
+0) 
2(1 — »v,) l1—»\Xs 
(1 — v»?)R(1 + 5)? 
G = * + 
11 — Vo~)pSey 
j(1 — 2y,)R(1 + s) | k(1 +s) ]) 
{ 2(1 — ») ~ Py § i 
‘i ; (13) 


tanh (ah /r 
T,. = tanh R(ah/d) 


and s is related to o,/E, by 


a ry — Ve + 2rp(1 — »v,”) ( S ) 14 
- = (14) 
Ey F — ve + 2rp(1 — vyr,) 1+ », 


Here, the special values for v, and », (1/3; and 0) have 
not yet been inserted, and the notation follows refer- 
ence 3. The value of o,/, corresponding to a chosen 
e, can be found by trial when p and ¢ are given. The 
“exact’’ points obtained are shown as crosses (x) on 
the figures, each labeled with the value of r it belongs 
to. These crosses are thus check points that serve 
to show the correctness or incorrectness of curves ob- 
tained from less complete but simpler theories. The 
calculations required for these check points were too 
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Chart for critical wrinkling stress and wave length. 





Fic. 6. 


heavy to permit the plotting of complete ‘‘exact”’ 


curves. 


(5) A WRINKLING CHART 


The wrinkling stress given by a minimum point 
such as M in Fig. 4f depends on p and sometimes also 
on r. It will be designated by o,,. It occurs at a 
definite value of e,, the abscissa in Figs. 4a-4g. Such 
a point can be transferred to a diagram such as Fig. 
6 with oy,/E; as ordinate and e, = mh/X = re, as ab- 
scissa. If other minimum points for the same p but 
different values of r are also obtained and recorded on 
Fig. 6, the locus of the points is a curve for this value 
of p, along which are points with different values of r 
attached. We then have one of the curves in Fig. 6 
which appear as almost horizontal straight lines, the 
value of p being indicated at its right-hand terminus. 
The intersecting curves are curves of fixed 7, and the 
values of r are indicated along the top and along the 
right-hand edge of the diagram. Fig. 6 contains all 
wrinkling (minimum) points shown in Figs. 4a to 4g 
with others added (from computations for the minimum 
only) to complete and refine the chart. Its purpose 
is to show, by inspection of this one chart, whether 
wrinkling can occur or not and, if it can, at what 
stress and wave length. Thus, being given p and 7, 
we look for the corresponding curves in Fig. 6 and read 
off from their intersection point, on the left-hand and 
If there is no 
no wrinkling 


base axes, the values of o,,/H, and éy. 
e.g., p = 1,000, r = 0.05 
(minimum) point exists. This corresponds to the 
form of the curve for r = 0.05 in Fig. 4f; it shows no 
When no wrinkling can occur, the 

that is, the curve of the type 
must be of the general form of 


intersection 


minimum point. 
curve of o,/ Ey, versus e, 
shown in Figs. ta—4+g 
those in Fig. 4b or that in Fig. 4g marked r = 0.05. 
The ‘‘horizontal’’ curves for given p in Fig. 6 are 
truly horizontal except for the slight droop at the left- 
hand end. Except here, therefore, a change in the 
value of , moves the wrinkling point horizontally, 
changing the wave-length-thickness ratio but not the 


wrinkling stress ratio o7,/H,. This corresponds to the 
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Fic. 7. Several simple wrinkling stress formulas. 
independent wrinkling of the two faces on a relatively 
The droop at the left- 
hand end indicates that, for the corresponding range of 


thick core discussed above. 


r values on any one ‘horizontal’ curve, the core is 
not so thick that each face can be regarded as a sheet 
1,000, for 
instance, the wrinkling of the faces is independent, 


on a core of infinite thickness. For p = 
as though the core were infinitely thick, only when + is 
about 0.025. 
wrinkling stress can fall short of the infinite core value 


less than For larger values of r, the 


by as much as 15 per cent. 


(6) WRINKLING STRESS FORMULAS 


On the straight parts of the “horizontal’’ curves in 
Fig. 6, the ordinate is independent of e, and r and, 
hence, of e, and 7. Thus, for the region covered by 
these straight parts, oy,/H, is dependent on p only. 
A simple formula can be found to cover this region. 

Regarding o,,/H,, in general, as a function of e, 
and r, independence of r means that the core is suf 
ficiently thick to behave as a semi-infinite foundation. 
This corresponds to putting 7° = tanh e, 
(9). 


| in Eq. 
As noticed at the end of Section (2), it is not 


necessary to consider the correction factor (10). We 
have then 
OF é? 
at “ie 
E,; 3(1 — v/) 
I 4 l 2p, Le? 4 l vy 
e e,* 
l VY, ; (1 — vy.) pe, . 
15) 
(3 — 4y.)(1 + »v,) 
(1 vy") + pe; 
l— yp, 


To find the wrinkling stress, we have to minimize 
this with respect to \—-that is, with respect to e,. 
It may first be simplified, however, since examination 
of the minimum points in Figs. 4a to 4g shows that they 
occur for pe, greater than 50, except where the minimum 
point is above the practical range for o,//, (0 to 0.08). 
This justifies neglect of the last term in the numerator 
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of Eq. (15) and the first term in the second denominator 
leaving 


1 — 2», 
° l + ey + er” 
Or er” l — ¥% 
Ey 3(1 — p,’) 1+ »y, 
(3 — 4») pe; 
—"s 
or, with », = 0, », = '/3 
3 ,, lteter 
E, Ss ; Spe, 
1/. 


This has a minimum at e, = 0.726p" 


OT fu ee 2 
= 0.655 ag 
E, “ ( 


’, yielding 


0.51p”? + sa 
(16 


p+ 0.39 
which can be replaced to a good approximation by 
Crw/ EL, = 0.815p 


This is shown on Fig. 7. 

Several other formulas of this kind, obtained by 
other approximations and other basic methods, have 
been proposed. They are as follows, changed from 
the original notation to the notation of this paper (the 
distinction between plane strain and plane stress for 
the core vanishes since we take vy. = 0): 

Gough, Elam, and de Bruyne! 

Trw/Ey = 0.65p 
Hoff and Mautner!! 
Ope /E, = 0.75p 
- 3 . . ‘ 
(adapted from o,, = 0.91 V E;E.G, 


Cox’ 





Sry/ LE, = 0.629 


6 
[adapted from o,,/H; = 0.69 Y E,E,°/E,;* correspond- 
ing to Eq. (33d) of reference 9} 
All of these formulas are plotted in Fig. 7. 
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Allowable Compressive Stresses in 


Aircraft Structures 


B. E. GATEWOOD* ano E. L. WILLIAMS 
U.S.AF. Institute of Technology and Curtiss-Wright Corporation 


SUMMARY 


This paper shows procedures for determining the allowable 
buckling and maximum compressive stresses in both the elastic 
and plastic regions for structural materials for which compressive 
stress-strain curves are known. In particular, curves are shown 
that are applicable for the aluminum and magnesium alloys. 

Detail procedures for obtaining the allowable compressive 
stresses in both the elastic and plastic regions are indicated for 
the following types of aircraft structure under compressive 
loads: (1) compressive buckling of flat plates; (2) compressive 
maximum of flat plates; (3) shear buckling of flat plates; (4) 
bending buckling of flat plates; (5) column stress for structural 
elements with stable cross sections (no local buckling); (6) 
interrivet buckling of flat plates; (7) compressive buckling of 
formed and extruded sections with local instability; (8) com- 
pressive maximum or crippling of formed and extruded sections 
with local instability; (9) column stress of formed and extruded 
sections with local instability; (10) bending buckling of formed 
and extruded sections with local instability; (11) bending maxi- 
mum of formed and extruded sections with local instability; 
(12) compressive buckling of flat plates with stiffeners; (13) 
compressive maximum of flat plates with stiffeners; and (14) 
column stress of flat plates with stiffeners. 


INTRODUCTION 


gene AND MELCON?® discussed the use of non- 
dimensional buckling curves to determine allow- 
able buckling and column stresses for plate elements 
in aircraft structures. It is the purpose of this paper 
to extend and modify the procedures of reference 2 
to include the latest results on the buckling and failure 
of plates, formed or extruded sections, and stiffened 
plates in both the elastic and plastic regions. If the 
method of reference 3 is extended to stiffened panels 
and incorporated into the procedure of reference 2, 
then the maximum stresses (crippling) of many formed 
and extruded sections, as well as many stiffened panels, 
can be read directly from curves without having to 
compute laboriously the stress for each element and 
take the average of all the elements. For the stif- 
fened panels it is not necessary to compute the effective 
skin. 

To determine the suitability of the curves resulting 
from the above procedure for analysis of component 
parts of the aircraft structure, the results of approxi- 
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mately 1,700 tests were compared with values predicted 
by the curves. The scatter appears to be less than 
that usually encountered in compression tests of air- 
craft components. 


NONDIMENSIONAL STRESS-STRAIN CURVES 


References 6 to 10 show nondimensional stress- 
strain curves for a large number of stress-strain tests. 
The curves are made nondimensional by plotting 
a/oy against eH /o,,, where oy is the yield stress defined 
by the secant line of slope 0.7E. There is little variation 
in the nondimensional curves, even though the dimen- 
sional curves vary widely. In fact, if the nondimen- 
sional compressive stress-strain curves for 75ST, 
24ST, and R301T in references 6 to 10 are plotted 
together, there is little variation between the alloys. 
References 11 and 12 show that these nondimensional 
stress-strain curves, as well as those for other materials 
such as steel and magnesium, can be represented by 


the equation 


l eR o ssa" 
_= = l += - (1) 


B? Oey Tw 7 \Oey 


where is the stress-strain shape parameter given in 
references 2, 11, and 12 and where B is a parameter 
for constructing curves. Fig. 6 of reference 11 and 
Fig. 14 of reference 12 show graphs of Eq. (1) for 
selected » from 2 to ~. A study of the average stress 
strain curves in references 6 to 10 and of Table 1 in 
reference 2 indicates that almost all of the aluminum 
alloys ordinarily used in aircraft structures can be 
approximated by one stress-strain curve: = 10 
in Eq. (1). References 13 and 14 also indicate that 
the magnesium alloy extrusion of J-1 and O-1HTA 
can be approximated by the same curve, » = 10. 
Some of the extrusions with a high yield point, such 
as 75ST, have a value of m greater than 10. However, 
if m = 10 is used for them, the results will be conserva- 
tive up to the yield point; since these materials do 
not reach much beyond the yield point in elongation, 
the stress-strain is usually cut off at, or slightly above, 
the yield point. Fig. 1 shows Eq. (1) for m = 10 with 
the cutoffs marked for alloys with high yield points. For 
design, the o,, defined by 0.7E is close enough to the 
usual o,, defined by the 0.2 per cent offset line so that 
either can be used. However, in the nondimensional 
plot, the strain associated with 0.7/ must be used. 
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NONDIMENSIONAL BUCKLING CURVES 


Reference 


buckling stress of columns and plates in the plastic 
range by introducing a factor y so that nF is the 
n depends upon the 


effective modulus of elasticity. 


cross section and length of the section and upon the 


tangent modulus /;,, and the secant modulus F,... 
Since the stress-strain curve as represented by Eq. 
(1) corresponds to n = E,../, then the nondimensional 
buckling equation for any 7 can be written as 


nu—1) 


l Cras Fee (Ger/ Oey) [1 + (3/7) (Ger/ Oey) 


B? "tw nE/Esec 


(2) 


where o,, is the buckling stress. Note that the buck- 
ling curve can be obtained for any stress-strain curve, 
whether the equation is known or not, by dividing the 
strain by E/E ec. 

Reference 15 gives equations for n for plates, flanges, 
and columns in compression; reference 16 gives 7 for 
shear buckling; reference 17 gives y for composite 
sections; reference 18 gives a comparison of the method 
of calculating y in reference 15 to other methods of 
obtaining the buckling curves. Reference 49 gives a 
method of calculating 7 different from reference 15, 
but it does not check tests for plates as well as that of 
reference 15. It appears from these references that 
one value of 7 can be chosen for most plates, formed 
and extruded sections, and stiffened plates provided 
they fail by local instability or crippling, while another 
value can be chosen for columns. These values are 
approximations, but the theory itself is also an ap- 
proximation so that the assumption of two values of n 
for aircraft structures seems justified. In fact, refer- 
ence 2 uses one value only, » = /j,,/£, for plates as 
well as columns; reference 4 uses n = F,,,/ for plates. 
The approximate values selected are 


mise (; 1 /t 3 us) f - | 
= ton) 
, | eS. oN | 4 Bos shal ci 


n = Etan/E (for columns) 
(3) 
Reference 48 shows that this value for columns with 
n = 10 is accurate enough for design. 
If Eq. (3) be substituted into Eq. (2), then 


5 - n 1 
>) cr 
1 + “( ) 
l es Car 6\Gcy 


Tey I “tan 


Le. od 
oy +s5 


4 
“sec 


9 


for plates and sections in compression and plates in 


‘ n—1 
l Ser Ser 3 [Ger , 
_= = 1 + =n (5) 
B , No cy Ocy ‘ Tey 


for columns, where 


shear and 
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15 gives a method of calculating the 
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a 1 + (3/7) (Ger/Ocy)"~' 

" = - n - (6 

Pigec 1 + (3/7)n(oer/oey)"~ | ’ 
Eqs. (4) and (5) are plotted on Fig. | for n = 10, 


By using Eqs. (4) and (5), the standard buckling 
and column equations can be extended to the plastic 
range 


o 


cr (= (‘) a. : 
NO cy - Ca ) b} = BB M 
Cc =)( p ) | 
= = (S 
NO cy Cig J VE B? 


Thus, the buckling stress, Eq. (7), for any plate, 
formed or extruded section, or stiffened plate can be 
read from Fig. 1 provided o,, of the material is known 
and A can be obtained. Any column stress can be 
read if ¢,, and L’ are known [Eq. (8)]. It should 
be noted that Johnson’s parabola for short columns 
differs little from the column curve, Eq. (8). In fact, 
from Eq. (8), 


O- l l 
= = (9) 


NT ce 4[1 — (o¢/ac) | B* 
For sections that do not fail by local instability (crip- 
pling), O¢ = Gey; otherwise, o;¢ is the crippling stress 
of plates, formed or extruded sections, or stiffened 
plates. Eq. (9) is shown on Fig. 1. 
Paragraph 14.9 of reference 5 discusses the effect 
of eccentricities on columns for a stress-strain curve 
with » = o-—that is, the stress is elastic up to a maxi- 
mum stress at which failure occurs. By use of Eq. 
(14.39) in reference 5 and Eq. (8) above, one can 
write 
I Te i Oe/ Cec 
= = (10) 


B* NO ce t | ea - 
,cos—! 
| (Gcc/o,) — 1]I 


where o,- is the yield stress o,,, or the crippling stress 
Fcc, OF SOME Cut-off stress; where e is the eccentricity; 
and where a = J//Ac with J/c section modulus and 
A section area. Fig. 1 shows Eq. (10) for several 
values of e/a. 

From reference 2, the formula for interrivet buckling 


can be written as 


O- for? ie =e) l , 
= = ( ) 
NO cy a ( S B* 


cy 


with s as rivet spacing. o,/o,, can be read on Fig. | 
from either Eq. (8) or (9). 

It should be noted that the sheet must be formed in 
the O-condition and heat-treated so that it will have 
the same o,, throughout. On basis of references 3, 
19, and 20, extrusions have practically the same ¢, 
throughout. If formed in the T-condition, o,, will 
vary, giving a higher stress in the corners than on the 


flats. The most satisfactory procedure for forming in 
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Fic. 1. Allowable buckling stresses. 
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Fic. 2. Calculation of B. 


T-condition is to use the test buckling curve from Fig. 
14 of reference 21 for sections formed in the T-condi- 


tion (see Fig. 4). 


K VALUES FOR BUCKLING STRESSES 


The N.A.C.A. has shown that Eq. (7) can be used 
for the buckling of formed sections, extruded sections, 
It is only necessary to determine 
Then Fig. 1 can be 
K values for dif- 


and stiffened plates. 
the correct A for the section. 
used to obtain the buckling stress. 
ferent sections can be calculated from data and pro- 
cedures in references 23, 24, and 25. A values for 
flat plates are given in Fig. 14.25 of reference 5. A 
values for compression of channels and zees are in 
Reference 19 shows A for com- 
K for H-sections in 
5 shows 


references 19 and 26. 
pression for rectangular tubes. 
compression is in references 19 and 26. Fig. 
K for compression of angles as calculated from refer- 
ences 23 and 25. Fig. 6 shows A for bending of chan- 
nels and zees as calculated from references 24 and 25 
K for buckling of compression panels with Z or T 
stiffeners is in reference 28. 

The N.A.C.A. has made 
stresses on formed and extruded sections and for stif- 


many tests for buckling 
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Fic. 4. Allowable crippling stresses 


fened plates. If these test values are compared with 
the plate buckling curve in Fig. 1, agreement within 
the test scatter is obtained. Test ae values 
are given in references 13, 19, 20, 21, 31, 32, 33, 27, 
35, 36, and 18 for most of the aluminum and mag- 
nesium alloys used in aircraft structure. 


MAXIMUM ALLOWABLE STRESSES 


The maximum allowable stress or ultimate allowable 
stress is the allowable crippling stress for local failure 
of a section and is the allowable column stress for 
column failure in either the long or short column range. 
In most cases the allowable crippling stress, which may 
be higher than the buckling stress because of the fact 
that the edges and corners of the section may be able 
to carry higher stresses than the flat elements and which 
may involve yielding in the corner, must be deter- 
mined from empirical formulas or curves obtained 
from tests. Reference 29 gives a theory for flanges 
which checks test data with little scatter. 

On the basis of tests, several procedures have been 
set up to determine curves that may be used to obtain 
crippling stresses for design of any type section. 
Reference 2 gives a procedure of constructing the design 
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TABLE 1 
M Values 
M 


Extrusion or formed 
in O-condition 


Formed in 
T-condition 


b/t of Web 


0-25 0.85 0.80 
40 0.75 0.67 
75 0.67 0.60 


value for any section and material by use of a combina- 
tion of loads carried by the flat and curved elements 
of the section. The structures manuals of many air- 
craft companies give curves based on tests to build 
up the design values from flat elements only. Refer- 
ence 3 sets up curves from tests on zees and channels'.*! 
for crippling or maximum stress based upon the yield 
stress and the buckling stress. 
procedure for the effective area to use on stiffened 
panels. 

The method of reference 3 can be modified and ex- 
tended by the use of Eqs. (4) and (7) above to give 
the maximum stress for any section for which A can 
be obtained. Reference 3 shows that for a certain 
range of values of o,, the maximum stress omar. can be 
expressed in the form 


Reference + gives a 


o Tma Ser 
c=—"=M \ (12) 


cy N0y 


where M depends upon the geometry of the section. 
Fig. 4 shows a plot of Eq. (12) for selected M values. 





“M"and"K" Values for Angles in Compression 
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cond. formed in"O" 
cond. 








0.80 



































b,/b, 











Fic. 5. M and K values for angles in compression. 
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SCIENCES OCTOBER 19051 
The maximum stress is taken as equal to the buckling 
stress to the left of the intersection of Eqs. (12) and 
(4). For zees, channels, and H-sections, the A/ values 
as correlated with 6/t were taken from Fig. 11 of 
reference 3 (see Table 1) 

For plates use WW = O48, 
interpolated on Fig. 4. 


These values may be 
M can 
be evaluated from several tests or from analogy to the 


For other sections, 
zee and channel sections. For angles in compression, 
M as correlated to the longer leg of the angle should 
correspond to about one-third of the }/t for channels 
or zees (see Fig. 5). For bending of channels and zees, 
it appears that 6/t of the compressed leg should be used, 
with M being that for a channel or zee with the same 
b/t. For stiffened panels in compression, the test 
data of reference 27 indicate that 1J = 0.80 gives 
suitable results regardless of b/t of the skin (for stiffener 
sizes used). 

For sections for which A is not available, the maxi- 
be obtained from the 


elements of the sections as in reference 2. Fig. 3 


mum or crippling stress may 


shows curves for obtaining B for the elements of the 
section for any alloy. This B is to be used for buckling 
of the elements on the curve for sections in Figs. | 
or 4. Fig. 3 is constructed from the curves for 24SRT 
of reference 2 except that buckling for sections in Fig. 
| is used instead of the tangent modulus. A compari- 
son of Figs. 2 and 3 shows that there is no fixed value 


of K for the elements of the sections as there is for the 
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Fic. 6. Mand K values for channels in bending. 
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COMPRESSIVE STRESSES 


section as a Whole. This would seem to indicate that 
Fig. 3 should not be used for sections that differ ap- 
preciably from those tested to obtain Fig. 6 of reference 
9 Note that any set of test curves for elements of a 
section for a given alloy may be converted to a form 
similar to Fig. 3 by using Fig. 1. 

Section 14.14 of reference 5 points out that columns 
subjected to local crippling failure follow approximately 
Johnson's parabola [Eq. (9)], where o,- is the crippling 
stress. Thus, when o,, is obtained as the maximum 
stress in Fig. 4 or as the average stress from the ele- 
ments in Figs. 3 and 1, then it is used to obtain a B 
for reading the column stress from Johnson's parabola 


in Fig. |. 


ALLOWABLE MAXIMUM STRESSES FOR STIFFENED 
PLATES 


Because of riveted joints, different o,, for stiffeners 
and plates, and eccentricities, the allowable maximum 
stress for a stiffened panel is more difficult to deter- 
mine than that for a formed or extruded section. 
Reference 47 shows the large scatter present in panel 
tests and shows the effects of the size and spacing 
of the rivets in reducing the maximum average stress 
for the panel. Although the curve of reference 47 for 
the reduction in stress is for stiffened panels that fail 
by crippling, a few tests at Curtiss-Wright Corporation 
show a reduction in strength in the column range for 
stiffened panels. For local instability or crippling of 
stiffened panels use 


Tmax. = RiGee (13) 


with o,, from Fig. 4 or Eq. (12) and with MZ = 0.80 
and R, from reference 47. For column action of 
stiffened panels use the above omar. in place of oa, in 
Figs. 1 or 4 and for long column use 


Tet Ror’ E l (14) 
= ed 2 == : 
Omar. Omaz.(L /p)° B? 
For short columns use Johnson’s parabola 
— B? l ub 
~“j-—si< (15) 
Tmar. { 4( oer Tmar ) 


where B is defined in Eq. (14). A few tests indicate 
R; may be about 0.75 for flush rivets and close to 1.00 
for round-head rivets closely spaced. It also appears 
that R, depends upon the location of the neutral axis 
of the panel. If the neutral axis is near the joint of 
stiffeners to skin, then R» is larger than when the 
neutral axis is near the center of the stiffeners. Let 
Rx, be the value of Rs for neutral axis at skin stiffener 
joint and let Ro» be the maximum value of R» for the 
neutral axis near the center of the stiffeners. For 
the variation from R» to Re» use J as calculated for the 
stiffener about the skin stiffener joint until o,;, reaches 
the cutoff value given by Re». 
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There is a further complication in using Eqs. (14) 
and (15) on stiffened panels tested in a test machine. 
For panels with large deep unsymmetrical stiffeners 
closely spaced the neutral axis is not parallel to the 
plate. This introduces eccentric loading on the end 
stiffeners. This eccentric loading will cause the panel 
to fail at stresses below those given by Eqs. (14) and 
(15) [see Fig. 1 and Eq. (10)]. 

To obtain o,, for stiffened panels when KA is not 
known, use can be made of an effective skin area acting 
with the stiffener. 


Reference 4 gives the effective 


width 2w in terms of the stiffener crippling stress as 
os, = KE(t/2w)? (16) 
where A is for the panel. This does not help when K 


is unknown. Reference 5 recommends the effective 


width as 
Os, = 2.89E(t/2w)? (17) 


Note that both Eqs. (21) and (22) can be extended to 
the plastic range by using nF for E. With » for the 
effective skin the same as 7 for the stiffeners Eqs. (16) 
and (17) can be written as 


ow ? f. | E cr 
= = vk ( ) 8. ( ) (: ) (18) 
b, b, Tey/ st\Fee/ st 


with B,, > Beyrog at intersection of o,,< curve and cut- 


off line in Fig. 4; for By, < Beton, use Beutog- Eq. 
(18) is shown in Fig. 4 against 


B «= V (Gee) s1 KE(6, ts) 
ByV (oer Tey) st 
Note that 
BuiV (Ger ey)s1 = | 


in elastic range and that any stiffener o,, may be put 
on Fig. 4 to obtain (¢,,),, and B,,. 

To compare values given by Eq. (18) to test values 
and to values obtained from Eq. (12) or Fig. 4 use 


2Qw we A, A, 
> we ate (1 + “) _— (19) 
b, Ost byt, byt, 


where (o,-), is test value or value from Fig. 4, and 
A,, is stiffener area. Comparisons of values of 2w/b, 
for several stiffened plates (Z stiffeners in reference 35) 
as given by Eq. (18) with A = 2.89 (Eq. 17), A from 
reference 4 and A from reference 28, and as given by 
Eq. (19) with (¢,,), from Fig. 4 and from tests show 
that Fig. 4 checks the tests better than the other 
methods. Eq. (18) with A from reference 4 is high 
[Eq. (16) is also high]; Eq. (18) with K = 2.89 is 
usually low, but, along with Eq. (17), checks well in 
the elastic region; Eq. (18) with A from reference 28 
checks for certain stiffened plates with large 6,/t,. 
Hence, use Fig. 4 when K is known, and use Eq. (18) 
with A = 2.89 when A is not known. 
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Fic. 7. Test values for maximum stresses. 


TEST VALUES FOR MAXIMUM STRESSES 


If the K and .W be obtained for the section and if 
the yield stress o,, is specified, then the buckling stress 
and maximum local failing stress can be read from 
If the structural element is a long 
if the structural element is a short 


Figs. 1 and 4. 
column, use Fig. 1; 
column, use Fig. | with o,, from above in place of o;,. 
For stiffened panels use Fig. 4. 

Using the above procedure, allowable maximum 
stresses were calculated for about 1,700 test specimens 
and were compared to the actual test stresses. Fig. 7 
shows the frequency of tests in intervals of 0.05 of 
Frest/%eal. aS plotted against opes¢/¢cai. Included in 
Fig. 7 are the crippling stresses of 827 extruded H- 
sections, extruded and formed Z-sections, extruded and 
formed C-sections of 75ST, 24ST, R-308T, O-1IHTA 
magnesium, 14ST, 17ST; 1% !9—?1,31—33,38,51 the crippling 
stresses of 87 angle sections of 75ST; ** °° the maximum 
stresses of 813 stiffened plates (stiffened by Z’s, T’s, 
C’s, and angles) of 75ST, 24ST, 14ST, and ZK60A 
mg. !*:°7,54~ 86,4450 Ror the stiffened plates, in Eqs. 
(13), (14), (15), Ri = Re. = 1.00 except for the 24ST 
panels of reference 35; in reference 35, R; = 1.00, 
Ry = 0.75, and Re = 0.95. Also, for the panels of 
reference 36 with deep stiffeners closely spaced, the 
eccentricity on the end stiffener was used in Fig. 4 
for the panels with large L/p; for low L/p, the amount 
of eccentricity used was reduced. It should be noted 
that this problem of eccentricity arises because of the 
method of testing and ordinarily will not occur on the 
actual airplane panels. 


SUMMARY OF PROCEDURES FOR ANALYSIS 
In conclusion, the procedures for compressive 
allowables for different types of aircraft structure that 
discussed above may be summarized as 


have been 


follows: 
(1) Compressive Buckling of Flat Plates. 
curve in Fig. 1 with 


Use plate 
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B = Vo;,/KE(b/t) 


(K on page For L/t > 2, K 
0.385 for one edge free, 3.62 for simple support, and 


372 of reference 5). 


6.35 for clamped edges. 

(2) Compressive Maximum of Flat Plates. 
Fig. 4 with 1J = 0.48; B and K as for compressive 
buckling. 

(3) Shear Buckling of Flat Plates. 
of o,, (for aluminum alloys can approximate o,, as 
0.6 oy) in procedure for compressive buckling (+ being 


Use 


Use oy, instead 


smaller dimension of the plate) with A for shear on 


page 394 of reference 5. For L/b > 4, K = 4.9 for 
simple support and 8.0 for clamped edges. 
(4) Bending Buckling of Flat Plates—Use same 


procedure as for compressive buckling with A’ for 
bending on page 395 of reference 5. For L/b > 0.6, 
K = 21.6. 

(5) Column Stress for Stable 
Sections.—Use Euler's column and Johnson's parabola, 


or eAc/I curves for eccentric loading, in Fig. 1 with 


Sections with Cross 


B= WV Cu Kad ‘or (L’ p) 


(6) Interrivet Buckling of Flat Plates—Use Euler's 
column and Johnson's parabola in Fig. | with 


B= WV e6./E(1./2/2)(s/0) 


(7) Compressive Buckling of Sections with Local 


Instability.—Use plate curve in Fig. 1 with 
B = WV ¢q/KiE(/t) 


K in Fig. 5 and references 19 and 26. 

(8) Compressive Maximum or Crippling of Sections 
with Local Instability.—Use Fig. 4 with M from Fig. 5 
and from Table 1. If A is not known, use Fig. 3 to 
obtain B for each element of the section; use each B 
on the plate curve of Fig. 1 to obtain the stress for the 
element; average all elements (}>¢A/A) to obtain 
o-- for the section. 

(9) Column Stress of Sections with Local Instability. 

Use Euler's column and Johnson's parabola, or 
eAc/I curves for eccentric loading, in Fig. 1 with 


B = V o,-/E(1/mr)(L'/p) 


Occ = Crippling stress from item (8) above. 
(10) Bending Buckling of Sections 
Instability._Use plate curve in Fig. | with 


with Local 


B= Vow/KE(b/t) 


K in Fig. 6. 

(11) Bending Maximum or Crippling of Sections 
with Local Instability —Use Fig. 4 with M from Fig. 6. 
If A is not known, use Fig. 3 on selected (?) elements 
of the section and proceed as for compressive crippling 
in item (8) above. If the section has lateral instability 
then the above value must be reduced. 
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(12) Compressive Buckling of Stiffened Plates.—Use 
plate curve in Fig. | with 


B = Vo.,/KE(b/t) 


K in reference 28. 

(13) Compressive Maximum of Stiffened Plates. 
Use Fig. 4 with M = 0.80. If A for stiffened plate is 
not known or is not applicable, obtain ¢,, for the stif- 
feners as described in items (7) and (8) above and read 
the effective skin 2w/d, acting with the stiffener from 
Fig. 4 with 

V (Gcc)s2/KE(6,/ts) 
Ba (Bis) Center 


B= 


where A = 2.89 for middle stiffener* and A = 1.44 
for end stiffener (pages 374 and 375 of reference 5). 
Then 

(Occ) p 1 — (2w/bd,) 


esse, 
an 1 + (A,g,/dgt;) 


which is shown in Fig. 4 with A,,/b,f, as parameter. 
Omar. for panel is Rio,,, where R; is rivet factor. 

(14) Column Stress of Stiffened Plates.—Use Euler's 
column and Johnson's parabola, or eAc/J curves for 
eccentric loading, in Fig. 1 with 


mar l q," 
B= q¥— ( ) 
RE wr\p 
R, = rivet factor. 
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Summary of Recent Experimental 
Investigations in the N.O.L. Hyperballistics 
Wind Tunnel’ 

PETER P. WEGENER?t 
UR Nast Obie Labontny 


SUMMARY 


The N.O.L. 12- by 12-cm. continuous blow-down hyperbal- 
listics wind tunnel is briefly described. Limitations of present 
supersonic wind-tunnel techniques as applied to hypersonic wind 
tunnels are discussed. Measurements of pressures in the nozzle 
show that air condensation occurs if flow pressure and tempera- 
ture fall below the saturation line of air. This condensation can 
be eliminated by preheating the supply air of the wind tunnel. 
Tentative experimental criteria on minimum supply tempera- 
tures to minimize air-condensation effects on the flow in the test 
section are given. The highest Mach Number achieved to date 
at N.O.L. on this basis is 8.25. The effects of air condensation 
on various standard measurement methods for the determination 
of Mach Number are given. It is found that only static-pressure 
and shock-angle measurements are sensitive to air liquefaction 
effects. Results of a thermodynamic treatment of the condensa- 
tion process are presented. Finally, some measurements on 
boundary layers and possible flow separation at high Mach 


Numbers are given. 


INTRODUCTION 


er YEARS HAVE WITNESSED a rapid increase in 
speed of supersonic aircraft, missiles, and pro- 
jectiles. In order to stay one step ahead of these de- 
velopments, research and development facilities to 
provide basic aerodynamic information on _ bodies 
moving at hypersonic Mach Numbers are required. 
One such facility is the hypersonic wind tunnel, where 
the word ‘‘hypersonic’’ is used here loosely to describe 
the Mach Number range above 5. Producing hyper- 
sonic flows suitable for model testing entails difficulties 
uncommon to conventional supersonic wind tunnels. 
To explore these new problems, a 12- by 12-cm. hyper- 
ballistics tunnel to study hypersonic flows under a 
wide range of Mach Numbers, temperatures, and 
pressures was designed and put into operation in May, 
1950, at the Naval Ordnance Laboratory. The re- 
search program presently being carried out in this 
facility primarily concerns physical problems such as 
air condensation, Mach Number determination, effect 
of temperature on tunnel operation, boundary layers, 
and instrumentation. 

Presented at the Hypersonic Aerodynamics Session, Nine 
teenth Annual Meeting, I.A.S., New York, January 29 
February 1, 1951. 

*The work presented here represents the joint effort of 
Messrs. E. Stollenwerk, S. Reed, G. Lundquist, and H. Staab 
of this laboratory. 

t Chief, Hyperballistics, Subdivision II] 


This paper summarizes some of the observations 
made so far; sample data characteristic of particular 
phenomena occurring in the tunnel are given. 


APPARATUS AND OPERATIONAL PROBLEMS 


The N.O.L. 12- by 12-cm. hyperballistics tunnel 
and some of its auxiliary components are shown in Fig. 
1. Filtered dry air (dew point —50°C. or better) is 
stored in tanks at 3,000 Ibs. per sq.in. pressure. This 
air is bled via pressure-reducing valves into a settling 
tank on top of the vertical working section. The sup- 
ply pressure of the tunnel (up to 8 atmospheres for the 
investigations to be described) is held constant during 
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Artist’s sketch of the N.O.L. 12- by 12-cm. hyper- 
ballistic wind tunnel 
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Fic. 2. Schematic expansions beyond air saturation line. 


the run. The supply air also passes through an elec- 
tric heater (SO kw.) located in the settling tank. 
Automatic power regulation controls the heating of the 
supply air at any value from room temperature to 
500°C. during the run. Then the air flows from top 
to bottom through the vertical working section into 
vacuum pumps via a 2,000 cu.m. vacuum sphere com- 
mon to all the N.O.L. tunnels. The opening of a fast- 
acting slide valve between settling tank and working 
section initiates the runs. Pressure ratios up to 50,000 
across this valve are available. The hypersonic flow 
is established in a simple nozzle with arbitrary curva- 
ture at the throat and two 50-cm. long plane diverging 
walls, ending in a square area of 12 by 12 cm. By 
means of setscrews the area ratio of this nozzle can be 
adjusted before the run to an arbitrary Mach Number 
setting 1/,. The useful operating range is 1/, = 4 to 
11. After a short parallel-wall section, the air enters 
a simple diffuser. Since supply air-mass available 
and vacuum capacity exceed the tunnel requirements, 
the operation of the tunnel is continuous. 

Some of the operational problems occurring with 
this tunnel are the following: If the tunnel operates at 
high supply temperatures, the nozzle throat gradually 
closes because of the heating of the nozzle walls. 1/, 
is extremely sensitive to small throat-area changes at 
high Mach Numbers. Before and after each run the 
throat distance is measured with feeler gages, and J/, 
is interpolated linearly with respect to time between 
these two points. Knowing the time at which a par- 
ticular measurement is made in the test section, it can 
then be approximately related to the actual .1/,. Glass 
windows are also heated gradually from the inside out, 
causing distortion of optical images. In addition, 
the schlieren system is pushed to its limit of sensitivity 
because of the low density in the test section at high 
Mach Numbers. The same is true for conventional 
pressure instrumentation, and oil, as well as McLeod, 
manometers must be used. All tunnel components 
must resist high temperatures, necessitating the use 
of copper tubing instead of rubber, silicone rubber 
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gaskets, and heat-resistant paints. Since the working 
section has no cooling system, the time history of data 
taking becomes important in contrast to supersonic 
wind-tunnel techniques. 


On AIR CONDENSATION 


Fig. 2 shows schematic wind-tunnel isentropes of air. 
Pressure and temperature drop with increasing Mach 
Number. Isentrope I of this figure intersects the vapor- 
pressure curve of air' at some point; we will call the 
Mach Number at this intersection ‘‘/,, the Mach 
Number of equilibrium condensation.’’ Here the 
danger of air liquefaction commences. An increase in 
supply pressure (isentrope II) decreases /7,, and an in- 
crease in supply temperature (isentrope III) increases 
M,.. We will call that supply temperature that 
maintains the air in the test section just at /,, or 
saturation, 7). For example, at 1, = 10 and fp = 
100 atmospheres, 7). = 700°C. Tunnel operation at 
such temperatures and pressures results in a number of 
technical difficulties. On the other hand, calculations 
based on Becker, Doering, and Volmer’s theory of 
nucleation” predicted extensive supersaturations of air, 
thus lowering the minimum supply temperature for 
successful tunnel operation at high Mach Numbers. 
This .conclusion was in agreement with experimental 
experience on water vapor in the air of supersonic wind 
tunnels, where large supercoolings can be observed. 
In fact, if the temperature gradient in the nozzle is 
high, as is the case in hypersonic tunnels, water vapor 
could be supercooled up to 100°C. before condensing 
spontaneously.* However, experiments with air do 
not show this supersaturation. C.H. McLellan demon- 
strated first that air condensation in hypersonic tun- 
nels begins near the saturation point.‘ This fact was 
also borne out by the investigation described here, as 
well as those by Stever® and by Nagamatsu.® 


PRESSURE MEASUREMENTS WITH AND WITHOUT AIR 
CONDENSATION 


Table 1 gives a collection of data for different ./, 
measured at the nozzle exit at supply temperatures of 
the order of room temperature, called ‘“‘cold runs.” 
Pitot pressure (fo’) and static pressure (p) are meas- 
ured and divided by the supply pressure (fo). The 
corresponding Mach Numbers are found in an isen- 
tropic flow table;’ they are also determined directly 
from the Rayleigh pitot-tube formula’ (/o’/p). These 
measured ‘“‘cold run’? Mach Numbers are in disagree 
ment. The pitot Mach Number is somewhat lower 
than .1/,, as expected from the presence of the boundary 
layer. However, the static-pressure Mach Numbers 
and, in particular, those from the Rayleigh formula are 
considerably lower. This conventional evaluation 
assumes the isentropic expansion of a perfect gas in the 
nozzle, an assumption that is not true if a part of the 
air liquefies. Therefore, the discrepancy in formal .\/ 
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evaluations from pitot and static pressure measure- 
ments may be regarded as an arbitrary measure of the 
magnitude of the air-condensation effects. The latter 
increase with increasing .1/,. These cold-run data are 
in agreement with those given by Nagamatsu for the 
C.L.T, hypersonic tunnel, which is nearly the same size 
as the one at N.O.L.°® 

Fig. 3 shows dimensionless static pressure measured 
at various distances on the centerline of the nozzle wall. 
Such data are given for a ‘“‘cold run” and for a ‘“‘hot 
run.” In the latter case the supply temperature was 
raised to the order 7%... Cold-run and hot-run static 
pressures diverge in the vicinity of the air saturation 
point in the nozzle of the cold run, with the hot-run 
static pressures being lower at the nozzle end. The tri- 
angles in Fig. 3 show the results for a cold run, where 
the ratio of supply density to atmospheric density 
(po/Patm.) iS equal to that for the hot run in order to 
minimize Reynolds effects. (Reynolds Number itself 
is a dubious parameter if liquefaction sets in.) Be- 
cause both the low- and high-density cold-run static 
pressures are nearly equal, the discrepancy between 
cold runs and hot runs cannot be attributed to a dif- 
ference in boundary-layer thicknesses for the cold and 
hot case. 

Fig. 4 shows “‘formal’’ Mach Number determinations 
by use of isentropic-flow tables from pitot and static 
pressure measurements at the nozzle exit as a function 
of supply temperature. Four different Mach Number 
settings are shown. Only when 7) nears 7%, does 
agreement of the Mach Number derivation from pitot 
and static pressure measurements become apparent. 
The pitot pressure is nearly insensitive to liquefaction 
effects throughout the range of test temperatures, 
while static pressure changes considerably with supply 
temperature, as shown by the derived 7 on Fig. 4. 
The change of J/, during the runs due to throat closure 
by heating and the saturation line of air are also indi- 


cated on Fig. 4. 
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Formal Mach Number Determinations from Pressure Measure- 
ments at Nozzle Exit fer Low Supply Temperature 
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Fic. 3. Static pressure measurements on nozzle-wall centerline 


THERMODYNAMICS OF THE CONDENSATION PROCESS 


Scrutiny of Fig. 3 shows that the expansion with air 
condensation is characterized by a gradual deviation 
of static pressure from that obtained for the hot run 
without condensation. It can be shown that this 
deviation begins about the time that J/, is reached in 
the nozzle. The exact location of beginning condensa- 
tion is difficult to determine because of the inaccuracy 
At J,, flow pres- 
sure and temperature are equal to equilibrium satura- 


of static pressure measurements. 


tion pressure and temperature of air. Our seemingly 
gradual air condensation would be unlikely to occur 
without a sufficient number of condensation nuclei 
available at or near the air saturation point. These 
nuclei provide surfaces for the air to condense on. It 
was estimated that these nuclei are, in fact, available 
in large numbers because of the previous condensation 
of the rest of unavoidable water vapor and of carbon 
dioxide in the air, upstream of the air saturation point. 
With this assumption we suppose that the condensation 
process is the following: The air follows a perfect-gas 
isentrope up to the saturation line. The air then 
enters the region of coexistence of liquid and vapor 
phase and follows there an isentropic saturated ex 


pansion. With the knowledge of the vapor-pressure 
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Fic. 4. Mach Number determined from pitot and static 
pressure measurements at nozzle exit as a function of wind-tunnel 
supply temperature 
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Fic. 5. Preliminary experimental criteria for supply tempera- 
tures to minimize air-condensation effects in the test section 


equation, equation of state, and experimental data on 
the static pressure with condensation (cold runs in 
Fig. 3), all other flow quantities of interest can be 
computed by use of the thermodynamics of flow sys- 
tems. The list below summarizes the effects of con- 
densation on flow parameters. All parameters are 
made dimensionless relative to conditions upstream of 
the saturation point—e.g., conditions in the settling 
tank. Each parameter for condensed flow is then 
compared to that obtained at the same M, with the 
supply temperature raised to 70, so that no condensa- 
tion occurs. Then one finds, for the dimensionless 
parameters, that: 

(a) Pitot pressure remains nearly unchanged. 

(b) Static pressure does not drop to low values. 

(c) Average density (liquid plus vapor) remains 
nearly unchanged. 

(d) Temperature does not drop to low values. 

(e) Velocity remains nearly unchanged. 
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(f) Actual-flow Mach Number does not increase to 
high values. 

(g) The amount of air liquefied at the nozzle end js 
about 10 to 20 per cent of the total air mass for the 
N.O.L. pressures, temperatures, and Mach Numbers, 

(h) Evaluating 7 from isentropic-flow tables is a 
formal procedure leading to incorrect results. 

The above list provides a tentative guide to the selec- 
tion of proper methods of investigation of condensation 
effects. Quantitative data, computed from this ther- 
modynamic condensation model, however, depend in 
their accuracy on the measured expansion and, in 
particular, on available vapor-pressure data.' The 
latter are obtained by extrapolation from measure- 
ments made of pressures and temperatures that are 
higher than those found in the wind tunnel, thus intro- 
ducing considerable uncertainty.° 


PRELIMINARY CRITERIA FOR FLOWS PRACTICALLY 
FREE OF AIR CONDENSATION 


Examination of Fig. 4 shows that general agreement 
of Mach Numbers derived from pitot pressure and 
static pressure measurements with an isentropic-flow 
table occurs before 7) reaches 7). An investigation 
of this apparent agreement reveals four factors. 
(1) Any increase in 7) increases A/, (Fig. 2). (2) If 
the expansion is started at larger specific volumes, as 
is the case for higher 75, inspection of a p-v diagram 
shows that the rate of liquefaction during the expansion 
is smaller." (8) The existing measuring accuracy 
does not permit detection of small amounts (order of 5 
per cent or less) of liquid air in the flow. Therefore, 
it is assumed that supply temperatures sufficient to 
result in agreement of Mach Number determinations 
by different measurement methods on the basis of 
perfect-gas expansions, are sufficient to permit survey 
model testing. (4) Because of the inability to locate 
the beginning of condensation exactly, some super- 
saturation is not impossible. The resulting saving in 
supply temperature is shown in Fig. 5 which is ap- 
plicable only to a tunnel of the size of the N.O.L. 
tunnel and a measuring accuracy and fo-range like that 
of our tests. The supply temperatures that minimize 
air condensation effects are then about 100°C. below 


Te: 
OBSERVATIONS ON SIMPLE BODIES 


Only small changes in Mach angle and shock-wave 
angle on bodies take place if the Mach Number of the 
flow changes from, say, 8 to 10. In general, the shock- 
wave pattern around bodies is rather insensitive to M 
changes in the hypersonic range. However, qualita- 
tive data can be obtained from shock angles as shown 
on Fig. 6. Three schlieren pictures of a i0° wedge 
are shown at M, = 7. The supply temperature 
ranges from low temperatures to 7»., while the supply 


density is kept constant to minimize changes of the 
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INVESTIGATIONS IN 
refractive index in the test section. The increase in 
M due to elimination of condensation by preheating 
can now be seen directly as a decrease in shock angle 
on the 10° wedge by noting the wave position with 
respect to the (accidental) bubble in the glass wall. 
The top curve of Fig. 7 gives these shock angles de- 
termined by two observers for both sides of the wedge 
from several such pictures. No further decrease in 
shock angle or corresponding increase in Mach Num- 
above about 7) = 150°C.—the 
temperature sufficient to minimize 
liquefaction effects at this J/, (Fig. 5). The Mach 
Numbers determined from the shock angle for a theo- 
retical 5° flow deflection'! at high supply temperature 
without liquefaction are lower than those given for 
This is due to 


ber can be noted 


that was termed 


example by static pressure (Fig. 4). 
the boundary layer on the wedge which increases the 
“effective’ flow deflection angle. The lower curve 
in Fig. 7 presents an attempt to eliminate this bound- 
ary-layer influence. The pressure is 
measured by one pressure tap on each side, halfway 
Since the static 


wedge-surface 


between tip and end of the wedge. 
pressure ahead of the wedge is also known, the pres- 
sure ratio across the oblique shock can be determined. 
From the relationship of this pressure ratio to Mach 
Number and shock angle for plane flow deflections,* 
M can be determined without knowing the deflection 
angle. Since both wedge sides show systematically 
different shock angles, a '/3;° correction is applied to 
compensate for either flow inclination or misalignment 
of the model. Again, this is a “formal’’ evaluation 
(bottom, Fig. 7) which is correct only at the high sup- 
ply temperatures. That the Mach Number at high 
supply temperatures is still lower than that found by 
static pressure (dashed line) may be attributed to the 
fact that the pressure on other points of the wedge 
surface is different from that measured at the center- 
This fact points out the need of detailed studies 
of viscous effects at hypersonic Mach Numbers. 
Shock angles on blunter bodies (above about 40° nose 
angle) were found to be rather insensitive to condensa- 
tion or M changes. In particular, no difference in the 
shock pattern about spheres could be noted with or 
It is suspected that strong 
pressure 


point. 


without condensation. 


shocks with their considerable increases in 
and temperature cause a vaporization of some of the 
liquefied air. This phenomenon has previously been 
observed for condensed water vapor by the method of 
light scattering by Head.'* In the N.O.L. air case, 
light scattering is also visible if a strong beam of col- 
limated light traverses the condensed flow. This light 
scattering disappears behind strong shocks on bodies, 
as well as in the boundary layer. (It also disappears 
if the air is heated.) 

Fig. 8 shows a 40° cone-cylinder body at M = 8.25. 
Air-condensation effects are negligible at 7) = 300°C. 
(Figs. 4 and 5), and this is the highest Mach Number 
at N.O.L. The schlieren system 


achieved to date 
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Fic. 6. Schlieren picture of head shock wave on a 10° wedge as a 
function of supply temperature 
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Schlieren picture of 40° cone-cylinder body at J = 8.25 
(po) = 7.8 atmospheres, 77 390°C 
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Fic. 9. 
nozzle-exit center for flow with and without air condensation 


shows only the head wave (static pressure in the test 
section is about '/. mm. mercury), and the picture is 
further obscured by air convection on the outside of the 
“hot” tunnel. New methods of flow visualization 
must be developed for such cases. However, it can 
be seen that all changes in flow conditions due to the 
presence of the body are confined to a narrow zone near 


the body. 


OBSERVATIONS OF BOUNDARY LAYERS AND FLOW 
SEPARATION 


Fig. 9 shows two pitot-pressure profiles taken from 
nozzle wall to centerline at right angles to the flow at 
the nozzle exit center. In both cases the supply den- 
sity is the same. One case is with, the other without, 
air condensation. That condensation does not affect 
pitot pressure in the center of the flow was shown on 
However, in the boundary layer the pitot 


On the other 


Fig. 4. 
profiles also show no marked difference. 
hand, tentative stagnation temperature profiles, meas- 
ured with a shielded total head temperature probe 
traversing the flow like the pitot tube, are different 


with and without liquefaction. For hot runs, the ratio 


Pitot-pressure profiles at right angle to nozzle wall at the 


SCIENCES—OCTOBER, 1951 

of stagnation to supply temperature reaches a maximun 
at the edge of the boundary layer, as also observed by 
McLellan at about the same Mach Number." With 
condensation, this maximum of the ratio reaches 
higher values and is moved into the middle of the 
boundary layer. Before quantitative data can _ be 
given, more development and calibration work on such 
shielded probes is required, since they must operate 
free of errors because of radiation in the immediat: 
vicinity of the heated nozzle walls. 

Some observations tentatively identified as flow 
separation might be of further interest. Fig. 10 shows 
boundary-layer profiles taken as the ones shown be 
fore. Only hot runs (to eliminate condensation) are 
given for various Reynolds Numbers, where Reynolds 
Number is based on centerline free-stream conditions 
at the nozzle exit for the 12-cm. width of the nozzle. 
With decreasing Reynolds Number the flow apparently 
detaches from the nozzle wall, as evidenced by the in 
The effect can 
On the 
nozzle-wall side of the profiles near the maximum, the 


creased peak size in the pitot profiles. 
also be observed with the schlieren system. 


pitot pressure fluctuates pointing to an unsteady back 
and-forth motion of the separation point. The same 


(Continued on page 682 
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Skin-Friction and Heat-Transfer 
Characteristics of a Laminar Boundary 
Layer on a Cylinder in Axial Incompressible 


Flow 


R. A. SEBAN# ano R. BOND? 
University of California at Berkeley 


ABSTRACT 


A solution is given for the case of the laminar boundary layer 
of an incompressible fluid of constant properties on the exterior 
of a cylinder with flow parallel to the cylinder axis. This case 
differs from the Blasius solution for flow along a flat plate by 
considering the effect of the curvature in a plane transverse to 
the flow direction. The local skin-friction heat-transfer 
coefficients for a Prandtl Number of 0.715 are 
compared to the similar magnitudes for flat plate flow, and the 


and 
evaluated and 


effect of the curvature is shown to be significant in some practical 


cases. Recovery factors are evaluated, and this quantity is 


found to be insensitive to the effect of curvature of the boundary 


NOMENCLATURE 


= cylinder radius, ft 


? = heat capacity at constant pressure, B.t.u slug. F. 
H = enthalpy, B.t.u./slug. 
k = thermal conductivity, B.t.u./(hr. °F. ft.) 
= heat flux, B.t.u./(ft.? hr.) 
= radial coordinate, ft 
] = temperature, °F. 
u = axial velocity component, ft. per sec. 
= radial velocity component, ft. per sec. 
x = axial coordinate, ft. 
! = kinematic viscosity, sq.ft. per sec. 
p = fluid density, slugs per cu-ft. 
a = Prandtl Number, uc,/k 
T = shear stress, lbs. per sq.ft 
im = absolute viscosity, lb. sec. /ft.? 
y = stream function 
é = modulus giving relative thickness of boundary layer, 
4(x/a)V/ v/uyx 
ugx /v = Reynolds Number 
cy = skin-friction coefficient 


0 = free-stream condition 
w = wall condition 
F.P. = flat plate value for same Reynolds Number 
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INTRODUCTION 


WwW THE INCREASED FREQUENCY of prototype 
flight and model tests under low-pressure condi- 
tions, the characteristics of a boundary layer, the 
significant dimension of which is comparable to the 
model dimension, have become of significance. It 
can be noted superficially that, for the case of a bound- 
ary-layer flow over a slender body, the momentum and 
energy rates within the boundary layer will be in- 
creased over the corresponding magnitudes for a flat 
surface because of the curvature of the boundary 
layer in a plane transverse to the flow. If the bound- 
ary-layer thickness becomes comparable to the thick- 
ness of the body, it is to be expected that this effect 
will become significant. 

The simplest model upon which an analytical in- 
vestigation of this effect can be carried out is that of a 
laminar boundary layer of a fluid of constant properties 
on a slender circular cylinder with axis parallel to the 
flow. To make the analysis possible, it is necessary 
to postulate the beginning of the boundary layer with 
zero thickness at the leading edge of the cylinder. 
While this is an abnormal requirement, assuming a 
suction flow through the interior of the cylinder of such 
magnitude as to make the stagnation stream lines of 
the flow coincident with the cylinder generators, it will 
be noted later that the applicability of the solution 
may not be essentially restricted by the necessity of 
this assumption. 

This problem, which differs from the Blasius solu 
tion for a laminar boundary layer on a flat plate only 
in the of the to the 
direction of flow, has been noted by Young,’ Jakob 
Their considerations involved 


curvature surface transverse 
and Dow,’ and others. 
the assumption that the velocity distribution within 
the boundary layer and the growth of the boundary 
layer upon the cylindrical surface were the same as 
upon a flat plate. Application of the integral form of 


the momentum equation under these assumptions 
yielded an expression for the local skin-friction co 
efficient, the magnitude thereof being greater than 
that for the flat plate because of the proportionately 
the boundary layer on the 


that the skin 


greater cross section of 


curved surface. Young found local 
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friction coefficient for the cylindrical surface could be 
related to that for the flat plate for the same length 
Reynolds Number by the expression 


Cy/ Copp, =] 1 g 
where 
& = 4(x/a) V y/ UX 


In the present investigation the problem is considered 
more exactly, without requiring assumptions of the 
type used by Young. In addition, the solution of the 
momentum equation has been applied to the energy 
equation to specify the rate of heat transfer and the 
effect of frictional heating for a fluid such as air, 
having a Prandtl Number of 0.715. 


ANALYSIS 


For the case of an incompressible fluid of constant 
properties flowing with uniform free-stream velocity 
parallel to the axis of a right circular cylinder, the 
boundary-layer forms of the equations of continuity, 
momentum, and energy can be written as follows: 


Ou 1 O(rv) 
+ 


= 0 1) 
ox r Or ( 
[Ou Ou 1 O(rr) 
fu ot 4 pH) () 
106 or r Or 
| oH ea oT} ou n 1 O(rq) 9) 
2 hy Ox =| ai or , ( 
with 
T = p(Ou/Or) (4) 
g= k(OH/or) (5) 


Eq. (1) is satisfied by the introduction of a stream 


function, such that 
u = (1/r)(OW/Or) (6) 
v = —(1/r)(OW/Ox) (7) 


The remaining equations, Eqs. (2) to (5), are trans- 
formed by a change of variables—-the independent 
variables (x, r) to (&, 7) and the dependent variables 


(vy, H) to (f, g). 


The transformation is defined by 


—& = (4/a)V vx/uUo (8) 

n = Viuo/vx[(r? — a?)/4a] (9) 
y = av vxumf (10) 

H = Hy + wg (11) 


Eqs. (6) and (7) then become 


u = (t/2)(Of/On) (12) 
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la y" ( of ) Qv df 
v= elf le (13) 
2¢ x "Oe . r O€ « 
Eqs. (2) to (5), for momentum and energy, become 


or orf orf 
Dal Ct + mf) nf +s55 + 
mk. 


On? On? 
c of of oF) “i 
QE On? dn dyve/ C4) 


Lo Og Og 
=f +e iain 


1 st) (= Og Of 22) 
1 + n€ ~==)=0 (15 
4 (1 + nf) ee + €& dE On On dE QO (15) 


; - 
lr Itty O7 
{i— f (16) 


T = ~ —pMo i 
4a V vx On? 
lrk |. luo Og 17 
- uo” VJ (17) 
2aC, \ vx On 
The boundary conditions considered are u = 0 and 


v = Oforr = aand uw = m~ and A = Ay forr = o 


or x = 0. In the new variables these boundary condi- 


tions become 


f=0, =—=0, =—=0 (ory = 0) (18) 


Of/On = 2, g=0 (forn7 = ~) (19) 


Taking the wall value in Eq. (16), the local skin- 
friction coefficient is 


—— 1 |» ed " 
Cr = = ZU) 
. (p/2)uo” 2View On’/,=0 ‘ 


A boundary-layer displacement thickness is defined 


ae We .. 
6, = ( — ) dr = \ lim (2n — f) 
Ja Up/ a U0 ya 


(21) 


as 


The method of solution adopted for the determina- 
tion of f and g from Eqs. (14) and (15) is based on 
expansion in a series of functions of 7 with each function 


multiplied by a power of &. The functions thus become 


S(E, n) = foln) + Efi(n) + Efo(n)... (22) 
g(é, 7) = go(n) + Egi(n) +... (23) 


Substitution of these functions in Eqs. (14) and 
(15) and equation of the coefficients of like powers of 
& to zero yield a series of differential equations as 


follows: 
f0'? tf. ff! = f (24 


dd at fof" een he Ay aa Df," f, 4 “_ 2 f te =e (25) 
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fl”! + fof2”” — 2fo'fe’ + 3fo' "fe + Eq. (24) is the well-known Blasius equation, and 

nh’ + fi’ + 2ffi'’ — fi = 0 (26) its solution has been taken from the calculations of 
Cope and Hartree.* Approximate solutions of Eqs. 
(25) and (26) were then obtained by numerical inte- 


/ ” Pa # /A\ £142 97) 
(1/a)go" + Jogo’ + (1/4)fo* = 0 aad gration using first-order approximations with finite 
7 oh «tas i 4 gy!) + intervals of 0.1 in the independent variable. Curves 
7 Jo&i ~ Jo 8 T \nBo 80 of the functions f;’ and f2’ so determined are shown 
on Fig. 1, and these functions specify the velocity 


2figo’ — 4 fo"? + Sfo"’fi” = 0 (28) distribution. The most important results are in the 
values of the functions at the extreme values of 7. 


The boundary conditions associated with fo, fi, fe These are 


™ fo’’(0) = 1.328, fi’’(0) = 0.704, fo’’(0) = —0.095 (29) 
fo(0) = fo'(0) =O fo’ () = 2 
fi(0) = fx'(0) = 0 fi'(@) =0 lim (2n — fo) = 1.721, fa(@) = 0.051, 
=0 ver fol) = —0.093 (30) 


f(0) = fe'(0) = 0 fr'(~) = 
aes ars The availability of the functions f; and f2 enables 
For the functions go, gi, ..., the only prescribed the solution of the energy equations, Eqs. (27) and 
(28), in a similar way. This has been done for a 
Prandtl Number of 0.715, leading to the following 
results: 


boundary conditions are 
go(°) = 0, 21(7) = 90,... 


while the values of the function and its derivatives : 
. ; 29° (O) = 0.249 — 0.590g0(0) 31) 
at 7 = 0 are adjustable to accommodate a choice of go (0) D0 go \ 


thermal boundary conditions at the wall. g:'(0) = 0.1457 — 0.339g.(0) — 0.7824g,(0) (32) 


The skin-friction coefficient and boundary-layer displacement thickness as defined by Eqs. (20) and (21) can 
be evaluated from the known functions. The skin-friction coefficient is 
Cy = 0.664V v/xuo(1 + 0.53 — 0.071£? .. .) (33) 
The first factor on the right of Eq. (33) is the skin-friction coefficient for a flat plate as given by the Blasius 
solution. Designating this by c,,,, the result given by Eq. (33) can also be written as 


C)/Crpp = 1 + 0.53E — 0.0718? ... 


The displacement thickness is 


61 = V/ vx Uo lim (2n - fo) = Efi(@) on E°f2( @ ) os ]| 


6: = 1.721/vx/uo(1 — 0.029F + 0.0548? .. .) | 


(34) 


and this can be written in reference to flat plate values as 
51/61 pp = 1 — 0.029F + 0.0542 ... 


The expressions of Eqs. (33) and (34) are approximations permissible only if — is small. An estimate of this 
approximation can be gained from Fig. 2, where the representation of Eq. (33) reveals that the use of the first 
two terms of that equation expresses the skin-friction coefficient within 10 per cent at values of £ large enough 
to make c,/c;,, equal to 1.50. 

The boundary-layer thickness as given by Eq. (34) is seen to be initially less than that for a flat plate. 
subsequent behavior is not clearly defined because of the poor convergence and limited number of terms in the 


Its 
series as obtained. 

The heat-transfer rate at the wall is given by Eq. (17) and the results as given by Eqs. (31) and (32). This is 

gq = (1/2)(R/c,)Uo?V Uo/ vx} 0.249 — 0.59g0(0) + [0.1457 — 0.34g(0) — 0.782g¢,(0) JE + . } (35) 


This result is unwieldy and simple results are obtained only for some special cases. 
For an insulated wall, g = 0 for all £, and two equations result from Eq. (35), from which go(0) and g;(0) can 


be calculated. Then 
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g(t, 0) = 0.422 + 0.0034€ ... (36 
and the surface temperature is (with /7 = c,7) 
T(E, 0) = To + (uo?/2c,)(0.844 + 0.007 .. .) (37) 

The last terms, which give the recovery factor, are seen to depend but slightly on é and are practically equivalent 
to the recovery factor of 0.844 resulting from the Pohlhausen solution for the flat plate for the Prandtl Number 
of 0.715 as considered. 

For the case of low-velocity flow, in which the frictional heating effect may be neglected, Eq. (35) assumes the 
simpler form 


q = (1/2)(k/cp)to?V uo/ vx | —0.59g0(0) + [—0.34g0(0) — 0.782¢:(0) JE... | (38) 
If, in addition, the temperature is uniform along the wall, then 
go(O) = (Cy Uo?) (Tv —_ To) 


and g;(0) must be zero. A heat-transfer coefficient can then be defined as 


h = —q/(Te — To) = 0.295kV/Uo/vx(1 + 0.575 +... (39) 
and the Nusselt Number is given by the equation 
hx/k = 0.295+~/xuo/v(1 + O.575E + ...) (40) 


In this equation the first factor on the right side is the expression of the Nusselt Number for the case of flow over 
an isothermal flat plate. 

If the velocity is large enough so that frictional heating must be taken into account, it is not possible to express 
rigorously the heat-transfer coefficient in terms of an effective temperature difference. For the case of uniform 
wall temperature, g:(0) = 0, and the heat rate is given as 


g = (1/2)(k Cy) Uo? Uy vx[—0.59(1 + 0.575€) (c/o?) (Ti, — Ty) + 0.24911 + 0.5858) 
This is, approximately, 


luo ae ion Wu 
q=k —0.295(T,, — To) + 0.12455 


\ Vx 2b» 


9 
07 


9 
07 


(1 + 0.589) 


<Cp 


| 
a jUol _. a u 
= —().295k T, — (7, + 0.8445 


\ vx 


(1 + 0580 


Since {7 + O.S44(a%?/2c,)| is approximately the adiabatic wall temperature, the definition of the effective 
temperature difference in the usual way, as the difference between the wall temperature and the adiabatic wall 
temperature, will give, from Eq. (41), a heat-transfer coefficient and Nusselt Number similar to that given by Eqs. 
(39) and (40) for the case of nondissipative flow over an isothermal plate. 

It is of further interest to note that the combination of Eqs. (40) and (33) yields the expression 

hx /k cpl + 0.575€ 
=! “sg = 0.715 (42) 
(Mx v)a( I 3) 2? 1 + 0.58E 

It is seen from this expression that the ratio of the heat-transfer to skin-friction coefficients is but slightly 

affected by the curvature of the wall. 


DISCUSSION frictional heating. In the case of the energy equation, 


a . ? , , , the results, given for o = 0.715, were not carried out to 
rhe results of the present theoretical investigation ; ; : 

. : . so good an approximation as for the momentum 
reveal that, while the effect of curvature on the skin- 


eee iP : ' equation, but they are of sufficient accuracy up to the 
friction coefficient for the case of laminar boundary- é : : : Bes os 

; : point at which the heat-transfer coefficient is increased, 

layer flow parallel to a cylinder is not so great as : ae ; 

: 5 by the effect of curvature, to one and one-half times its 


estimated by Young, it is still significant in many 
: 6 6 ~ flat plate value. 


possible cases. The solution has been extended to 
give wall temperatures and heat rates for the cases of The effects of curvature on the skin-friction and 
an adiabatic and isothermal wall, with and without heat-transfer coefficients are similar, so that the ratio 
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Fic. 1. Functions defining the velocity distribution The 
function fy'(m) defines the flat plate velocity distribution, and the 
velocity distribution on the cylinder is given by 























u/uy = */o[fo'(n) + Efi'(n) + E%fe"(n) 
2.0) 
ae 
5+ a 
Cr 
— 6 
Ctep 
as; 
% 05 1.0 15 
4% Vv 


a Vuox 
Fic. 2. Ratio of local skin friction on the cylinder to that on 
plate at same Reynolds Number 
(a) ¢s/cy,, = 1 + 0.538¢ 
=1+4 0.53 


E — ().071£2 
(b) cy Sep 3s 
of these coefficients remains almost unaltered from its 
flat plate value. The recovery factor is relatively 
unaffected by the curvature, and this lack of sensitivity 
is probably due to the almost equal increase in the 
heat-transfer and skin-friction coefficients. 

The relative magnitude of the effect of curvature 
can be demonstrated by the consideration of some 
typical cases. If the circumstances under which the 
local skin friction is increased by 10 per cent over the 
value that would obtain on a flat plate at the same 
Reynolds Number are considered, then at this point, 
from Eq. (33), € = 0.19 and 


x @19 UoX _ |uox 
= = 0.05 | 


a t\Vv Vv 
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or 
x/a = 25 X 1074(uga/ v) 
For a |-in. diameter cylinder, with uw = 100 ft. per 
sec. and vy = 1.8 & 10~‘ sq.ft. per sec., typical of air 
conditions near ground level, this gives x/a = 4580, 


so that for most practical cases the effect would be 
insignificant. If, however, with all other conditions 
remaining the same, the pressure is reduced to 1/100 
of an atmosphere, the skin-friction coefficient (and the 
heat-transfer coefficient) would be 10 per cent above 
the flat plate magnitude at x/a = 5.80. 

For a cylinder 4 ft. in diameter, in flight under 
standard atmosphere conditions, the local skin-friction 
coefficient exceeds the flat plate magnitude, c, = 
1.10cyp p, at a point 20 ft. from the leading edge when 
the Mach Number is 1.3 at an altitude of 200,000 ft. 
and when the Mach Number is 3.5 at 240,000 ft. 

It was noted initially that the determination of the 
solution required the assumption of zero boundary 
layer thickness at the leading edge of the cylinder, 
a condition impossible to fulfillment in practical 
cases. While little is known in general of the develop- 
ment of a laminar boundary layer on an arbitrary 
body of revolution, the development of a laminar 
boundary layer on a half body has been studied by 
Scholkemeier.t- He demonstrated that on this body 
the boundary layer assumes the characteristics of a 
boundary layer of a flat plate at a distance of 1'/; 
diameters from the stagnation point, the Reynolds 
Number for a flat plate case being taken on the basis 
of the meridianal distance from the stagnation point. 
The similarity of the half body to a cylinder with a 
hemispherical nose suggests that this type of evalua- 
tion, combined with the methods outlined herein, will 
correctly establish the boundary-layer characteristics 
on a cylinder of this type. 
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Analysis 


of a Friction Damper for 


Clutch- Type Servomechanisms 


HERBERT K. WEISS* 
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SUMMARY 


The stabilizing action of a flywheel, coupled by coulomb fric- 
tion to the output shaft of a clutch-type servomechanism is 
analyzed. It is shown that such a damping device is capable of 
partly compensating for the undesirable effects of time lag in the 
control. Methods are presented for quickly determining the 
transient and steady state of both servomechanism and damper 
The analysis of the damper-servo 
is accomplished by 


for any set of initial conditions. 
system, which has two degrees of freedom, 
tracing pairs of trajectories in the phase plane. 
for a given design of damper, the amplitude of the steady-state 
oscillation is proportional to the square of time lag in control func- 
tioning, while the period of the oscillation is proportional to the 


It is shown that 


time lag. 

Graphs are presented to show the optimum ratio of friction 
torque to clutch torque for minimum amplitude of steady-state 
oscillation. The greater the flywheel inertia, the smaller the 
amplitude of steady-state oscillation can be made for a given 
Such damping is always obtained, how- 
although 


time lag in the control. 
ever, at the expense of maximum servo acceleration, 
maximum servo velocity is not reduced by the friction-coupled 


flywheel. 


INTRODUCTION 


QO”™ OF THE SIMPLEST and most useful types of 
damping devices for suppressing hunting of 
servomechanism is a flywheel, coupled to the output 
shaft of the servo by coulomb friction. Known as a 
“Lanchester damper,’ such a device is also used 

a vibration damper."? The analysis offered earlier’ 
for the vibration damper is not, however, directly ap- 
plicable to the servomechanism case, since it assumes 
that the motion of the oscillating shaft is a sine wave, 
regardless of the damper’s action. 


When coupled to the usual torque-producing element 
of a servomechanism, the friction damper affects the 
manner that linear methods of 
In the present 


response in such a 
analysis are no longer strictly applicable. 
study, advantage is taken of phase-plane representation 
to obtain analyses of the transient and steady-state 
motion of a clutch-type servomechanism coupled by 
coulomb friction to an inertia element. This approach 
permits a graphical solution to be readily obtained for 
the motion, even though nonlinearities are introduced 
driving and damping 


by the abrupt changes in 


torques. 
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CLUTCH-TYPE SERVOMECHANISM 


Descriptions of a clutch-type servo are available in 
the literature.*4 This type of servomechanism has 
been found useful in applications where the position of 
a shaft is to be determined accurately but where high- 
frequency oscillations of small amplitude about the 
desired position are not undesirable. A constant-speed 
motor constitutes a source of power to which the output 
shaft is coupled by friction clutches for rotation in 
either direction. When a clutch is closed, it is assumed 
that a constant torque is transmitted across the clutch 
faces until the output comes up to drive speed—i.e,, 
until relative motion between the clutch faces ceases. 
For present purposes, the simplest representation of 
the servo elements will be employed, as indicated in 
Fig. 1, according to the symbolism originated by Hazen.* 
It is assumed that when servo output differs from in- 
put, the proper clutch is closed to apply a torque in a 
direction tending to reduce the difference to zero. 
Constant delay time ¢, in the action of the clutch is 
assumed. The load is assumed to consist of inertia 
alone, the only damping being supplied by the friction- 
coupled flywheel. Finite size of the inactive zone be- 
tween the clutch-actuating contacts is neglected. The 
frictional torque exerted across the friction coupling is 
assumed to be independent of the magnitude of relative 
velocity across the coupling faces and to depend only 
upon the sign of the relative velocity. 

The methods to be used are an adaptation of those 
developed in an earlier study of relay servomechanisms 
by the author.® 


SYMBOLS 


The parameters associated with the components of the servo- 
mechanism in Fig. 1 are identified by the following symbols. 
Im = 


moment of inertia of the servo plus load, referred to the 
output shaft 


Ig = moment of inertia of the flywheel damper, referred to 
the output shaft 

Tm = torque applied by the driving clutch, referred to the 
output shaft 

Ta = friction torque applied by the damper, referred to the 
output shaft 

0m = angular position of servo output shaft 

64 = angular position of damper 

6; = angular position of servo input shaft 

t = time 

la = constant time delay in clutch operation 
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Fic. 1. Schematic of clutch-type servomechanism with friction 


damper 


The following parameters are derived from those above. 


r; = Ig/Im = ratio of moments of inertia of damper 

to servo plus load (1) 
r = |Ta/Tm\| = ratio of damping torque to driving 

torque (2) 
ém = Om — 0; = servo error (3) 
ba = 0a — % (4) 
a = |Tm/Im| = servo acceleration without damper (5) 
r = time required for the servo without damper to 

move !/, radian from a rest position—ti.e., 1 = 

2ar? (6) 
s = t/r = dimensionless measure of elapsed time (7) 
Sa = ta/t dimensionless measure of delay time (8) 


EQUATIONS OF MOTION 
Summing torques on the servo output shaft 


Bg — 0) — Ty sgn (Omn’ — 0a’) (9) 


= —T,, sen (Om 
where ‘‘sgn’’ means “‘sign of.'’ Assuming constant in- 
put velocity 0; and using Eqs. (2), (3), and (4), Eq. 
(9) is reduced to 


6," = —asgn 6, — r, a sgn (6m’ — 64’) (10) 


For the damper 


Ta" = V4 sgn (Am’ “2 64’) (11) 
whence 
ba” = (ri/ria sgn (bm! — 5a’) (12) 


The time scale is now changed from ¢ to s by intro- 
ducing Eq. (7) and noting that 


¢, = d%b,,/ds? = 6,,"(ds/dt)* = 6,,"r* (13) 


Similarly 
Gq = 64"7? (14) 

and also 
* N , * , / ~ 
Tm = OmT; t = 6a T (15) 
Sa = b,: Ge = &, (16) 


Then the equations of motion for servo and damper 
become 


26m = —SEN Om — 7, SEN (Om — Ga) (17) 


2¢4 = (7./rs) sgn (Gm — Ga) (18) 
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S.Lip-STICK CRITERION 


Consider the behavior of the servo with damper, 
wlien both servo and damper are initially at rest and 
a control contact is suddenly closed. The driving 
torque 7, is suddenly applied. A maximum torque 
T, can be transferred across the friction coupling to 
the damper flywheel. Hence, the maximum acceler- 
ation of the flywheel is 7,/Jg. If this is greater than 
the servo acceleration 6,,, the damper “‘sticks’’ to the 
servo, which acts as if it had a total inertia Jg + J». 
In this case there is no damping action. In order to 
obtain damping, the damper flywheel must slip, rela- 
It is therefore required, for 
But if the damper slips 


tive to the servo shaft. 
slipping, that 7y/I¢ < 8n. 

6m = (Im — Ta)/Lm (19) 
whence the requirement that the damper slip is estab- 
lished by 


Ta/Ta < (Tm — Ta)/Lm (20) 
or 
S < 1.00 (21) 
where 
S = 7,/(r{1 — r,)] (22) 


For the simple servo assumed, if the damper is ad- 
justed according to Eq. (21) so that it slips, it slips as 
long as the servo is in motion. 


PHASE PLANE AND TRAJECTORIES 


From Eqs. (17) and (18), the motion of damper and 
servo can be solved simply as long as the signs of ¢,, and 
(Gm — Ga) do not change. By plotting the “‘trajecto- 
ries’’ of the motion of servo and damper in the “‘phase,”’ 
or position-velocity plane, the changes in sign of these 
quantities can be easily monitored. The phase-plane 
approach has been extensively studied by Minorsky’ 
for problems in nonlinear mechanics, and its applica- 
tion to the analysis of servomechanisms has been made 
earlier by the writer,® and others.* 

In the present problem, the motion of two relatively 
moving shafts—the shafts of servo and damper—imust 
be simultaneously monitored in order to observe the 
sign of the difference in velocities. This is accomplished 
by plotting both c,, and o, against o,,, the latter serving 
in place of a coordinating time reference. Eqs. (17) 
and (18) may be written 


Gm = k,/2 (23) 


G4 = k 2 (24) 


where k;, k, depend upon the signs of ¢,,, and the rela- 
tive velocity of servo and damper. As long as o,, and 
(Gm — Ga) do not change sign, Eqs. (23) and (24) may 
be integrated to give 


om — Gmo = kis /2 (25) 
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Ca — Gao = kos /2 (26) 
Tn — Om Tmo S + (1/4)Ris? (27) 
where om, 0, and gg are values at time s = 0. 
Combining Eqs. (25) and (27), 
on — Omo = (Cm? — Cmo*)/R: (28) 
Combining Eqs. (26) and (27) 
om — mo = (ki/ke”) toa — [G40 — Gmo(ko/hi) |}? — 
Omo?/ky (29) 
and from Eqs. (25) and (26) 
in — Ome = hi/ke(Ge — Ges) (30) 


From Eq. (28), the trajectory of the servo in the 
phase plane will be identified as a parabola with its 
axis horizontal, and its vertex at 


tt, " Cy Tmo" ky; Fn = 0 (31) 


Similarly, the trajectory of the damper in the phase 
plane is recognized from Eq. (29) to be a parabola with 
axis horizontal and vertex at 


9 


—_ Tmo~ ky 


— Gnoko/k; (32) 


Tm = Tmo 


oa = Oa 


It appears, therefore, that all of the necessary tra- 
jectories are parabolas of the form 


Tm = Fm” ky 


ee ‘ (33) 
Tm = Far (ky k»”) 


which are shifted in the phase plane to satisfy initial 
conditions, but which always keep their axes horizontal. 
The values of k to be used depend upon o,,, upon the 
relative velocity of servo and damper, and upon the 
constants r,and 7;. 

Suppose that servo and damper are initially at rest 
when a constant velocity of input is applied. Assume 
also that the damper has been adjusted to slip. Both 
servo and damper accelerate in a direction to follow the 
changing input, but since the damper slips, its velocity 
lags the servo velocity by an increasing amount. When 
the servo has reduced the difference between output 
and input to zero, the driving contact opens and the 
opposing contact closes. After a delay time fj, the 
driving torque changes sign. The servo then deceler- 
ates, the damper torque aiding the deceleration. 

The servo velocity and the 
velocity increases, until the relative velocity is zero. 
The relative velocity then changes sign, and the servo 
continues to decelerate, but less rapidly, since the 
This 


decreases, damper 


damper is now opposing the change in velocity. 
sequence of events continues to repeat. 

Fig. 2 shows a typical steady-state set of trajectories 
(limit cycles) for servo and damper. The values of 
ki, ko for each trajectory segment are indicated in 
Table 1. 

Because of the symmetry of the phase-plane picture 
given for r, and r;, only four trajectories are needed to 
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TABLE | 
Segment ky ke 
1-2; 1-2’ l r re/Ts 
2-3; 2’-3 —l-r re/ti 
3-4; 3-4’ -l+pr —r./r; 
4-1;4’-]1 l+pr —1/T; 


construct a graphical representation of the motion, 
following any initial conditions. 
phase-plane plot can be constructed, it is necessary to 
establish the contour of torque interchange for con 
stant time lag. 


TIME LAG 


Let subscript a refer to the instant the servo trajectory 
crosses the o,, axis and subscript 2 refer to the point on 
the trajectory sq later, at which torque change occurs 
From Eq. (25) 


F y.2 i Sma oad kiSa 2 (34) 
and from Eq. (27) with o,,. = 0 
Cm2 = GmaSa + (1/4)Rise? (35 


whence the contour of torque interchange is determined 
to be the straight line described by 
4h) RSa? (36) 


Tm. = CmSa — (1 


CONSTRUCTION OF PHASE-PLANE GRAPHS 


Sufficient information is now available to construct 
phase-plane plots of the transient and steady-state 
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Phase plane trajectories for limit cycles of servo and 
damper 
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motion of servo and damper for any desired values 
of the system parameters and initial conditions. 

The transient motion of a clutch-type servo with 
friction damper and the steady-state “‘limit cycles’’ 
eventually attained are shown in Figs. 3, 4, and 5. For 
all three graphs, constant values of dimensionless delay 
time sg and inertia ratio 7; were assumed. The co- 
efficient of friction coupling r, between servo and damper 
was, however, progressively increased. 

It will be observed that the amplitude of steady-state 
oscillation is least for the intermediate value of 7,, 
shown in Fig. 4. It is to be expected that there will 
be a ‘‘best’’ value of 7,, since there is no damping action 
ifr, = Oorif r, = 0.50. In the latter case, the damper 
“sticks,’’ and there is no relative motion and, hence, 
no damping. 

The phase-plane solutions show that there is a limit 
cycle of steady-state oscillation which the servo attains 
and in which it hunts endlessly at a definite amplitude 
and period. It is of interest to obtain a mathematical 
solution for the amplitude and period. 


MATHEMATICAL SOLUTION FOR STEADY STATE 


The procedure will be to establish the equations for 
trajectory segments, determine end relations between 
trajectory segments, including the effect of time lag, 
and finally eliminate from the system of equations those 
phase-plane coordinates that are not desired. The 
result is to yield amplitude and period relations in 
terms of time delay, damper inertia, and friction. 


The notation of Fig. 2 will be used. By symmetry 


m3 = —Omi (. 
O m3 — O81 (38) 
Ga3 = Onm3 = —Omi — Cal (39) 
Between points (1) and (2), by Eq. (28) 
Sa — Gy = Cae — Ca )/ i — %) (40) 
Between points (2) and (3) 
Ooms — Ome = (Cm3” — Cm2*)/(—1 — 7;) (41) 
Combining Eqs. (40) and (41). 
2om2 = (Gm27 — Gm) (1 — 9)/C — 7,2) (42) 
where 
7 = (1 — 7,)/(1 + 7,) (43) 
Writing Eq. (36) in the form 
2(1 — 7,) Ome = AG moSax 2s,,* (44) 
where 
gs. = (1/2Z)sA1 r1) (45) 


Sar 


Then employing Eq. (30) for the damper, and intro 


ducing the end conditions from Eq. (39), 
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coupling 


Om? — Om = (1/S8)(Ga2 — Emi) (46) 
Cm2 + Om = —1/(nS) (642 + Omi) (47) 

and eliminating ¢ from these expressions 
Tmi/ Cm. = —S(1 + y)/(2 + 0S — S) (4S) 
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Variation of maximum velocity in limit cycle with 
friction torque and inertia ratios. 
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Let 


b = —S(1 + )/(2 + nS — S) (49) 


It is now possible to solve for ¢,2. by combining Eqs. 


(42), (44), and (48). The result can be expressed as 


ax? — 4x +2 =0 (50) 

where 
Xx = Ome/Saz (51) 
a= (1 —yn) (1 — 8?) (52) 


From Eq (50) 
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x = 2/a + [(2/a)? — (2/a)]” (53) 


It can be shown that in Eq. (53) only the root obtained 


by taking the positive sign in the sum is of physical 
significance. 


AMPLITUDE AND PERIOD OF OSCILLATION 


The maximum velocity reached by the servo jp 
The 
maximum amplitude is attained at point 0, indicated 
in Fig. 2. Between points ? and 0, from Eq. (28), 


steady-state oscillation is given by Eq. (53). 


Tm2 — Fm0 = Gmo"/(1 — rz) (54) 
Substituting Eq. (44) into Eq. (54) 
—y = (x — 1)* (55) 
where 
X = Om/Saz \ a 
Y = onl — r:)/Saz?h (00) 
To obtain the period P, note that 
P/2 = sio + Seg (57) 
Applying Eq. (25) 
Gm2 — Gm = (1/2) (1 — ris (58) 
m3 — Gme = —(1/2) (1 + 71) 532 (59) 
and combining Eqs. (58), (59), and (48), 
P/sq = 2x[(1 + n) — O(1 — n)] (60) 


or also 


P/sq = 4x(1 — br,)/(1 + 7,) (61) 


CHOICE OF OPTIMUM 1; 


For a given 7, and 7;, therefore, the amplitude of os- 
cillation is proportional to the square of the time lag, 
while the period and maximum velocity are both pro- 
portional to the time lag. In order to show how these 
characteristics depend upon 7, and 7;, Figs. 6 to 8 have 
been prepared, giving the amplitude, maximum ve- 
locity, and period of steady-state oscillation as func- 
tions of r, and 7;. 

Fig. 6 shows clearly that there is a value of r, for any 
given 7; for which the amplitude of steady-state oscilla- 
tion is a minimum. It is interesting to note that this 
optimum value depends only upon 7; and not upon the 
value of the time lag. 

The curves for 7; = © are the same as those that 
would be obtained if the input to the servo were sta- 
tionary and a friction damping torque were obtained 
by locking the flywheel of Fig. 1 to ground. This case 
was treated by Hazen,® who commented that ‘from 
the point of view of making the amplitude of oscillation 
small, a servo with Coulomb friction acting on the out- 
put is preferable to one with viscous friction,’’ but that 
this applied ‘‘only to the condition with a substantially 
stationary input.”’ 
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ratios, and corresponding amplitudes and periods of limit cycles. 


The use of the flywheel coupled to the servo by cou- 
lomb friction preserves the advantage of dry friction 
damping even in the case of moving input. Indeed, 
as shown by the curves of Fig. 6, as long as the optimum 
value of 7, is not exceeded, the friction-coupled flywheel 
provides a damping action for the servo with constant- 
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velocity input which is nearly as efficient as the same 
friction torque applied to a servo with stationary input 
in the manner considered by Hazen. 

Fig. 9 shows the relationship between r, and 7; for 
minimum amplitude of steady-state oscillation and the 
way in which the maximum amplitude and velocity of 
oscillation and the period then depend on r;._ In gen- 
eral, the amplitude of oscillation can be made as small 
as desired by making the inertia of the flywheel large. 
Increasing the inertia of the flywheel while maintaining 
the optimum torque adjustment, however, reduces the 
maximum acceleration of the servo and also increases 
the settling time required for the initial transient to dis- 
appear. 

The curves of Figs. 6, 7, and 8 apply only to the limit- 
ing steady-state oscillation and do not describe the proc- 
ess by which the transient motion converges to this 
steady state. To show the fundamental difference in 
transient response when the damper is free or fixed, 
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phase-plane trajectories similar to those of Figs. 3, 4, 
and 5 were constructed for a servo initially at rest, to 
which a constant velocity input is suddenly applied. 

Fig. 10a shows the variation of maximum error ampli- 
tude with number of cycles of oscillation when the 
damper is free to move in response to the friction coup- 
ling to the servo. Note that the abscissa is number of 
cycles rather than time. 
tually varies with the amplitude, in general decreasing 
as the amplitude decreases. 

Fig. 10b shows the error amplitudes for the same 
servo with the damper fixed so that the flywheel cannot 
rotate. In the latter case, the amplitude increases to 
an extremely large value in the steady state, while, in 


The period of oscillation ac- 


the former, it reduces to a comparatively small value 
as the flywheel is brought up to the same average 
angular velocity as the servo input. In fact, with the 
damper fixed, it is necessary that the absolute velocity 
of the servo output change sign in order to obtain damp- 
ing. This requirement makes the fixed damper unac- 
ceptable for most servemechanism applications. The 
friction-coupled flywheel, on the other hand, is a useful 
device regardless of the nature of servo input 


CONCLUSIONS 


A method has been developed for the analysis of the 
transient and steady-state motion of a clutch-type 
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| 
It 18 
shown that a steady-state oscillation of stable ampli. 


servomechanism with friction flywheel damper. 


tude exists. The amplitude is determined by time ]y¢ | 


in the control, by the ratio of friction torque to drivin 


Ee 


torque, and by the ratio of damper inertia to serve jy 


ertia. Charts are presented showing the amplitu 


and period of the steady-state motion as functions oj 
the parameters describing the servo and damper. 





Summary of Recent Experimental Investigations in the N.O.L. 
Hyperballistics 


(Continued from page 670 


phenomena are noted with condensation at N.O.L. 


and in the C.I.T. tunnel.® 
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The Lift Distribution on Low Aspect Ratio 


Wings 


at Subsonic Speeds 


H. R. LAWRENCE* 


Cornell Aeronautical Laboratory, Inc. 


SUMMARY 


A systematic procedure is developed for obtaining approximate 
integral equations from the linearized theory of a lifting surface 
in a steady incompressible flow. This procedure leads to a sys 
tem of classification of the various theories currently employed 
in the calculation of lift distribution® * * 7 and to a new theory 
for the chordwise lift distribution over wings of low aspect ratio. 
A certain duality between the high aspect ratio and low aspect 
ratio approximations is exhibited by this classification system. 
The following conclusions are established for the special case of 
a rectangular wing, but the results are believed to be generally 
applicable 

(1) There does not exist a low aspect ratio analog to the 
infinite aspect ratio theory. ! 

(2) The low aspect ratio analog to the Prandtl lifting-line 
theory’ is the Jones low aspect ratio theory.” 

(3) The low aspect ratio analog to the Weissinger lifting-line 
theory’ is the present theory. 

(4) The Weissinger lifting-line theory’ and the present theory 
are unique among the theories under discussion in that they both 
approach the exact lifting surface theory when the aspect ratio 
goes to infinity or to zero. 

The theory developed for a rectangular wing is extended to 
more general plan forms. By reinterpreting the integral equa- 
tion, the case of wing-body combinations and more general 
configurations may be solved. However, the régime of validity 
for these extensions of the theory must be established by com- 
parison with test data 

The integral equation for chordwise lift distribution may be 
used in conjunction with the spanwise lift distributions of Jones 
and Spreiter to estimate pressures on the configuration only for 
symmetric spanwise loadings. In the antisymmetric case, an 
integral equation for the chordwise distribution of rolling moment 
may be used in conjunction with the antisymmetric Jones or 
Spreiter spanwise lift distribution to estimate pressures. The 
required integral equation for rolling moments is developed by 
a procedure analogous to the symmetric case. 

A numerical method for solving the integral equations arising 
in the present theory is discussed and applied to calculations of 
the chordwise pressure distribution, lift curve slope, aerodynamic 
center, and damping in roll of rectangular and delta wings for 
the aspect ratio range from 0 to 4. The results are compared 
with test data and show satisfactory agreement. 

The results of the present theory as applied to low aspect 
ratio wings are valid only for small angles of attack. In refer- 
ence 13, a method is proposed for taking account of the nonlinear 
lift curve of low aspect ratio wings. A viscous cross-flow drag 
term 2a? is added to the lift coefficient and is assumed to act at 
the plan-form centroid of area. The present theory may be 
extended to large angles of attack by making use of this pro 


cedure. 


Presented at the Supersonic Aerodynamics Session, Nine 
teenth Annual Meeting, I.A.S., New York, January 29—February 
1, 1951 

* Assistant Head, Aerodynamic Research Department 


INTRODUCTION 


Status of Low Aspect Ratio Theory 


FS posiunna FOR ESTIMATING AERODYNAMIC FORCES is 
necessary in the design of aircraft or missiles 
with low aspect ratio wings or tails. There exist, at 
present, four methods useful in predicting the air 
Bollay" has treated 
By 
making appropriate assumptions as to the location of 
the trailing vortices and the lift distribution, a non- 
linear relationship between lift and angle of attack 
which agrees with test data was obtained. This 
method does not predict lift distributions and is not 
readily generalized to other plan forms. Weinig'! has 
also treated the rectangular wing by a rather obscure 
By making use 


forces on low aspect ratio wings. 
the case of a rectangular wing of low aspect ratio. 


procedure involving cascade theory. 
of the assumption that the momentum of the stream 
intercepted by the streamwise projection of the wing 
is transferred to the wing, a nonlinear relation between 
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the lift and angle of attack was obtained. The results 
agree with test data over the entire range of aspect 
However, the lift distribution cannot be cal- 
culated, and the results do not appear to be readily 
generalized to more complex plan forms. A method 
for predicting the lift distribution on extremely low 
aspect ratio wings with unswept trailing edges has been 
developed by Jones.” For small angles of attack, the 
theory agrees well with test data. 


ratio. 


Since this theory 
predicts a linear relationship between lift coefficient 
and angle of attack, the results diverge from test data 


at higher angles. A method for estimating aero- 
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dynamic forces applicable to all plan forms has been 
Weissinger.* The derivation of the 
Weissinger theory from lifting surface theory has been 
clarified by Reissner.® 


developed by 


The spanwise lift distribution 
predicted by this theory is in good agreement with test 
data for all aspect ratios at low angles of attack and 
for high aspect ratios at unstalled angles. The chord. 
wise lift distribution does not agree with test data for 
low aspect ratio plan forms. Since the Weissinger 
theory leads to a linear relation between lift and angle 
of attack, the agreement with test data on low aspect 
ratio wings at high angles is not to be expected. 


Desired Improvements in Low Aspect Ratio Theory 


An investigation of methods for computing the lift 
distribution on low aspect ratio wings has led to a new 
avenue of approach to this problem. The difficulties 
with present theories may be summarized as follows: 
The theories of Bollay” and Weinig!! predict the non- 
linear effects correctly but do not produce lift distribu- 
tions and are restricted to rectangular wings. The 
theories of Jones’ and Weissinger® predict the spanwise 
lift distribution correctly but do not give accurate 
chordwise lift distribution and omit the nonlinearity 
in lift curve. A theory that predicts the chordwise 
lift distribution, is applicable to. a large class of plan 
forms, and takes account of nonlinearity is required. 


Summary of the Nonlinear Theory of Reference 13 


A summary of the history, motivation, and results 
of reference 13 are presented below as background 
material. The investigation started with a comparison 
of the Bollay’ and Weinig!! theories with test data, 
which indicated that the Weinig theory was more 
accurate. Closer examination of the Weinig!! theory 
established the fact that the cascade analogy and 
momentum terms were independent. The cascade 
analogy is used to establish a potential flow solution 
that gives a lift proportional to angle of attack. The 
momentum term, resulting from the assumption that 
a column of air equal in cross section to the streamwise 
projection of the wing is deflected downward, leads to 
a lift proportional to the square of the angle of attack. 
These considerations suggested the possibility of 
adding the nonlinear term obtained by Weinig to the 
linear potential theory approximations developed by 
Jones* and Weissinger.’ The results of this combina- 
tion of theories reproduced Weinig’s results for rec- 
tangular wings. Application of this method to more 
general plan forms gave excellent correlation with test 
data. Furthermore, the spanwise lift distribution 
could now be estimated for wings of low aspect ratio. 
Although the results obtained thus far were encourag- 
ing, Weinig’s method for obtaining the nonlinear term 
in the lift curve did not appear convincing from a 
physical point of view. An alternate method for ob- 
taining a nonlinear term was suggested by Allen’s 
work on slender bodies of revolution. 
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Allen’? assumes that the lift on a slender body con- 
tains a nonlinear term, resulting from cross-flow viscous 
drag. Investigation indicated that better agreement 
with test data was obtained if, in computing the cross- 
flow drag, the streamwise Reynolds Number was used 
in place of the cross-flow Reynolds Number as proposed 
by Allen. There appear to be sound physical reasons 
for stating that the boundary-layer state is determined 
by the streamwise, rather than the cross-flow, Reynolds 
Number. 

If the cross-flow drag concept is applied to low aspect 
ratio wings, the cross-flow drag coefficient must be 
determined from some hypothesis as to the viscous 
cross flow. A viscous drag coefficient of two, cor- 
responding to a plate of very low aspect ratio, was 
found to give precisely the Weinig nonlinear term. 
Further work of a fundamental nature is required to 
establish the physical validity of either the Weinig 
or the nonlinear term of reference 13 for low aspect 
ratio wings. 

For engineering purposes, the nonlinear effects ob- 
served when testing low aspect ratio wings may be 
predicted by adding a term 2a’ to the potential theory 
solution for lift coefficient as a function of angle of 
attack a and by assuming this lift to be uniformly 
distributed over the plan form. This leads to a great 
simplification in the search for an adequate theory for 
low aspect ratio wings. The Weissinger*® theory pro- 
vides a method for computing the spanwise lift dis- 
tribution, so that only a method for computing chord- 
wise lift distributions is required. 


Linear Theory of Low Aspect Ratio Wings 


A new method for computing the lift distribution on 
low aspect ratio wings is developed in this report. This 
method represents a generalization of the Jones theory 
and is applicable throughout the range of aspect ratio. 
The Weissinger* theory assumes the chordwise lift 
distribution corresponding to high aspect ratio and 
computes the spanwise lift distribution. The present 
theory assumes the spanwise lift distribution corre- 
sponding to low aspect ratio and computes the chord- 
wise lift distribution. The slope of the lift curve com- 
puted by these two widely divergent theories leads to 
almost identical results over the entire range of aspect 
ratio 0 to ». The choice of theory would then depend 
on whether accurate values for the spanwise or chord- 
wise lift distribution are required. If a complete lift 
distribution is required, the two theories can be com- 


bined in an appropriate manner. 


List OF SYMBOLS 


U = velocity of free stream 

u(x, y) = parametric lift per unit area 

x,y; = Cartesian coordinates of point of observation 
§ 7 = Cartesian coordinates of source point 

P(x, y) = 2oU%u(x, y) = lift per unit area 


= free-stream density 
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VARIATION OF LIFT SLOPE WITH ASPECT RATIO 
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¢(x, ¥, 2) = perturbation potential K'(x) = 2pU*k'(x) = Jones lift per unit chord 
—w(x, y) = local slope of airfoil in stream direction g'(x) = see Eq. (30) 
b = wing semispan G'(x) = 2pU%g'(x) = lift per unit chord 
c = wing semichord H(@B7r) = see Eq. (50) 
2 A(x) = local wing span H,,(B) = see Eq. (60) 
hy) = see Eq. ( 10) ; ees s = wing area 
H = 2oU*h = lift per unit span for wing of infinite AR. = aspect ratio 
aspect ratio i 4 ren 
/ e'(x) = see Eq. (45) 
Ky) = ¢ déu(é, vy) = parametric lift per unit span E'(x) = 2pU%e'(x) = Jones rolling moment per unit chord 
L(y) = 2pU*l(y) = lift per unit span m'(x) = see Eq. (46) 
k(x) = see Eq. (21) M(x) = 2pU?m'(x) = rolling moment per unit chord 


(1) INTEGRAL EQUATION FOR A LIFTING SURFACE 
Eq. (62) of reference | will be used as a starting point for this investigation. With the notation 
fuk - | + & - 9 + # 


Eq. (62) may be written 


| Fi jx —€& O| (x — é)(y — n) y-n ) 
(%,y,3) = — IS 1é dn u(é, n) - | ‘ a | 1) 
G(X, | = dé dyn ul&, n i» + oy r[s? + (y — n)?] a 22 + (y — 7)? ( ( 


The airfoil is assumed to lie in the plane z = 0 with local angle of attack w(x, y). ¢,(x, y, 2) represents the z 
component of the perturbation velocity at the point (x, y, 2), and u(é, 9) is proportional to the jump in the x 
component of the perturbation velocity on crossing the plane z = 0. The lift per unit area p(x, y) is given in 
terms of u(x, y) by the expression 

p(x, y) = 2plu(x, y) (2 


Eq. (1) may be rewritten in the following forms: 


l ¥ Oo f(x — E)(y — 7) (x — §)(y — n) y—-n 
OAS; 9; 2) = = IS dé dyn u(é, n) ‘ - a a = es i aa A | 
Qa J a Oy (r[z? + (x — £)?] r[2z? + (y — n)?] 2? + (y — »)*f 
(3) 
1 Oo a % yr (x — &)(r? + 2?) 
(x, y, 2) = — / / dé dn u(&, n) | a : 2 J + zs - - l 
27 Oy J« 2? + (y — n)’ { r[22? + (x — £)2)f 
For the case of steady flow, $,(x, y, 0+) = —w(x, y) on the lifting surface and u = 0 over the remaining por- 
tion of the plane z = 0. Upon substituting these boundary conditions in Eq. (3), there results the integral equa- 
tion 
1 oO -_ u(é, 7) V (x — €)? + (y — 7)? 
w(x, y) = | / ie 41 4 : — (4) 
2r Oy J Ra yn x—f& 


relating the unknown velocity u(x, y) to the known angle of attack of the lifting surface w(x, y). «is proportional 
to the lift per unit area and will be called the parametric lift per unit area. 


(II) INFINITE ASPECT RATIO 


The integral Eq. (4) may be solved for the case of infinite aspect ratio when the angle of attack is independent 
of the spanwise coordinate. Eq. (4) may then be written as follows: 


l " és a) | V (x — &=)2? + (vy — n)? 
w(x) = — / ae f dn u(&) + ae = 
i, = Only —n (x — E)(y — n) 
and, upon performing the integration over 7, 
l © dé u(é) . 
w(x) = 2 (5) 
T cx—€é 


Eq. (5) is the well-known infinite aspect ratio result. The solution of this equation for the parametric lift per 


unit span, h, is given in reference 4 as 
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k= / dé u(é) = 


LIF 


The lift per unit span, IT, is related to h by the expression 


H = 


(IIT) 


The case of a rectangular lifting surface is considered so as to simplify the analysis. 


integrated by parts and written in the form 


| *b a 
w(x, y) = - dn dé u,(é, 7) 
a = é 


The integral Eq. (S) will be satisfied on the average in the chordwise direction. 
Eq. (8) by a weighting factor and integrating over the chord. 
Multiplying Eq. (S) by Vc + x/Vc — x and integrating over the chord, there results the 


c+ x 
Ix 
c—x 


wing of infinite span. 


equation 


h(y) = 


where /i(y) is defined as follows: 


/ dx w(x, y) 


It will be noted that the definition of h(y) is consistent 
with the interpretation of h(y) as the parametric lift 
per unit span of a wing of infinite aspect ratio |[cf. 
Eq. (6) |. 

Now an approximation based on the assumption of 


CT x 


c-— F¢ 


h(y) = (10) 


large aspect ratio is employed to simplify the integral 
Eq. (9). 


J'= [ dx V + * [ 4 ¥ Ge - 2 — ”) 
. =.= x—& 


Consider the integral 





(1 | ) 

Since the aspect ratio is large, y — »| will, in general, 

be much larger than |x — &). Hence, as a first ap- 
5 s 

proximation, the radical V (x — &)* + (y — n)* may be 


replaced by |y — 7}. Then the integration over x 


may be performed, and Eq. (11) may be written 
Ji’ = a#(c + |y — 9}) 


Upon substituting this result in Eq. (9) and making 
use of the definition for parametric lift per unit span 


I(y), 
I(y) = / dé u(&, y) 


the integral Eq. (9) acquires the form 


(12) 


2pl*h 


y-7 


, > be a 5 
> dé u,(&, 7) ti 
2 sy—7 J - 
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“¢ c+ex 
dx w(x) (6) 
r = =F 
(7) 


PRANDTL LIFTING-LINE THEORY 


Then Eq. (4) may be 


+ 


l . Vie— sy + yy —-»)* 
oe 


(x — &) 


This is achieved by multiplying 


The weighting factor is chosen to be exact for a 
V (x 2. 
i +r , (9) 


h(y) = 


. *b 1 1 
= I(y) + = / —— (13) 
2 J-»y —ndy 


The lift per unit span L(y) is defined in terms of the 
parametric lift per unit span by the relation 


Liy) = 2pUl(y) (14) 


Eq. (13) is the Prandtl lifting-line integral equation 
for a rectangular wing. 

(IV) Jones Low Aspect RATIO THEORY 
The integral Eq. (4) may be solved for the Jones 
only the 
Then Eq. 


low ratio case. For convenience, 


rectangular plan form will be considered. 


aspect 
(4) may be replaced by Eq. (8) or 


l “> dn : . 
dé u,(&, nn) X 
2r . ‘I-% 


V (x — &)° + (v — n)? 
1+ = = (S) 
x—é€& 


w(x, y) = 


wtr 


a 
For sufficiently low aspect ratio, the radical 
V (x — &)? + (y — »)? may be replaced by |x — & 
Then Eq. (8) may be simplified to the form 
I b dy - 
Wx, Vv) = od, (xX, 7) (15) 
T T-F 
where (x, y) is defined as follows: 
"6 
d(x, y) = / dé u(é, v (16) 


. 
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The integral Eq. (15) is equivalent to the low aspect 
ratio theory of references 2 and 5. A solution of Eq. 
(15) is given by the formula® 


—1 > dn w(x, 7) Vb? — n° 
$,(x, y) = j (17) 
' rV b? — wt 5 ie 


In the case where w is independent of y, Eq. (17) may 
be integrated to give the familiar result of Jones’s 
low aspect ratio theory for a flat plate, 


(x, y) = w(x) Vb? — y? (18) 


The lift per unit chord K’(x) and the parametric lift 
per unit chord k’(x) are defined as follows: 


ik : 
k’ (x) a e = f dy u(x, y) 
dx —b 


dK 
K'(x) = - = 2oU?R’(x) 
dx 


(19) 


It is assumed that k is zero at the leading edge of the 
plan form. The relation between ¢ and k is given by 
the formula 


b b 
k(x) = / dy o(x, y) = — f dy y,(x, y) (20) 
- —b 


where ¢(x, Upon multiply- 
ing Eq. (17) by (—y) and integrating over the span, 
the parametric lift k(x) may be written in the form 


b : 
R(x) = y dy w(x, y) Vb? —¥? (21) 
-b 


It is important to note that the approximation made 
in the Prandtl lifting-line theory consists in replacing 
V (x — §)? + (y — n)? by ly — n| and that the ap- 
proximation of the Jones low aspect ratio theory 


+ (y — n)? by |x — é}. 


+b) is assumed to vanish. 








replaces V (x — §)? Hence, 
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these cases are closely analogous. It is rather curioys 
that no approximation exists in the low aspect ratio 
case that corresponds to the case of infinite aspect 
ratio. The first approximation to a low aspect ratio 
wing leads to the Jones theory, which corresponds to 
the second approximation in the case of high aspect 
ratio (Prandtl lifting-line theory). The shortcomings 
of the Prandtl theory and a method of improving the 
theory will be taken up next. 


(V) WEISSINGER LIFTING-LINE THEORY 


There is no reason to believe a priori that the high 
aspect ratio lifting-line theory can be applied at low 
aspect ratio. However, a lifting-line theory that is 
correct in the limit as the aspect ratio approaches zero 
would be expected to give fair results throughout the 
aspect ratio range. Such is the case for the Weissinger 
theory. A procedure for developing the Prandtl and 
Weissinger theories similar to that employed in this 
paper was presented earlier by Reissner.* Consider 
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first the limiting case approached by the Prandtl 
lifting-line integral Eq. (13) as the aspect ratio ap- 
Making use of the definition (10) for 


proaches 0. 
13) may be written as 


h(y), Eq. ( 


in =) ; 1 + (x/c) 
d w(x, y) = 
a c ; 1 — (x/c) 
Sf, dl 5 + ( c) 
dn 
2 b dn 
Upon transforming the integral on the left by the 


substitution A = Xx/c noting that |y — 7n|/c 
becomes negligible compared to unity as c approaches 


yn 
dah 


and 


infinity, the limiting form as c ~ © of the above 
integral equation may be written as 
“I 1+A l “6 dyn dl 
d(A)w(aAc, y) = / - 
, L—A 23-69 — Has 


If w is independent of x, the above equation becomes 


*b 
dn dl 

2rw(y) = / (22) 
J -py —ndn 


Eq. (22) represents the limiting form of the Prandtl 
lifting-line integral equation as the aspect ratio becomes 
small. 

The integral Eq. 
aspect ratio theory is rigorously correct in the limiting 


(15) representing the Jones low 
case of low aspect ratio. 

Making use of the relation /(y) = 
Eq. (15) may be written 


*b dn dl 
Tw(y) = 
—_by —ndn 


where it is assumed that w is independent of x. A 
comparison of the correct result for low aspect ratio 
given by Eq. (23) with the low aspect ratio limiting 
form of the Prandtl lifting-line theory given in Eq. (22) 
shows that the Prandtl integral equation is in error by 


o(c, y), the integral 


(23) 


a factor of two. 
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The Weissinger improvement on the Prandtl theory 
is found by consideration of Eq. (11). Instead of ap- 
proximating to the radical V (x — §)? + (y — »)? by 
|_y — n|, a more accurate approximation ¥V (y —»)?* +c 
is used. Evaluating the integral in Eq. (11) with the 
above approximation leads to the result 


J,’ = rc +V (y—n)?t+ c? | 
which may be substituted into Eq. (9) to give the 
Weissinger integral equation 
l ~~ dl \}crt Ve + (y — n)? 
h(y) = / dn ‘ (24) 
2 J+ dy yr” 
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As in the Prandtl case, the limiting form of Eq. (24) 


or 
1 ’b dl is easily obtained as the aspect ratio becomes small, and 
h(y) = Uy) dy — X ar , ; le 
) + S$ Js ’ dn the result is identical with Eq. (23). Thus, the 
c+ Ve +(y—n)?-ly—- n| 54:1) Weissinger theory leads to the correct limiting case 
aa ~ as the aspect ratio goes to zero. 


(VI) A Low Aspect RATIO APPROXIMATION 


In view of the background material presented, a method of improving the Jones theory is obvious. The weight- 
ing factor V b? — y* obtained from Eq. (21) will be used to satisfy Eq. (4) on the average in the spanwise direction, 
The Eq. (21) ensures that the integral equation will be satisfied exactly in the spanwise direction as the aspect 
ratio approaches 0. Multiplying Eq. (4) by Vo? — y® and integrating over the span, the following result is 


obtained for a rectangular wing: 


| “b y, ra) i “6 u(é, 7) V(x — £)2 (y — n)? 
R(x) = / dy V b? — y? / dé / dn =i + ye . (25 
Le Ff ab Oy , ya y—n x—€& 


where k(x) is defined in Eq. (21). After integrating by parts with respect to y and interchanging the order of 


integration, Eq. (25) becomes 


| . b s "6 dy y l V (x — &)? vy — n)” 
k(x) = - [ af dn u(é, n) / — ( ) 1+ - =e ? (26 
aa r b bY bh? — y? ¥y —s ~~ 2 


‘ <4 > » lact 7; > y . TI rR IG -\9 9 +\ ¢ 1 9 
Consider the last integral on the right of Eq. (26), V (x — £)2 + (y — 9)? = Vx — 2? + 
iis dy y as an improvement over the ; imati 
on 2M x as ¢ | inent over the approximation 
J-b (y — 9) V5? — y’ 
Vix —§&)*+y—n’=Iix-—E 
V (x — &)? + (y — 2) 7 “ ' 
1+ (27) _ ‘ 
x—t hen the integral J may be evaluated as 
For sufficiently low aspect ratio, the radical Je = r} 1+ [V(x — £)? + B7/(x — r)]} (29) 
V(x — &? + (y — 7)? may be approximated by 
x — £|. Then J may be evaluated as Upon introducing the parametric lift per unit chord 
2 ; I 


. g’(x) and the lift per unit chord G’(x) defined as follows: 
Ji = nfl t+ [|x —eé]/(e—8]} (27.1) 


~~. 
and, upon substituting this result in Eq. (26), there g"(x) -f dy u(x, »| - 
) (e } 

results . 

G'(%) = 2pUe"(x) 
*b 

Rix) = In b(x, 7) 28) ite i . , ; 

b ai ata \ and substituting the result, Eq. (29), into Eq. (26), 


the following integral equation is obtained: 
Comparing Eq. (28) with Eq. (20), it is evident that 


this approximation leads again to the Jones result. 1 Vix — 2+ B? 
This analysis bears out the assertion that the Jones = &(x) = 5 / dé g'(é) | 1 + io (31) 
theory for low aspect ratio wings is the analog of the 7 . 
Prandtl lifting-line theory for wings of high aspect 1 c 
ratio. R(x) = g(x) + 3 [ dé g’(é) X 
It is natural at this point to follow along on the of ae 
analogy between low and high aspect ratio theories ee -otr~i~ 8 (32) 


by trying the approximation x 
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Here, k(x) is defined by Eq. (21). The resemblance 
between Eqs. (31) and (32) and the Weissinger integral 


Eqs. (24) and (24.1) is evident. 


It is of interest to compute the limiting form of Eq. 
(31) as the aspect ratio becomes large. For simplicity, 
the case where w is independent of x is considered. 
Then, upon integrating Eq. (21), the integral Eq. (31) 


may be written 


rbw | ie : \ Victi{[(x -—é pz! 
= lg 7 (€) , . 
a? fare jy + Ste Oar 


and for sufficiently large values of 6, Eq. (31) may be 


written as 


(32.1) 


(VII) ImMpROvED Low 


forms: 
™ dy y Ix — &| “ 
J= i l + + 
P b (yv —_ n)v b? —_ yy = é —b 
ix —€ dy y 
J=nr{1+ - 
e~§ 


Now in Eq. (35), the approximation 


/ 


V(x —éP + y—a’?~Ve- 


is introduced. It will be noted that approximation to the radical by V (x — ¢ 
For example, when 0 < y — n < 3, the following inequali- 


ferent result than the same approximation in Eq. (34). 


ties apply: 


LOW 


*b 
b (x — t)Vb? — y2 LV @& — &)* + (y — 9)? + x - 
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For a rectangular wing of sufficiently large aspect 
ratio, the lift per unit span is substantially constant. 
Hence, from Eq. (30), g’(x) = 2bu(x), and substituting 
this result in the integral Eq. (32.1) leads to the follow- 


l © dé u(é) a, 
W(x) = o¢ (3) 
2r c x—€ 


Upon comparing Eqs. (33) and (5), it is seen that the 
integral Eq. (5) differs from the limiting form of Eq. 
(31) as the aspect ratio becomes large by a factor of 
two. It is believed that one great factor in the success 
of the Weissinger theory is the fact that it is correct 
ratio. Hence, 


ing equation: 


in the limiting case of low aspect 
although Eq. (31) is probably an improvement on the 
Jones theory, it cannot be considered entirely satis- 
factory. An integral equation that eliminates this 
objection is developed in the next section. 


ASPECT RATIO THEORY 


In order to obtain a better approximation to the integral J of Eq. (27), the integral is rewritten in the following 


V(x — £)? ++ (vy — 9)? — |x -—€& ™ 
(x — EY — ®) ; 


dyy 


Vb? — sia 


S aa 
(35) 


rer 


+ 


£)? + 6’ in Eq. (35) leads to a dif- 


V(ix—&)?+07?-—|x-—t Vw-—t)? + (y— nn)? —-— |x -E 
> i 
y—-” y-é 
y~ 9 y~-9 
kt > 
Vix —'?+(y—7)*?+ |x -—€E Vix —é?+h+ ix fs 
With the above approximation, the integral expression for J may be evaluated as follows: 
x —€& 1b l I 
] w{il-+ t ae 3 % ¥ (36) 
x—€& 2 | V(x — £)? + BF + |x = &] | @ — &) 
V (x £)? + & ies 
(54) 


Upon comparing Eq. (37) with Eqs. (2 
of the Jones and the first improved approximation 


J; = 


(1/2)(i 4 


x—& 


27.1) and (29), it is seen that the approximation /; represents the average 


J.) 3S) 
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The expression, Eq. (38), for J; is substituted in the integral Eq. (25) to give the second improved low aspect 


ratio integral equation 


1 | ™¢ ; \ V (x 
R(x) = 5 g(x) + / dg g'() 1 + —— 


2 f 


\ (39) 


By comparing Eq. (39) with Eq. (32), it is evident that Eq. (39) reduces to the following form as the aspect 


ratio becomes large: 


l 


7 


2bw = 


when w is independent of x. 


fact, represent a correct limiting case as the aspect ratio becomes large. 


*e dég’(£) ‘ 
: (40 
cx—€ 


> 


From the comparison shown in Eqs. (32.1) to (33), it is clear that Eq. (40) does, in 


Eq. (39) will be taken as the improved 


Jones theory integral equation; Eq. (32) will no longer be considered. 


(VIII) 


MorRE GENERAL PLAN FORMS 


The integral Eq. (39) proposed as an improvement over the Jones theory for low aspect ratio wings has been 
@ o 5S 


developed thus far for rectangular plan forms. 


In the following section, this result will be generalized to other 


cases by means of further approximations where necessary. 


Consider a plan form with a trailing edge normal to the undisturbed stream. 
nondecreasing in the direction of the undisturbed stream. 


approximation as follows: 


This implies that the span is 
It is evident that Eq. (26) may be written without 


I ‘. “B(E) “BE dy y V(x — 2+ (vy — 7)? 
k(x) = J dé / dn u(é, 7) / = 4 - : (42 
Le pm — B(E) J-Be) V B(x) — yLyY — 0 (x — &)(y — ») 
The approximation used in going from Eq. (35) to Eq. (39) consisted of replacing V (x — &)* + (y — n)* by 
V(x — &)? + Bb in Eq. (35). It may be shown that the replacement of Vix — &)% + (y — 7»)? by 
V (x — £)? + B(x) in Eq. (42), after transforming the radical to the denominator, leads to the proper generaliza- 


tion. Then, upon performing the indicated integrations, Eq. (42) is replaced by the approximate integral equa- 


tion 


i 
R(x) = < g(x) + 4 dé g’( 


This generalization is equivalent to the assumption 
that each spanwise section acts like an element of 
rectangular wing with the local span 26(x) and chord 
20. 

The definition of k(x) is altered slightly to the form 


B(x) j 
R(x) = 7 dy w(x, y)V B(x) — yx? (44) 
— B(x) 


The application of this theory to more general plan 
forms may be carried out along the lines given above. 
However, for wings that are influenced by their wakes, 
even the Jones case becomes complicated. 

A more general procedure can be evolved from an 
interpretation of the results obtained thus far. Eq. 
(43) gives a relation between the lift computed from 
the Jones theory k(x) and the three-dimensional lift 
g(x). k(x) may be considered as two-dimensional lift, 
since it depends only on the geometry at the chordwise 


station x. If, now, k’(x) is redefined as the two- 


V (x — &)? + B(x) 


f)} 1+ - (43 


x —é 


dimensional parametric lift of the spanwise element, 
the integral Eq. (43) may still be expected to hold 
approximately. Now, however, the configuration may 
be a wing-body combination with the two-dimensional 
lift per unit span given by Spreiter® or a more general 
case. This idea is entirely analogous with the treat- 
ment of a wing nacelle as a region of low lift curve 
The method 
appears to be more useful in the present context, since 


slope in the Prandtl lifting-line theory. 


the two-dimensional lift may be computed by the 
Spreiter theory, whereas in the high aspect ratio case 
the two-dimensional lift was difficult to estimate. 
Of course, it will be necessary to determine how far 
this interpretation can be carried. 


(IX) IMPROVED Low AspEcT RATIO THEORY FOR 
ANTISYMMETRIC CASE 


The pressure distribution over a wing with a span- 
wise symmetric local angle of attack distribution may 
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be estimated by combining the parametric chordwise 
lift distribution g’(x) determined from Eq. (39) with 
a spanwise lift distribution estimated by the methods 
of Jones* and Spreiter.6 When the spanwise distribu- 
tion of local angle of attack is antisymmetric, the 
chordwise lift distribution vanishes and the above 
method for estimating pressures is no longer applicable. 
However, when the chordwise distribution of rolling 


Eq. (4) is multiplied by —(1/2)y+1/ 6* — y* and integrated over the span. 
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moment is known, the pressure distribution over the 
wing can be estimated for an antisymmetric local 
angle of attack distribution. 

The procedure for obtaining an integral equation 
for the chordwise distribution of rolling moment 
is so closely parallel to the symmetric case that 
the development will be pre- 


only an outline of 


sented. 


The resulting integral equation when 


simplified by the approximation employed in Section (VII) leads to the following formula: 


| *b | ( 
aa | dy yV Bb? — y? w(x, y) = / dé 
2 J a 


By making use of the definitions 


Dt / 
e(x) = — | dy yV b? — y? w(x, y) 


—b 


; dm 7 
m'(x) = = dy yu(x, y) 
dx ~b 


1+ 


V(x — £)2 + B? 
+14 


x—é& 


*b 
/ dn nu(é, n) 
J 


and performing an integration on & the above equation may be written as 


l l 
e(x) = rs m(x) + 


2 mg 


(45) 
(46) 
“¢ , v fe m= ¢)? + ? 
/ dé m'(&)| 1 + (47) 
P x—€& 


It will be noted that Eq. (47) becomes identical to Eq. (39) when k(x) and g(x) are substituted for e(x) and m(x). 


From the definition (46), m’(x) may be called the parametric rolling moment per unit chord. 
It may be shown that e’(x) defined in Eq. (45) represents 


2pu2m'(x) becomes the rolling moment per unit chord. 


the parametric rolling moment per unit chord given by the Jones low aspect ratio theory. 


Then M/'(x) = 


The corresponding 


rolling moment per unit chord £’(x) is given in terms of e’(x) by the formula E’(x) = 2pU%e’(x). 
As in the symmetric case, the generalization to plan forms with a straight trailing edge is found by replacing 


b by B(x) in Eqs. (45), (46), and (47), where 28(x) is the local span. 


Since e(x) represents the Jones low aspect 


ratio solution,” * © the integral Eq. (47) may be used for more general configurations, such as a wing-body case. 
As in the symmetric case, the range of validity of these approximations must be established by comparison with 


test data. 


(X) 


A method for solving the integral Eq. (43) is presented below. 


SOLUTION OF THE INTEGRAL 


SQUATION 


It will be observed that the same procedure 


is applicable to the integral Eq. (47) and that a large portion of the numerical calculations is common to both 


equations. 
' dé g’(£) 


4k(x) = 2g(x) + g(1) + B(x) Z 
a 2. = é 


where the semichord is taken as unity. 
The transformation x = cos 0, — = 
function: 


1k(0) = 2¢(0) + g(0) + B(A) / 


70 


where 


H(0, B, r) = [ V (cos r — cos 6)? + B°(6) — 8(6) | (cos r — cos 6) 


Vv (: 
+ dé g’(é) 
-1 x—£ 


= cos 7 reduces Eq. (48) to the following form, where g is considered a point 


7 
T 


cos tT — cos 6 


The integral Eq. (43) may be written in the form 


x — §)? + B(x) — B(x) 
(48) 


Ir g'(r) 
aT g v (49) 


+ [ dr g(r) H(0, r) 


(50) 
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A series expansion of g(#) which satisfies the Kutta-Joukowsky condition and vanishes at the leading edge js 
given below, where, by definition, Ay_1 = Ay = 0. 
\ | 
g(0) = (4 — 0)(Ao + Ai) + +» (A,_1 — A,s1)(sin 76 /r) (51) 


? l 


Upon substituting the expansion (51) for g(@) in Eq. (49), the resulting equation may be reduced to the form 
N —2 

(4/r)k(0) = [F\(0) — Fo(@)JAo + SO [Fyui(0) — F,-1(0)]A, 52 
, l 


where the functions F,.(@) are defined as follows: 


Fy(0) = +(20/r) + (6, B) — 3 
; 2 sin 70 B(@) sin ré (53) 
F.(0) = + . + H,(0, 8) r=1,2,.. | 
Tr sin @ 
l °m a 
H,(6, 8) = dr cos rrH(6, B, 7) (54 
Tr JO 
With the notation k, = k(nr/N) and F, , = F,(nr/N), the following set of simultaneous equations is ob- 


tained from Eq. (52) by collocation: 


N-2 
(4/x)Rn = (Fin — Fon)Aot+ D> (Frain — Frit n)An n= 1,2,...,N—1 (55) 
’ l 
For the wing alone case, k = (7/2)8°w. Upon substituting this result in Eq. (55) with w equal to 1/(2b°) 


where b is the semispan, there result the equations 
N—2 
(Fi, —_ Fo,)Ao + _ (Frat. = F. "a alts = (Br b)?, n= l, 2. eo N — | (56) 
r 1 ° 
If the A, satisfy the Eqs. (56), the lift curve slope C,, is given by the expression 
Ci, = Ir A.R.(Ao + A}) oe (57) 


and the center of pressure is located aft of the leading edge a fraction dCy/dC, of the root chord, where 


dCy ‘dC, = (1 2) —_ (1 +) [(Ao inte A») (Ao + A})] (5S) 
The parametric lift per unit chord dg/dx is given by the formula 
N-2 

dg/dx = A, tan (0/2) +2 >> A, sin r6 (59) 
r 1 


where the angle of attack w is taken as 1/(20°). 

Numerical evaluation of the parameters /7,,,(8,,) accounts for a large portion of the total time required to obtain 
a solution of Eqs. (56). The //,,, are seen to be the coefficients of the Fourier series development of H7(@, 8, 7) 
in terms of the variable +. As shown in reference 9, the trapezoid rule is the preferred procedure for numerical 
integration under these circumstances or 


l Ls 
fi BL) = [ dr cos rrH(6,, Bn» T) = 
tT Jo 
VU l 


1 | A7(6,, B,, 0) + (—)'H(6,, B,, 7) mr mr 
+ cos Hté,, &.. (60) 
M y M M 


m l 


The //,,(8) may be computed as functions of 6 and plotted. The results may then be used in all further cal 
culations. 


(XI) NUMERICAL RESULTS eters /7,,,(8)} were computed and are plotted as functions 

of 8 in Fig. 1. The slope of the lift curve, center of 

Calculations were performed using three points for pressure, damping in roll, and representative chordwise 
rectangular- and delta-wing plan forms. The param-_ lift distributions are plotted and compared with test 
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data and with other theories where applicable. In 
applying these results to the curves of lift and moment 
versus angle of attack, the theory of reference 13 is 
used. 

It will be noted that the chordwise lift distribution for 
low aspect ratio rectangular wings is negative over a 
portion of the chord, indicating that the three-point 
collocation procedure is not satisfactory. The results 
shown in Fig. 4, obtained from a five-point calculation, 
differ somewhat from the results from the three-point 


collocation procedure. 


REFERENCES 


Reissner, Eric, On the General Theory of Thin Airfoils for 
Nonuniform Motion, N.A.C.A. T.N. No. 946, August, 1944. 

2 Jones, R. T., Properties of Low-Aspect-Ratio Pointed Wings 
at Speeds Below and Above the Speed of Sound, N.A.C.A. Report 
No. 835, 1946 

Weissinger, J., The Lift Distribution of Swept-Back Wings, 
N.A.C.A. T.M. No. 1120, March, 1947. 

Munk, Max M., General Theory of Thin Wing Sections, 
N.A.C.A. Report No. 142, 1922. 

> Lomax, Harvard, and Heaslet, Max, Linearized Lifting 
Surface Theory for Swept-Back Wings with Slender Planforms, 
N.A.C.A. T.N. No. 1992, December, 1949. 

6 Spreiter, John R., Aerodynamic Properties of Slender Wing 
Body Combinations at Subsonic, Transonic, and Supersonic 
Speeds, N.A.C.A. T.N. No. 1662, 1948. 


LOW ASPECT RATIO WINGS 695 


7 Prandtl, L., Tragfligeltheorie. I. Mitteilung, Nachrichten 
der K. Gesselschaft der Wissenschaften zu Gottingen, Math.-Phys. 
Kiasse, pp. 451-477, 1918. 

® Reissner, Eric, Note on the Theory of Lifting Surfaces, 
Proceedings of the National Academy of Sciences, Vol. 35, No. 4, 
pp. 208-215, April, 1949. 

® Milne, William Edmund, Numerical Calculus, p. 302; 
Princeton University Press, 1949. 

1° Bollay, W., A Nonlinear Wing Theory and Its A pplication to 
Rectangular Wings of Small Aspect Ratio, Jb. d. Luftf.-Forschg., 
P.I. 139, 1939. 

" Weinig, F., Liftand Drag of Wings with Small Span, N.A.C.A. 
T.M. No. 1151, 1947. 

2 Allen, H. Julian, Pressure Distribution and Some Effects of 
Viscosity on Slender-Inclined Bodies of Revolution, N.A.C.A. 
T.N. No. 2044, March, 1950. 

18 Flax, A. H.,.and Lawrence, H. R., The Aerodynamics of Low- 
Aspect-Ratio Wings and Wing-Body Combinations, to be pre- 
sented at the Third International Joint Conference of the 
R.Ae.S.-1.A.S., Brighton, England, September 3-14, 1951 

14 Tosti, L. P., Low Speed Static Stability and Damping-in- Roll 
Characteristics of Some Swept and Unswept Low-Aspect-Ratio 
Wings, N.A.C.A. T.N. No. 1468, October, 1947 

18 Wind-Tunnel Tests on Cornell Aeronautical Laboratory 
Static Stability Missile Models, GALCIT Report 539, May, 1949. 

16 Ribner, H. S., The Stability Derivatives of Low-Aspect-Ratio 
Triangular Wings at Subsonic and Supersonic Speed, N.A.C.A 
T.N. No. 1428, September, 1947. 











no less than !/, in. high. 





Notice to Contributors 


Manuscripts submitted for publication must be double- or triple-spaced originals. There should be wide margins on 
both sides of the sheets and triple spacing around formulas to allow for the marking of diréctions for the printer. 


Original drawings (jet black India ink on white paper or tracing cloth) should accompany the manuscript. (Blueprints 
are not acceptable.) Photographs must be on glossy white paper. 


Only the simplest formulas should be typewritten; all others should be carefully written in pen and ink. 
between capital and lower-case letters and Greek symbols should be clearly distinguished. 


See Inside Back Cover of this issue for further details. 


The smallest lettering on 8- by 10-in. figures should be 


The difference 


























J of such integration is substituted into the momentum integra] M 
e Cum relation de 
sul 
Oo 6 F F a! 
— a) dy 
Ty, = dx Jo pu (My u) ay 3 : i 
RIEF REPORTS of investigations in the aeronautical sciences “ete the x-direction is parallel to the surface of the flat plate, j 

c oi So stig ; . 
and discussions of papers published in the JOURNAL are represents the total thickness of the boundary layer, and the = 
presented iu this special department. Publication will be completed subscript ; tears the pine ato to pee main free lat 
: : ‘ a ath am. T re i ‘ >W é e turbule ic as 
approximately 8 to 10 weeks after receipt of the material. The Edi- — . ; mi _ hen er wags . essa Sar She Sermetent ih eq 

torial Committee does not hold itself responsible for the opinions ex- Ee AE Ree Saee eeeee ae oN 
pressed by the correspondents. Contributions should not exceed F(A) 2\'" ly, CeR , (: 42> I wi 
800 w 7 aad = ; =-—)\In Cr,X; — win M,? j+ : 
800 words in length. ly — )M,2/2]'"“\ Cr, k\ A\ 3 wi 
1 inl (2a — 1) (y — 1)M,? }] co 
n + ae i 0 
+ 2 211 + (y — 1)M,2/2) Jf P 
su 
where in, 
ba 


Effects of Density Fluctuations on the 
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Cr, = f, 7, ax / = piu,;2? x = mean skin-friction coefficient 
| as A FULLY DEVELOPED turbulent flow problem, it is im- = 
portant to estimate the Reynolds stresses due to the eddy- i wine te Sp 
ing motions of the fluid masses. In the compressible flow regime, 17 
these Reynolds stresses consist of three major parts—viz., (1) w = log (u,/m)/log (T,,/T)) 
that due to velocity fluctuations, (2) that due to density fluctua- caf ee } y uo Jo 
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a reasonable estimate of the effects of density fluctuations at 
high supersonic speeds 























Extending von Karman’s similarity considerations and 
Prandtl’s mixing length concepts to the compressible flow regime, 
we are able to derive the following fundamental differential equa- 
tion representing the turbulent friction on a flat plate: 
du dp\ , , di 
= (> dy sli a idl dy (1) 
O 
In the derivation of Eq. (1), we have assumed that the flat plate 
is lying parallel to the main stream. @ and p designate the mean 5 
velocity and the mean density of the flow within the turbulent A 
boundary layer; they are essentially functions of y, which is the I 
distance in the direction perpendiculir to the flat plate. 1, is 
the shear stress at the surface of the flat plate. & is the incom- 
pressible mixing length constant. @ is a proportionality con- I 
stant that is defined as the compressibility mixing length con- o1 
stant. The term containing a@ represents the contribution from tl 
density fluctuations. | | Ee 
Crocco® has shown that the differential equations for compres- | | | | | as 
sible flow, including viscosity and heat conductivity, have a Retin” 
particular integral of the following form: : | ] | ] 1 | | | in 
| . 
(u?/2) + CpT = constant (2) | | | | | | | | hh 
provided that no heat is added from external sources and that 
the Prandtl Number is unity. We shall assume that Eq. (2) aie . - : . > 7 . . are 
holds for the mean flow in the turbulent boundary layer on an , M 
insulated flat plate. We can then integrate Eq. (1). The result Fic. 1. 
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M, < 2.5—we may take a = 0. On the other hand, when the 
density variation is pronounced -for instance, in the higher 
supersonic range—we have to use some other values of a. Thus, 
a has to be regarded as a function of the free-stream Mach Num- 
ber, M;. We shall assume 
a = A” (6) 

Eq. (6) states that as 4, —~ 0, a0, andas M, > ~, a—1; the 
latter statement appears reasonable, b2cause it is plausible that, 
as M, > @, the effects of density and velocity fluctuations are 
equally important 

In Fig. 1, several curves representing the variation in Cr,/Cr;* 
with M, at R; = 1.1 X 10’, computed on the bisis of Eq. (4) 
with various values of a, are presented. It is noted that a = 0 
corresponds to Van Driest’s (and Wilson’s) curve, a = 0.5 corres- 
ponds to Ferrari's analysis, while a = 1 gives Clemmow’st re- 
sults. We think that the effects of density fluctuations are vary- 
ing with the free-stream Mach Number; thus, a theory on the 
basis of a fixed value of a perhaps may not be adequate for the 
entire speed range that we are interested in (0 < AM, < 10, say). 
A simple variation of a with A/,, as proposed in Eq. (6), is there- 
fore likely to be more useful. Some intensive experimental 
program will be conducted in the higher supersonic range by 
utilizing the Army Ordnance-GALCIT Hypersonic Wind Tun- 
nel to verify and improve the present theory. 
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On Subsonic Conical Methods 


Seymour Lampert 
Ames Aeronautical Laboratory, N.A.C.A., Moffett Field, Calif 
June 18, 1951 


iy THE COMMENTS! by Theodore R. Goodman which appeared 
in the Readers’ Forum of the June, 1951, issue of the JouURNAL 
on the paper? by the author, it was pointed out that, in general, 
the conical solutions to Laplace’s equation do not in themselves 
represent physically possible subsonic flow fields since the bound 
ary condition of zero disturbance at infinity will not be fulfilled 
In reference 2, the example of the uniformly loaded semi 
infinite sector was inserted to illustrate this point, which, how 
have been sufficient clarification. It is clear, 


ever, may not 


however, that by the proper superposition of these conical solu- 


*Cr;, = CF.) M4, =0 


T There is apparently a misprint in Clemmow’s formula for F(A) 
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tions, as was done in reference 2, the boundary conditions at the 
wing and infinity can be satisfied so that a physically possible sub- 
sonic flow may be represented. 
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Small Pitot Tubes with Fast Pressure Response 
Time* 


W. S. Bradfield and G. E. Yale 

University of Minnesota, Rosemount Research Center, 
Rosemount, Minn. 

June 25, 1951 


— NEED FOR A SMALL PITOT TUBE appears fairly often in 
for example, exploration of thin 
The production of 


aerodynamic investigation 
boundary layers is one common application. 
small probes suited to this application has been well understood 
for some years. The experiment of Stanton and Marshall! in 
volved one of the first attempts at this technique. 
the Stanton tube, while small, was unlike the small probes that 
It differed principally in that it utilized as its 
The distance of the 


However, 


are in use today. 
lower surface the actual surface of the plate. 
“tube opening” from the surface of the plate was obtained by 
raising or lowering the top surface of the tube from within the 
plate. Examples of the recent application of small pitot tubes 
of conventional types are found in references 2 and 3. 

Particularly in an intermittent flow wind tunnel, emphasis 
must be placed on fast response as well as on small size. In con 
nection with an experiment?! relating to the determination of ve 
locity profiles within the boundary layer at high speed, the authors 
succeeded in developing a technique for producing a pitot tube 
with an opening 1 mil high which will come to equilibrium in a 
period of time between 5 and 20 sec. 

Analysis of the optimum configuration of a pressure-sensing 
instrument involving a pressure orifice, tubing, and an end organ 


* This work was done as a portion of research sponsored by the Office of 


Air Research, Air Materiel Command 
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Fic. 2. Poor total pressure probe entry 
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Photomicrograph of typical total pressure probe. 
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Fic. 4. A total pressure profile through laminar boundary layer 


at Mach Number 3. 


(such as a manometer or capsule type electrical pickup) has re- 
sulted in a clear specification of conditions to be met.*® Fig. 1 
shows schematically the optimum configuration of a fast response 
pressure-sensing system. The salient point is that the opening 
into the pressure measuring system—that is, the pressure orifice 
itself It is 
clear that, if it is a passageway, the Reynolds Number will be 
small because of its small height. 


must be an orifice; it must not be a passageway. 


The friction coefficient, then, 
will be correspondingly high. In general, the other features of the 
design are simply structural and aerodynamic considerations that 
make certain that the tube between the orifice and the manom 
eter diverges in such a manner that extremely high friction co- 
efficients are avoided. The length of connecting tube is mini 
mized in order that the volumetric capacity of the system shall 
be as small as possible. However, the most important single 
feature is the production of an orifice at the probe opening. At- 
tention to this particular detail made a difference in response time 
of a factor of six to ten in the present investigation. 

The most satisfactory methed of construction was found to be 
the following. Heat gently a piece of !/32-in. brass tubing ap 
proximately 3 in. long by passing a gas torch flame back and forth 
The tube 
will elongate a certain amount and become smaller in diameter, 


across the tube while the tube is placed under tension. 


finally breaking and yielding a tube of considerably less diameter 
The wall of the tube is still thicker than necessary, and the thick 
ness is reduced to approximately 0.002 or 0.003 in. by etching with 
a solution of 1 part copper sulphate and 2 parts water (distilled). 

Obtain a piece of sheet steel 0.001 in. thick by 0.010 in. wide, if 
necessary, by etching a thicker piece with the following solution: 
34 parts sulphuric acid (concentrated), 42 parts ortho-phos 
Slightly flatten 
probe and insert steel approximately !/; in. and flatten probe with 


phoric acid, and 24 parts water (distilled) 
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hammer for about '/i) in. from end. 


with a fine wet stone and withdraw steel insert. 


Polish surface of probe 
View probe 
under a suitable microscope to verify probe dimensions and 
smoothness of surfaces. When the size and smoothness of the 
This 
can be done easily on the bench using a U-tube with the probe 


tip are satisfactory, the response time should be checked. 


connected in series with one leg and a pressure (or vacuum) soure 
with the other leg. If the response time is unreasonable, it js 
most probably a result of too much tip flattening. In this case 
the passageway will appear as in Fig. 2. The error may be cor. 
rected by clipping about '/¢; in. off the end of the tube and re. 
opening the tube by filling or stoning. The clipping procedure 
should be repeated until the response time is satisfactory 

A photomicrograph of a typical total head probe produced by 
the foregoing technique is shown as Fig. 3. The dimensions of 
this probe were checked by the use of a 100-power microscope 
The response time was found to be from 5 to 7 sec. for 21 in, of 
mercury. 

Tubes produced in accordance with the technique outlined in 
the foregoing paragraphs were used to obtain velocity profiles in 
a laminar boundary layer on a cone at Mach Number 3. The 
boundary-layer thickness was as small as 16 mils (0.4 mm 
The tunnel was an intermittent tunnel having a running time at 
this Mach Number of approximately 50 sec. After bench testing 
of the probes used to determine their response time, they were 
applied to the actual determination of the total pressure in the 
By careful and repeated observa- 
tion of the manometer during the run, it was established that a 
definite point of equilibrium was reached during each run and 


boundary layer of the model. 


that the time required to reach equilibrium was not noticeably 
different from that required on the bench. 
being run varied from approximately 0.7 to 50 in. of mercury. A 


The total pressure 


typical total pressure profile at a relatively low Reynolds Number 
is shown in Fig. 4. 

It should be easily possible to make probes of even smaller 
physical dimensions having almost equally good response. The 
metal deposition technique has been suggested as a possibility 
in this regard. 
fibers has been developed to the point where tubes having inside 


For example, the deposition of metal on nylon 


diameters as small as 2 mils have been made.’ 
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Cooling Required To Stabilize the Laminar 
Boundary Layer on a Flat Plate 


E.R. Van Driest 

Aerophysics Laboratory, North American Aviation, Inc., 
Calif 

July 6, 1951 


Downey 


T° HIS REMARKABLE PAPER,! Lees predicted that it would be 


possible, by cooling, to completely stabilize the laminar 
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Fic. 1. Cooling required to stabilize the laminar boundary 


layer on a flat plate regardless of Reynolds Number. 


boundary layer on a flat plate at all Reynolds Nuimbers in super 
sonic flow. In fact, his calculations showed that radiation alone 
would be sufficient to maintain complete stability of the bound- 
ary layer ona plate moving at Mach Number 3 at altitude 50,000 
ft. This was, of course, comforting to engineers engaged in 
high-speed aerodynamics. The assumptions made by Lees were 
that the Prandtl] Number was unity and that the viscosity was 
directly proportional to the temperature. 

Because of the significance of Lees’ conclusion to high-speed 
aerodynamics, the writer has recalculated the amount of cooling 
required to completely stabilize the boundary layer, using the 
more realistic Prandtl Number of 0.75 and Sutherland viscosity- 
The Crocco method was used to obtain the 


temperature law. 
Fig. 1 shows the results of 


properties of the boundary layer. 
the calculations. It seems that considerably more cooling is 
required than predicted by Lees and that, above about Mach 
Number 9, it is impossible to stabilize the boundary layer by 
any amount of cooling. Furthermore, it appears that when the 
layer can be stabilized by cooling, even full radiation (emissivity 
unity) is not sufficient. 

The question arises as to the reason for the radical difference 
between the results of Lees and those of the writer. In the first 
place, the writer’s lower Prandtl Number would require more 
cooling than predicted by Lees, but the writer does not expect 
this to be the major cause of the discrepancy. In the second 
place, it seems that Lees’ assumption of viscosity-temperature 
proportionality should show an even greater effect than obtained 
by the writer because of the greater thickness of the layer com- 
pared to that resulting from viscosities given by the Sutherland 
law. Therefore, it seems that Lees’ original calculations should 
have shown a much greater cooling required to completely sta- 
bilize the laminar boundary layer regardless of Reynolds Num 


ber than was indicated in reference 1. 
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+ 
Secondary Flow in Cascades 


J. M. Stephenson 
Farnborough, England 
July 13, 1951 


Spee RECENT PAPER by Squire and Winter'is the first to de 
scribe the gyroscope effect that occurs when the flow into a 
bend contains vorticity. However, the effect can be derived 
simply from the principle of conservation of circulation (or of 


angular momentum in the case of a gas). 
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Consider the cascade of Fig. 1, deflecting a: stream of incom- 


pressible inviscid fluid. The inlet velocity 7; has a profile of arbi 
trary shape. With the usual assumption that the number of 
blades is infinite, the velocity diagram is as drawn in Fig. 2. At 
every point upstream, there exists a vorticity vector at right 
angles to the stream direction of magnitude — = —0v,/Ox; simi 
larly, downstream is a vector 7 = —0Ov2/Ox. If we assume that 
the spanwise perturbations are small, the profile of the axial (<z 
wise) velocity component vq is unchanged through the cascade. * 
Therefore, if a is independent of the velocity (a common as 
sumption ), we have 
~= —sec a,-OV,/ OX 

n = —sec az2:Ov,/ONX 
To conserve the circulation in the system, an axial vorticity 


vector appears downstream of magnitude 
Esin ay — 7 Sin ag = —(tan a, — tan ae)-dv,/Ox 


which has a component ¢ in the direction of as, such that 
¢ = —(tan a — tan a2)-COS a2: Ov,/OX 

= (tan a; tan a@y)-COS a) COS ae-& (1) 
This vorticity implies a secondary rotation of the fluid as in Fig. 3 

Both the theory and the measurements of reference 1 only 
apply to a single passage in a cascade. But there is also present 
in the downstream flow the circulation shed from the trailing 
edges of each airfoil. This phenomenon is well known; the re- 
sulting vorticity is 


Q = —(1/s)-(O0K/dx) 


where A is the circulation on each airfoil and s is the airfoil spac- 


ing. 
But 
K = sv, (tan a — tan a2) 
Hence, 
Q = —(tan a, — tan a2)-dv,/Ox 


* If this assumption is not made, it can be shown that the subsequent 
analysis is unaltered, except that in Eqs. (1) and (2) tan a: must be replaced 


by C tan a», where 








2 + (tan*o: tan*a2) 
( 
2 — (tan’a: — tan%a2) 
4 
= 
aa!” 
—“ 
wn 
z v, =z 




















Fic. 1. 














700 JOURNAL OF THE 


AERONAUTICAL 

















Fic. 3. 


and 


Q = —(tan a, — tan a2)-cos a-& (2) 


The resultant vorticity in the direction of a2 is thus, from Eqs. 
(1) and (2), 


a+ =(1 


which has always the same sign as ~ but is zero when a2. = 0). 


— COS a») COS a (tan a, — tan az)-é 


Therefore, there is a minimum induced velocity field and second 
ary flow downstream of the cascade if the gas flow at outlet is in 
the axial direction. 

It is perhaps worth mentioning that in the most important ap- 
plication of cascade theory—namely, to the design of axial-flow 
compressors and turbines—the gyroscope effect plays no part. 
For, when annular cascades are closely spaced, only the axial 
velocity component at any point has a peaked profile. There 
fore, the vorticity vector — is not rotated on passing through a 
row, and the only extra component to appear is that shed from 
the blades. It is still true, however, that this component is a 
minimum when the outlet velocity is in the axial direction, but 
for a different reason. 

A small correction is necessary in Fig. 4+ of reference 1, where 
the ordinate scale should be multiplied by 0.9. 


REFERENCE 
1 Squire, H. B., and Winter, K. G., The Secondary Flow in a Cascade of 
Airfoils in a Nonuniform Stream, Journal of the Aeronautical Sciences, Vol. 
18, No. 4, p. 271, April, 1951 


Note on the Hypersonic Similarity Law for an 
Unyawed Cone* 


Lester Lees 
Department of Aeronautical Engineering, Princeton University 
June 4, 1951 


a TO THE HYPERSONIC SIMILARITY LAW," ? at very 
high Mach Numbers the pressure coefficient on the surface 
of a series of unyawed slender bodies of revolution with the same 
thickness distribution along their axes can be written in the form 


* After this note had been prepared, it was brought to my attention by 
S. I. Cheng that this problem had been treated previously by Shen,® who 
utilized Hayes’ suggestion? of considering the flow in a transonic plane as that 
generated by a radially expanding cylinder However, simple algebraic ex- 
pressions were not obtained, and the numerical solutions appear to differ 


somewhat from the present results, except in the limiting case K — « 


SCIENCES—OCTOBER, 1951 


ee. Sn 
(1/2)pi U1? 


é« v G(Ki, &) 
M,? 

Here, K; = M,/(//d) is the similarity parameter, / is body length, 
d is maximum body diameter, J, is free-stream Mach Number, 
and & = x/l, where x is axial distance measured from the nose 
Also, the drag coefficient based on maximum frontal area takes 
the form Cp = F(K)/M;,?. 
hypersonic similarity law for slender cones and ogival bodies 
applies not only for 11 > 1, as postulated in the original deriya- 
tions, but extends also to moderate supersonic speeds (\/ = 3 
The purpose of this brief note is to show that, at least in the case of 
the unyawed cone, the extensive range of applicability of hyper 


Recently, it has been found that the 


sonic similarity is foreshadowed by the analytical expressions for 
the velocity components of the conical flow behind the shock and 
by the boundary condition at the shock surface. In the course of 
the discussion, a relatively simple algebraic expression is derived 
for the drag function F(K) which is applicable when the shock is 
not too far from the cone surface. 


Taylor and Maccoll‘ showed that the continuity and momen- 
tum equations, plus the irrotationality condition, for the flow 
over an unyawed cone lead to the following differential equation 
for the velocity component, 4, along a ray from the cone vertex 


+ I 1 
[? —— yr? 7 - |e -(y — 1) UG = U4) 


: sal Bees ae ek 
(1 U2)U" cot w + yUU"? 4 US cotw =0 (1 
2 2 
Here, U = u/Wmar., Wmar. 1S the limiting adiabatic velocity, w is 


the conical ray angle, and the primes denote differentiation with 


respect tow. Also, V = U’, where V is the velocity component 
normal to a ray, taken positive in the direction of w increasing 
At the cone surface, U = U., U' = 0 and, from Eq. (1), U” = 


—2U., U'"’ = 2U. cot 6-, where 6. is the cone half-angle, and 
U** = —24a,, where 
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2 U.? 
= sot? 
a4 —— T ws - 1 = U2? 


Therefore the Taylor series expansions for the velocity com- 


aie from the cone surface outward are as follows: 
U ‘ l 
= ] (w — 0c)? + = cot Ac(w — Ac)? — ag(w —O-)* +... 
. ; (2) 
V = —Aw — 6.) + cot Oc (w — Oc)? — 4aw — Oe)? +... 
U. 


By substituting the Eqs. (2) in Eq. (3), one obtains 
: 2 ] ] 
y —1| sin? w, + 
Mi(w, — 6) = y¥y-1™, 
y¥+ 1 _ ——<$—<$_— 
sin 2w, a 
But 


2 K ,— 1 
K,? — (K 
: ( l 
sin? w = w* | 1 + 
3 M,? 
Therefore, the relation (4) becomes 
K (: + oe + ) + 
. 3 M2 oe ¥ 
l ¥ l 
K,-—-K#=e ; 
27+ 1 1K, —K 
rere 


Two important conclusions can be drawn from Eq. (6). 

(1) The hypersonic similarity law should apply for the unyawed 
cone so long as w, is less than about 0.4 rad., or w, < 24° (approx.), 
because, when this condition is satisfied, w, and @ enter only as 
the parameters K, and K and terms involving K,?/M? and 
(K, — K)?/M,? make only small contributions to the final result 
Thus, hypersonic similarity should begin to apply for the slender- 
est cones (@. ~ 0°) at M, =~ 2.5 —3, and the range of applicabil- 
ity increases rapidly as 4M increases, up to a value of 0 of about 
20° for air (y = 1.4) when 1, — © (see curve w, = 24.4°, Fig. 2) 

In Tsien’s original derivation,! the square of the perturbation 
velocity in the flight direction must be small compared with the 
square of the velocity normal to the flight direction in order for 
hypersonic similarity to exist. From Eq. (2) and the shock rela 
tions, it can be seen that this condition is satisfied when w? is 
sufficiently small—i.e., w < 24° (approx.), and 2.5 + 3.0 < V/< 
Also, under these conditions, the /ocal sonic velocity is small com 
pared with the local velocity parallel to the flight direction 


(a?/U? ~ w*), as required in Hayes’ analysis. 

(2) When, in addition, (Ks K)/K is small compared with 
unity—i.e., when the conical shock is not too far from the cone 
surface—Eq. (6) is approximated as follows: 


K K 


or 
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and the succeeding terms do not bring any new factors into this 


discussion. 


Now the boundary condition at the conical shock can be put in 
the form (see, for example, Kopal®) 


_ oe ~ = tal ws (3) 
y¥+1 Uz) 
I 4 
(ws 6.)? + , cot 6 ws 6.)* ay(ws — O-)* + 
i - (4) 
1 - cot O- (ae — Oc) + Zas(w.s 6.)? + 
K, — K)*}1 en 
‘ M,? 
x 2 ' (K, — K)? | a 
— K)? + ~ t...1+.. 
y¥—1 M,? 5) 
) ete., where K, = Mio, 
(: 4282 ry ) 1 ' (K, — K)? 
1 3 M2 ) Ky — 
(K, — K)? 6) 
(« . + | a - 
- 4. +O<. ? 
K )X.’ + \ 
¥ l 
+1 (vy +1) 2 
ae K ra? ) K? + (7) 
¥+3 7¥+3 7¥+3838 


This approximation, which amounts to utilizing only the first 
term in the expansions of Eq. (2), should be valid for K = 1.0 
(approx.) within the range of values of 1/7, and @. for which hyper 
sonic similarity applies. In other words, this approximation is 
valid for M, > 4 and 1/M < @, 
upper limit of the ‘‘zone’’ of hypersonic similarity. 
relation between conical shock angle and cone angle prescribed by 
Eq. (7) is compared with the actual values given by the Taylor- 
For K > 0.5, the dif- 


ference between exact values of conical shock angle and the pre 


6, where @ corresponds to the 
In Fig. 1 the 


Maccoll theory as computed by Kopal.® 
diction of Eq. (7) is negligible. 


When the approximate relation given by Eq. (7) is valid, a 
correspondingly simple expression is obtained for the pressure co- 


efficient on the surface of an unyawed cone. Now, 
Pe/ Pr = (ps/Pri) (Pe/ Ps) (8) 
where 
p 27 _ : , 
= VW? sin? w — — (8a) 
pi i 73 
and 
Pe/ Ps = (c*/a 1) (Sb) 
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Fic. 2. Comparison between approximate drag formula for 
unyawed cone and exact values as computed by Kopal. The 
curve w, = 24.4° shows the approximate limits of applicability 
of the hypersonic similarity law. 


Here, p- and a- are the static pressure and sonic velocity, re- 
spectively, on the cone surface, and p, and a, are the static pres- 
sure and sonic velocity just behind the conical shock. From the 


energy equations, 


a* + — 1 a a? u- \? v, \? 
=]¢ - - 1 — + (9) 
a? 2 a;* a," u a 


so that, with v./u. ~ —2(w, — 0-) and u,/ue ~ 1 — (ws 6.)?, 
a-* ay) ° ae 
1 + (¥ Li Um K)? 
a,° a, 
(10) 
p a," . 
<= 1+ y7— (A, — K) 
Ds a,* 


From Eqs. (8a), (8b), and (10), 


F(K) =~ (2 ~ ) CoM? = 
Y \b: 


4 y7¥+ 1 
(K,2—1)+2(K, — Kk)? - : 
ly —1+ (2/K.2)§ 


y+1 


(11) 
and K, is related to K by Eq. (7). The drag function F(K) is 
plotted in Fig. 1. Again, this approximate formula is expected 
to be valid for K 2 1, or for \/, > 4 (approx.), and 1/M, < @. < @. 
For comparison, the drag function Cpl/J,;? = 2K? In (2/K) ob- 


tained from the ‘“‘slender-body”’ approximation to the linearized 


supersonic theory® is also shown, for K <1. In this approxima 
tion, the shock or Mach cone must be far from the body, while in 
the present approximation the shock cone is close to the body 

In Fig. 2, the drag coefficient for an unyawed cone given by Eq 
(11) is compared with the ‘‘exact’”’ values as computed by Kopal 
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The approximate limits of applicability of hypersonic similarity 
are indicated by the dashed curve labeled w,; = 24.4°. For K >| 
the present approximation gives values of Cp within 7 per cent of 
the exact values. An interesting result is that, for the cones of 5 
and 10° half-angle, the von Karman approximation and Eq. (11 
taken together give an approximate analytical representation of 
the drag coefficient over practically the entire range of supersoni: 
velocities. At very high Mach Numbers, Cp — 2.08 sin? 6, for 
y = 1.405, which is close to the Newtonian value. 

Recently, hypersonic similarity has been shown to apply also 
for the slender yawed cone at moderate supersonic speeds,’ |) 
principle, there would seem to be no difficulty in extending the 
present considerations to this case, and perhaps also to ogiyal 


bodies. 
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A Limitation on Shock Position 


E. W. Graham 
Consultant in Aerodynamics, Douglas Aircraft Company, Inc 
Santa Monica, Calif. 


May 18, 1951 


HE FOLLOWING DEVELOPMENT APPLIES in the steady adiabatic 
gpm of a compressible fluid past a body if the flow upstream is 
uniform. The Mach Number and body fineness ratio are not re- 
stricted, but the results should be of greatest interest for blunt 
bodies and Mach Numbers slightly greater than unity. 

The maximum mass of fluid which can flow through a unit 
cross-sectional area in unit time is p*a* (where p* is the density 
and a* is the fluid velocity existing when the speed of sound is 
reached locally). Fora steady flow p,*a* (the value of p*a* in the 
uniform free stream ahead of a body) cannot be exceeded at any 
other point in the flow. The value of a* is constant, and p* be- 
haves in the same fashion as stagnation pressure and can only de- 
crease as the entropy is increased by viscous effects and shocks 

Since an upper limit for the mass flow of fluid is known, it is 
possible to establish an inner limit on shock position based on con- 
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tinuity. For the two-dimensional body shown in Fig. 1, the 
maximum possible flow across the line AB is p,*a*/. The flow 
across the line AC is pim(1 cos g + ys) and by continuity 


p\*a wi | . piu(L cos ¢ + VB) 
Therefore, 
l= 1 ¥B [(pi*a */piui) — COs ¢]} 


By finding / for different points on the body and different angles 
¢, an inner limit for the shock position is established 

A similar analysis can be made for bodies of revolution, assum- 
ing that the line AB lies in a conical surface connecting the circle 
of body points at one radius with the circle of shock points at 
another radius. However, a more positive statement about 
shock position can be made by considering the flow through the 
surface of minimum area connecting the two circles. In the two- 
dimensional case, AB was a line element in a plane, and the plane 
was the surface of minimum area connecting the specified end 
points. For the body of revolution the corresponding conical sur- 
face is not the surface of minimum area connecting the circle of 
shock points with the circle of body points. The minimum area 
surface is obtained by rotating the catenary 


cosh {(x C.)/C,} 


y= C, 


about the axis of the body, C; and Ce being determined by the 
end points 

This makes the geometry considerably more complicated fort 
the rotationally symmetric case. The distinction between the 
conic surface and the surface of minimum area is probably most 
important as the Mach Number approaches unity. 

The actual shock position lies outside of the limit calculated 
above for two reasons: (1) The actual p*a* behind the shock is 
less than the limiting value p;*a* because of entropy increases; 


(2) These actual values of p*a* are, in general, not attained at 
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every point on the surface of minimum area. It can be concluded 
from (1) that the limiting shock position could be far from the 
actual position at high Mach Numbers where the entropy in- 
creases are large. However, at Mach Numbers slightly greater 
than unity, the agreement might be good, depending upon the 
shape of the sonic line. 

In Fig. 2 the significant portion of the limiting shock position is 
shown for two-dimensional bodies and bodies of revolution at 


several Mach Numbers. 


Comment on “Vibration of Rectangular 
Plates’”’ 


Charles Massonnet 
Université de Liége, Belgium 
June 8, 1951 


I NOTE IN THE READERS’ Forum of February, 1951, the paper by 
Harold Lurie entitled ‘‘Vibrations of Rectangular Plates.’ 
Nearly the same demonstration is given in my book entitled Les 
relations entre les modes normaux de vibration et la stabilité des sys- 
témes élastiques, published in 1940. This book formed parts | and 
2 of the Bulletin des Laboratoires d’essais des Constructions du 


Génie Civil de I’ Université de Liége. (See pages 273-287.) 


Author’s Reply 


Harold Lurie 
The RAND Corporation 
June 15, 1951 


~ A READERS’ ForRUM NOTE (February, 1951), I derived an 
I expression relating the natural frequency of vibration of a 
rectangular plate and the plate buckling factor &. It has subse- 
quently been brought to my attention that Massonnet had pre- 
viously derived an expression relating the natural frequency to 
the buckling load. My note does, however, present the results in 
terms of the N.A.C.A. buckling factors, enabling actual numeri- 
cal calculations of frequency to be made by direct use of the 


N.A.C.A. tables. 


Shear Flow in a Thin-Skin Tapered Beam 


H. W. Sibert 


Professor in Aeronautical Engineering, University of Colorado 
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SUMMAR\ 


It is proved that a simple formula involving the shear force carried by the 


webs. which is sometimes used for the shear flow in a thin-skin tapered beam, 
is not valid in general but will give correct results under certain specified 


conditions 
SYMBOLS 


1 — moment of inertia of the cross section under investiga- 
tion about its x axis 

WV = external bending moment about x axis of the cross sec- 
tion (+ when compression at + values of 





JOURNAL OF THE 


total shear force on the cross section (assumed to be 
parallel to z axis and + when in +2 direction) 

portion of V carried by the stiffeners of the cross sec- 
tion (+ when in +: direction) 

portion of V carried by the webs of the cross section 
(+ when in +2 direction) 

area of ith stiffener of the cross section 

normal force on A; (+ when compression) 

distance from x axis of cross section to centroid of 
area A; 

dz;/dy 

difference in shear flow in the two webs of the cross 
section on either side of area A; 

number of stiffeners intersecting the cross section 

Norte: Subscript z will be used to indicate a particular stiffener. 
Subscript & will be used instead of 7 when a summation is to be 
taken for all the stiffeners of a cross section. 

ASSUMPTIONS 
THE FOLLOWING ANALYSIS it will be assumed that: 

(1) The x axis is a principal axis of the cross section under 
investigation, and the y axis is perpendicular to this cross sec- 
tion. 

(2) The skin of the beam is thin enough to assume that only 
the stiffeners resist the bending. 

(3) The cross-sectional area 
throughout its length, 

(4) Each stiffener is sufficiently far from the x axis for 2;24, 
to be its moment of inertia about the x axis to a close approxima- 


of each stiffener is constant 


tion. 
(5) The curvature of each stiffener is small enough to assume 


that the axis of the stiffener is a straight line. 
Although an analysis can be made when the curvature of the 
stiffeners cannot be neglected, it would be difficult to draw any 


general conclusions from the resulting formulas. 


FORMULAS 


Consider a portion of the beam of longitudinal length dy. The 
difference between the normal forces on the two ends of the 7th 
stiffener is dP;, and the difference in transverse shear flow in the 
two webs of the cross section on either side of area A; is gi = 
dP;/dy. Since the bending stress at area A; is Mz,/I, P; 
A;M:z;/I. SincedM/dy = Vand dz;/dy = sj, 


qi = (VaiAi/T) + (Ms:Ai/T) — (Msi Ai/I2)dI/dy (A) 


Instead of Eq. (1), the following simple formula for g; is sometimes 
used: 


qi = VwriA¢ I (2) 


with V, the shear force carried by the webs of the cross section. 
The validity of Eq. (2) will now be investigated. 

Since the correct value of g; is given by Eq. (1), a necessary 
and sufficient condition for Eq. (2) to give the correct value of g; 
for each stiffener can be found by equating the right sides of 
Eqs. (1) and (2). Since V — V, V, by definition, the result 
ing equation when multiplied by /?/ \/z;A; is 

(V.I/M) + (Is;/2:) — (dI/dy) = 0 (3) 


Since the normal force Py, on A, is parallel to the y axis, the 
z-component of the resultant force on A, is Py dzx/dy = Psy 
Msy2xAx/I. Then V,, the shear force carried by all the stiffen 


ers, is 


n 


(M/T) SS sezeAs 


k 1 


n 
Vy, = > P ysy 
k= 1 
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Since it was assumed that only the stiffeners resist the bending 
and since dz,/dy = s,, it follows from Assumption (4) that 
n 


n 
Zz 2°Ak, dI/dy = 2 Zz. SEERA & (5) 


k= 1 k=1 


[= 


The substitution in Eq. (3) of the above values of V,, 7, and dJ ‘dy 
gives 
n n 
(S; ‘2;) > oK7A E —_ y SERA = () (6) 
k=1 k=1 


= @,° . ‘ . 
as a necessary and sufficient condition for Eq. (2) to give the 


correct value of g; in the 7th stiffener. 

Now consider the lines formed by projecting the axes of the 
stiffeners upon the yz plane. Since all of these projection lines 
-annot be parallel in a tapered beam, at least two of them must 
meet at a point (call it 0). If » is measured from point 0, it 
follows from Assumption (5) that the equations of these projec- 
tion lines are: 


Sk = be + sey (7) 


in which },; is a constant for each projection line and 6, = 0 for 
the projection lines that meet at point 0. 

Since Eq. (6) must be satisfied for each stiffener if Eq. (2) is 
to give the correct value of g; at each stiffener, a necessary con- 
dition for the validity of Eq. (2) can be based upon the projection 
Then Eq. (6) for 


lines that meet at point 0, for which 2; = s;y. 
these particular projection lines becomes 


n n 


(1/y) y> (bi + sey)*Ak — > Su bK + SEY)Ak = 
k 1 k 1 


which reduces to 


n n 
(1 y) p> by? An + >. bist: At = 0 
k 1 kj= 1 


and A, are constant for each stiffener and, 
Since Eq. (9) must be valid for all 


Note that bz, sz, 
hence, are independent of y. 
cross sections of the beam (that is, for all values of y for which 
there is a cross section of the beam), each summation in Eq. (9) 


must separately equal zero. Hence, it is necessary that 


n 

> oxtAn = (10) 

k 1 
if Eq. (2) is to give the correct value of g; at each stiffener. Since 
neither 5,2 nor A; can be negative, Eq. (10) will not be satisfied 
if each b; is not zero. Hence, each b; = 0 is a necessary condi- 
tion for Eq. (2) to give the cerrect value of g; at each stiffener. 
Conversely, if each b; = 0, Eq. (9) reduces to zero. It follows, 
then, that: 

A necessary and sufficient condition for Eq. (2) to be valid for a 
thin-skin tapered beam composed of essentially straight stiffeners 
is that the lines formed by projecting the axes of the stiffeners upon 
the yz plane shall all meet at one point. 

It should be noted that this condition rules out the use of Eq. 
(2) in the important case of a wing with straight stiffeners which 
is tapered in both plan form and thickness 

Note that from Eqs. (4) and (5), d//dy = 2V,1I/M. 
stitution of this value of d//dy in Eq. (1) gives 


+ Ms;]A;/I 


The sub- 


di [((V — 2V,)z; (11) 


If the condition mentioned above in italics is not satisfied, Eq. 
(11) should be used instead of Eq. (2). Note that Eq. (11) is 
only slightly more difficult to apply than Eq. (2), because in us- 
ing Eq. (2) it is necessary to calculate V, before Vy (which equals 


V — V,) can be found. 











